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In this short note we give some comments tacent paper-by М; . Bokhari, 
o NS OF CEE A r 57 
A. Qadir, and Н. Al-Attas, On Gauss-type quadrati rules; Ruiter ипе. Anal. 


——— 


Optime. 31 (2010), 1120-1134. Their polynomials аге a spectat case of the Jacobi 
polynomials on (0,1). In addition we construct orthogonal polynomials 
T, (x),n =0,1,..., and the corresponding Gaussian quadrature rules with respect 


to the linear B-spline (as a weight function) and give some numerical examples in 
order to illustrate an application of such quadratures. 
2010 Mathematics Subject Classification : 66D30, 65D32, 33C45 
Key Words and Phrases : Orthogonal polynomials, Jacobi polynomials, weight 
function, three-term recurrence relation, Gaussian quadrature rule, linear B- 
spline, nodes, Christoffel numbers. 

1. Introduction. Recently Bokhari, Qadir, and Al-Attas [1] considered 


Gauss-type quadrature rules based on polynomials p, (4) orthogonal on (0,1) with 
respect to the linear weight function œ(t):=1-t. They discussed a development of 


р,(1) via Gaussian hypergeometric differential equation, narrated some of its 


properties, derived the three-term recurrence relation for the monic polynomials 


»ь(0-(1:22:1 11, кра төд. (1.1) 


where p,(t)=1 and р, (@)=0, and considered several numerical examples of such 


“kind of quadratures. 
In this short note we show that these polynomials p,(¢) are a special case 


of the well-known Jacobi polynomials on (0,1). In Section 3 we construct orthogonal 


polynomials л, (x),n=0,1,..., and the corresponding Gaussian quadrature rules 


with respect to the linear B-spline (as а weight function). Finally, in Section 4 we 
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TT 
_ 
give some numerical examples to illustrate and application of such quadratures. 0 
7 р ү | | 
2. Jacobi Polynomials on (0,1). Let w(x) = (1-х) (1+) ,0В> 1, and | 
rm (9) be a sequence of the corresponding monic Jacobi polynomials, which : 
satisfy the three-term recurrence relation (cf. (8, pp. 131-140]) 
PEA (x)= (xà) BP (x) B, P (om 40.» 
where Бі (x) - 1, P9 (x) - 0; end : 
нг 
5 (Qn+a+f)(2n+a+B+2) 
(є 
‚ __An(nta)(n+B)(n+e+B) (n21) 
(2n+a+f) |(2n+0+B) -1| 
except when ; then . і 
By a change of variables x=2t-1 we get the monic orthogonal polynomials ( 
pe? (= 2" Pi» (2; -1) orthogonal on (0,1) with respect to the weight function 
o(t):- (1-1) t^,a, B» -1. The coefficients in their three-term recurrence relation 
( 
рї (t)= (2-a, ) Pr ()-8,Р (0), л=0,1,.., 23) 
аге 
1, .. (2n+a+B+1) -(1+07-p? 
E or "o, (I 
2|(2n+a+B+1) -1| 
| I 
| m nin a )(n- B)Un 4 oc 
| Egon 
R 
| For о =] and = 0, this relation (2.1) gives orthogonal polynomials р, (t) а 
discussed in [1], which satisfy the relation (1.1). : 
| A t] 
| N 
| ow, we list the parameters о (n 2 0) and B, (n 21) for some other special 
| cases: 
| (à  a-0p-i. S 
| І 2(n+1) _ n(n+1) E 
Te € 


{ a = 
| iTO JA E 
| п+1) (2n +3)" 4 d 
| | СС-0. GuruKul Kangri ооо (9134) An eGangotri Initiative e 


| (b) a=1,B=1. 
| . 
| C cut W2) а 
2 4(2n -1)(2n +3)’ 
(c) в = 2,p 20. 
п? -3n«1 л (n2) 


CUN Ээ гд 


л 


2(n*1)(n-2) ^ 4(n+1) (2451) +3), 
(d) а= 2,В =1. 


2(n+1)(n +3) п(п+ 3) 
Qa, Sa ыйкы соро 
(2n«3)(2n-5) " 4(2n«3) 
(e) 0 22,p-2. 
sd n(n+4) ; 
Sa = gha “Тол к CE 
(£f) а = 3, B = 0. 
2n* +8п + 3 п (n3). 
Qa, Se = nN EN 
(2n 4 3)(2n 4 5) 4(n+1)(n+2)(2n+3) 


(g) 00-43: haus 
"m (n+1)(n+4) a _ n(n+1)(n+3)(n+4) 
| " 2(n«2)(n-3) " 4(п-2) (2n+8)(2n+5) 


(0222-0552 
| 1 1 1 
| ш- = — |) = — >2 5 С. 
| e, 2 ‚В, 8 ‚В, 16 (л ) еїс ; 


Remark 2.1. The relation (2.1) can be obtained taking w(x)-|x| (1-27) , with 


and the corresponding generalized Gegenbauer polynomials , where , which were. 
introduced by Lascénov [7] (see, also, [2, pp. 155-156] and (8, pp. 147-148]). Their 


three-term recurrence relation is 


a a, а) = 
WEP (a) = Wi (х) - B,WEP (x), Wo? (x) = LW C9 (x) =0, 


with recursion coefficients 
Pun n(n+a) вг eG Gee 
3 ^ (2nco*B)(2noa--B*1)' ni ^ (2n«o*B-1)(2n*o.* D) 


xs except a-rg- -1; then. (ST) Кайап Collection, Haridwar. An eGangotri Initiative 


"mmm — =—3hS—h 


: n (-1,1), using Theorems 2.2.11 and 2.2.12 


: i ion is even 0 i 
EN e lynomials orthogonal with respect to the 


from [8, pP- 102-103], we get (2.1) for po 


щй = (1-0, with o =B; =(B+1)/(a+ß+2), Gi. Jb, EB, 


2m: l9 


weight c(t) = 


p, = Byq1Bans™ 21. 


3. A Gaussian quadrature formula. Sometimes in applications it could 


1 = .., and the 
be of some interest to construct orthogonal polynomials п, (x), n = 0,1,..., and the 


corresponding Gaussian quadrature rules with respect to the linear B-spline (as a 


weight function) 


1+х, -1<х<0 


ш(х)= В.(х)=11-х 0<х<1 
0, otherwise, 


[,f() Bi (х)ах = |70) (2+8) а = > Ad (а) +R,(f), (3.1) 


where R,(p)=0 for all polynomials p of degree at most 2n-1. 
This weight function is an even extension of o(/) =1-/ from (0,1) to (-1,1). 
The coefficients in the three-term recurrence relation for the corresponding monic 
orthogonal polynomials л, (х), 
Tasi (x) = XT, (x) ij 74777 (x),n = 352 ос 
are 
B, = 1,8, =1/6,B, =7/30,р, -57/245,6, =683/2793,8, = 207725 /856482, 
Ве = 286749501/1159331030, }, = 286268318986 / 1164429355245, 


Bs = 272609711230510/1097298927604497, etc. 
For example, : 


(x) -1,л, (x) =x,m (x) =x" me) =x -Z n, (x) = хі 814719 


Б 49 ' 490" 
s 90x? 109x 4 2 г 
T(x) =x ^m UL es Ea AT . 5935 - ai 
98 15026 7513  631092' ^^ 
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5 
In numerical construction we use our Mathematica Package “Orthogonal 
Polynomials” [3]. In order to construct quadrature rules up to n points we need 


22 k 
the moment H, = ЇЇ: В, (x)dx,k =0,1,...,2n-1, which are in our case ШОЛУ, 


1+(-1)° 


= —_-—_—_ k>0. 
Нь (k+1)(k+2)”” 


By the mentioned Package, in this case, we can obtain the recursion coefficients, 
a, (=0 in this symmetric case) and , in a symbolic form for a reasonable n: 


(alpha, beta}=aChebyshevAlgorithm[moments, Algorithm->Symbolic]; 
taking, for example, the first 800 moments, 

moments= Table [(1+(-1) *k)/((1+k)(2+k)), {k,0,800}]; 

This enables us to construct parameters in Gaussian quadrature 


Q, (f) = МАГ ЫН 
ї 1 
up to n <400 nodes. For example, for n=100, with Precision ->40, the statement 
is the following 


{n100,w100}=aGaussianNodesWeigths[100,alpha, beta, 
WorkingPrecision->50, Precision->40]; 


where n100 and w100 are sequences of nodes x” ,k — 1,..., n, and Christoffel numbers 


А(" k=1,...,n, respectively. The last command implements the wellknown Golub- 


Welsch algorithm [6]. 

The corresponding software in Matlab was given by Gautschi [4], The first 
FORTRAN package of routines ORTHPOL was also developed by Walter Gautschi 
151 in 1994. 

4. Numerical results. In this section we consider three integrals 


I= [5 (х)В, (x)dx,k = 1,2,3, 


for the functions f, (х) = cos(nx), 


f 1 1 М) 
раа а 
СЭР EX тээ aoi 
(z+ i5) “100 Ё a 25 
and then apply the quadrature formula (8.1) fo 
is Smooth, the second is quasi-singular, and 


i in Fig. 4. th 
The graphics of integrands 5.(5)8, (x), k=12, are displayed in Fig. 4.1, and the 
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r their calculation. The first function 
the last is an oscillatory function. 


ee 


6 


graphic of the oscillatory function f,(x)B,(x) is presented in Fig. 4.2 


Using our Mathematica Package “Orthogonal Polynomials” [5] we 
constructed n-point Gaussian rules for n=5; 10(10)50; 100(50)400. Then, we applied 


these rules to integrals I,,k =1,2,3, and compared the obtained results with the 


exact values { 
l= 4/2 = 0.405284734569351085775517852838910555617435... 


І, = 92.3859074908351131640822746542985'112503464... 
‚= 1.15470053837925152901829756100391491 129520... 


-02 


Figure 4.1. Graphics of function x > B, (x)f, (x) (left) and x — B, (x) f, (x) (right) 


-05 
К! Fi T ^ 1 5 E 
| gure 4.2. Graphics of the function x B, (x)f, (x) 


In Table 4.1 i i 
| we give the relative error in Gaussian approximations 


q d 
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мэ? 
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Th 


Sere NY 


REL, ( h=1,2,3. 


-9.(5)-4, 
1 


k 


Numbers in parentheses indicate decimal exponents. 
Table 4.1. Relative errors REL, (п) in Gaussian approximations Q, (f) for k=1,2,3 


[pese pesas 


SL 


2. 718 5) 8. 356 1) 1.39(- 

4.01(-15) 1.25(-1) i 

1.01(-40) 1.81(-2) 1.19(-2) 
2.50(-3) 1.39(-2) 
3.26(-4) 3.69(-4) 
4.05(-5) 8.82(-4) 
4.22(-10) 2.68(-6) 
2.28(-14) 1.05(-6) 
9.20(-19) 3.63(-9) 
4.88(-24) 3.79(-10) 
6.49(-28) 1.49(-12) 
1.99(-32) 1.81(-14) 


3.40(-38) 7.64(-16) 
As we can see in the case of a smooth function f, the convergence is very fast. 


The relative error is about 10-49 for n=20 points. As we expected the convergence in 
the last case is very slow. We must use 400 points in the quadrature rule in order to 
obtain the so-called double precision result (with about 16 decimal digits). 


ACKNOWLEDGEMENTS | 
The work was supported in part by the Serbian Ministry of Education, 


Science and Technological Development (No.#O01 174015). 


[1] 


(2] 


[3] 


[4] 


REFERENCES 
M.A. Bokhari, A. Qadir and H. Al-Attas, On Gauss-type quadrature rules. Numer. Funct. 
Anal. Optim., 31 (2010), 1120-1134. 
T:S. Chihara, An Introduction to Orthogonal Polyno 
1978. 
А.8. Cvetkovic and G.V. Milovan 
Polynomials”. Facta Univ. Ser. Math. Inform., 19 (2 
W. Gautschi, Orthogonal polynomials (in Matlab). J. 


215-234. - 
W. Gautschi, Algorithm 726: ORTHPOL- A package of routines for generating e 
polynomials and Gauss-type quadrature rules. ACM Trans. Math. Software, 20 (1994), 


CC-0. Gurukul Kangri Collection, Haridwar. An eGangotri Initiative 


mials. Gordon an Breach, New York, 


ovié, The Mathematica Package “Orthogonal 


004), 17-36. 
Comput. Appl. Math., 178 (2005), 


| 


"n 


LT 


m 


rl 
_ 


8 


21-62. | 
[6] G. Golub and J.H. Welsch, Calculation of Gauss quadrature rules. Math. Comp., 23 | 
! 


(1969), 221-230. 
[7] K.V. LaSéenov, On a class of orthogonal polyno 
Inst., 89 (1953), 167-189 (Russian). 
G. Mastroianni and G.V. Milovanovié Interpolation Processes-Basic Theory and 


Applications. Springer-Verlag, Berlin-Heidelberg, 2008. 


mials. Ucen. Zap. Leningrad. Gos. Ped. 


[8] 


р ү 
{ СС-0. Gurukul Kangri Collection, Haridwar. Ап eGangotri Initiative 


" 


Га. m 


f 


—_— 0-6 


| 
! 
1 
| 
! 
А 


| 


~ 


Jnanabha, Vol. 44, 2014, 9-24 
HARDY AND COWLING-PRICE THEOREMS FOR A CHEREDNIK TYPE 
OPERATOR ON THE REAL LINE 
By 
Mohammad Ali Mourou 
Department of Mathematics, College of Sciences for Girls University of 
Dammam, PO. 1982, Damman 31441, Kingdom of Saudi Arabia 
E-mail : mohamed_ali.mourou@yahoo.fr 


(Received. : December 16, 2013; Revised : February 9, 2014) 


ABSTRACT 

This paper is aimed to establish Hardy and Cowling-Price type theorems 
for the Fourier transform tied to a generalized Cheredink operator on the real 
line. 
2010 Mathematies Subject Classification : Primary 33C45, 43A15, 43A32; 
Secondary 44A15. 
Key Words and Phrases : Differential-difference operator; Generalized Fourier 
transform; Hardy and Cowling-Price theorems. 

1. Introduction. In his 1933 paper [8]. Hardy obtained the following 
famous theorem : 


Theorem 1.1. Let 1< p,q < with at least one of them finite. Let f bea measurable 
function on R such that 

e* f e L (R) and e* 7 (f)e L (R), (1) 
for some positive constants a and b. Then 
• If ab>1/4, we have f=0 almost everywhere 


* If ab «1/4, there are infinitely many nonzero functions satisfying (1). 


Above F stands for the ordinary Fourier transform on R given by 


F, (0)09 = J, f 9e "dx. 


Later, Cowling and Price [4] obtained the following L^ version of Theorem 


1.1: 

Theorem 1.2. Let f be a measurable function on R such that 
| e" f e L (R) and e™ z (f) eL (R), 

for some positive constants a and b. Then 

• If ab» 1/4, we have (-0 almost everywhere. 


x 2 o =C en Яс ec: 
ч аб- 1/ 4, the function О Ehe form. fe war. An eGangotri Initiative 


(2) 


Бэй. 


d 
2) infinitely many nonzero functions satisfying (2). | disc 
here are шїї | 
$ Pus Ёс оо of Theorems 1.1 and 1.2 to new contexts have been | са 
i à For instance, these theorems have been obtained in [2] for semi-simple / 
nr : in (51 for the motion group and in [15] for Chébli-Triméche hypergroups. | T 
25 The intention of this paper is to establish analogues of Theorems 1.1 and diff 
1.2 when in (1) and (2) the usual Fourier transform 7, is substituted by a generalized 
Fourier transform F, on R associated with a the first-order singular differential- adn 
difference operator: 
df AE fe (- 
у= ANS) FG) 19) or (2), 
Nf) dx A(x) 2 
where 
A()- ki" В(8), о> -1/2, | 
яр | whe 
B being a positive C^ even function on R, and p=0. In addition we suppose that | 
(i) ^ Ais increasing on [0,-| and lim, „ A(x) ==, 
Ab 
(ii) АЈА is decreasing on 10,» [ and lim, ,. A’ (x)/ А(х) = 2р; 
(ii) There exists a constant $> 0 such that the function е" (A’ (x)/ A(x) -2p)is 
bounded for large x>0 together with its derivatives. 
Notice that the differential-difference operator : Мо 


2 


which is referred to as the Jacobi-Cherednik operator (se 
as Awith 


Daaf (a) =F [0ана ваъ | 82-9 арал) (с) | Re 


e [7] is of the ише {уре | by 


A(x) = (sinh x)" (со8вх)"а»р»-1/2: 
р=0+В+1; 822. 


The one-dimensional Cherednik operator (see [3]) is a particular case of 


D, s. Such operators have been used by Heckmann and Opdam to develop a theory 


| generalizing harmonic analysis on symmetric spaces (see [9,12].) For recent 
| important results in this direction we refer to [13,16,17]. (ii) 
The author [11] has initiated a quite new commutative harmonic analysis | 
E on Pe real line related to the differential-difference operator A in which several 
analytic structures on R were generalized 

4 CC-0. Gurukul VETE Ta dien ed. The tools ас tually required for the 
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discussion in the present paper, are essentially the Fourier transform and the 
Gaussian kernel tied to A. 


2. Preliminaries. In [11] we have shown that for each AEC 


the 
differential-difference equation 
Ли = и, и(0)=1, 
admits a unique C” solution on R, denoted Ф, given by 
З Дог == : Ёс 
Ф, (x) TEL (x) ifA#-ip, 
$7] ш 3 
1+ “A(t)dt Pas (3) 
A(x) ! ( ) 1 1р, 
where ф, denotes the solution of the differential equation 
Au = -(A* +р?)и, u(0)=1, u'(0) 2 1, Ж 


A being the second-order singular differential operator defined by 


1а а 
села 5 
ER (24). (5) 
Moreover, Ф, (x) is entire in А. 


Remark 2.1. For A(x) =(sinh|x|) ^ (cosha)^ ^ a2 8 > 1/2, the differential 


operator A reduces to the so-called Jacobi operator. The eigenfunction @, is given 


Ф, (x)= ,R, (Bita OBI ата) 
2 2 
Where „F, is the Gauss hypergeometric function [10]. 
From [14] recall the following estimates. 
Lemma 2.1. (i) For every хє, 


(6) 
el <@,(x)<1. 
00) There is a constant C>0 such that 
а" n (шах -Р)Х| (7) 
Ice secret 
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for all xeR, AEC and n=0,1.... 
Analogous estimates for Ф, (x) 

Proposition 2.1. There is a constant C>0 such that 


аш р |, 


A” o, (x) sci poda lat 


ал 


for all xe R, Ae C and п= 0,1... 
Proof. By (3), 


are provided by the next statement. 


(8) 


Ф E pM, 2) 
T o Ge ree е Mle 


As by (4), 
d A +p” үн | 
AA 090) |, e @А(@а, (9) | 
we obtain 
ал i а | sgn(x) çh d" 
Tu = [A ЛА( 
ааа) 38 А(х) Ji E (iA +р)о, (0) JA( 


The result follows now from (7) and Leibniz formula. 
Note 2.1. For a function f on R, write f,(x)-(f(x)*f(-x)/2 and 
f(x)» (f (x)- f (-x))/2 respectively for its even and old parts. We denote by 


e S(R) the space of (= function f on R which are rapidly decreasing together 
with their derivatives, i.e., such that for all m,n=0,1,.. 


Pan (f)=sup(1+22)" A 


art (x) 


The topology of S(R) is defined by the semi-norms P._.m.n= 0,1,... 


myn? Э 


« со 


5,(R) (resp. 5, (R)) the subspace of S (К) consisting of even (resp. odd) functions. 
• S’ (R) the space of (* functions fon R such that for all m,n=0,1 


OO <<< M M— M — NNNM 


2 
The topology of & UR) is efin 


ukul efined Ву, dhesemie MABE ahi miva = 0,1... 


— ÁN : 
| 13 Р 
| | 


. S? (R) (resp. S; (R)) the subspace of S? (R) consisting of even (resp. odd) 


functions. 
e т the map defined by Zh(x)- |А), х ER. 


Remark 2.2. (i) By (6) we see that S’ (R)c S(R). 
(ii) Itis easily checked that S” (R) is invariant under the differential-difference 


operator A. 
(iii) Due to our assumptions on the function A there is a positive constant Ё 


such that 
A (x) ~ ke” Pi ав ә е. (10) 
The following technical lemma will be useful. 


! Lemma 2.2. The map Z is a topological isomorphism from S; (R) onto 4 (R). 


| | 
9 | Proof. It is sufficient to show that Z maps continously S, (R) into 5, (E). Let I 
7 ! 
feS?(R). Clearly Zf is a C^ even function on R. For n= 52. | 
@©„„(ЖЁ)= @„„-\(/). Moreover, as by (9), Q, is decreasing on [0,<2[, we get | 
т 20. | 
(Lese Gym) 0+) m [roa | 
nd | 
«| 89) ) e (0) |f (dt | 
| - dt | 
< EEN | 
er = &m+1,0 (f). (1+2) 2 | 
: | 
! 
Hence Q,,,9(Zf) < 28. (7). This ends the proof. | 
The generalized Fourier transform of a suitable function fon R is defined | 
Бу | 
ns. ene А(х)ах, XER. | 
E | 


Remark 2.3. According to (7), (8) and (10), the generalized Fourier кане | 
| 


is well defined on S? (R). | 


M 


Proposition 2.2. For all f € 5° (В), if 
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ANAA- Fad (fo)(), (11) 

where F stands for the Fourier transform related to the differential operator д, 
^ 


defined on S? (R) by 
5,(8)0)- f. h(x)o; (x) A(x)dx, AER. 
Proof. If f є S? (R), identity (11) is obvious. Assume f € S; (IR). By using (3),(4),(5) 


and by integrating by parts we obtain 


^ (0097 ph OAG) 
mdi 
dup 


| Tf (x) (A (х); (x)) dx 


1 


= Л: Tf (x) Ag, (x) A (x)dx 


= (р) [. Tf (x), (х) A (x)dx 


= (ip) F, (Jf) (A), 
which completes the proof. 


Remark 2.4. For A(x) =(sinh|x|)" (cosh x)?" ,02 B» -1/2, the transform F 


3 


coincides with the Jacobi transform of order (0,8) (see [10]). 


Theorem 2.1. The generalized Fourier transform F, is a topological isomorphism 


between S*(R) and S(R). Moreover, 


F (ENEA eNe ЗАО 


ofall g e S(R). 


Proof. By [14] we know that the tranform F, is a topological isomorphism from 


2 
S; (R) onto S, (В). Then the result follows from (1D, Lemma 2.2 and the fact that 


the map f > Af isa topological isomorphism from S. (R) onto S, (R). The identity 
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L) 


io 


m 


at 


above follows easily from (11). 

Note 2.2. We denote by 

e D, (R),a 70, the space of (7° functions on R supported in [-a,a], provided with 
the topology of compact convergence for all derivatives. 


• D(R) U,,, D(R) endowed with the inductive limit topology. 


« H,a>0, the space of entire, rapidly decreasing functions of exponential type a; 


, 


that is, f € H, if and only if, f is entire on C and for all m=0,1, | 


Pm (f) = sup|(1 + 202 f (A)e 9 205 


Н, is equipped with the topology defined by the semi-norms p,,m=0,1,... . 


m, 


• H=U,,,H,, equipped with the inductive limit topology. 


Another standard result for the generalized Fourier transform F, is as 
follows. 


Theorem 2.2. (Paley-Wiener) The generalized Fourier transform Л, is a 
topological isomorphism from D(R) onto H. More precisely, f € D, (Е) if, and 
only if, F, (f)e H,. 

According to [11] the inverse generalized Fourier transform £,' may also 


be expressed as follows. 


Theorem 2.3. For all g e S(R), 


7; (&Y()- | а (0). (ех) ао (2), 


with 
Х-1р ах 
а40)-(559 | TP (12) 
^ 2 2л|е(А) 
where с(8) is а continuous function оп ]0,co[ such that 
ce(s)” 27 2 5 o as s— oe, (13) 
(8) ~ hs, as s 0, 


for some k,,k, ЄС. 
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Remark 2.5. (i) The tempered measure o is called the spectral measure associated | | 
е .5. 

with the differential difference operator A. 
(ii) For A(x) = (sinh IM (cosh xp" ‚®>В > -1/2, we have 2-3 


egre P[(asP+i+is)/2\r[(a—B+1+is)/2] 


5> 0. 


The next statement provides a Plancherel type formula for the generalized 


Fourier trasform F, . 
Theorem 2.4. For all f,g¢S’(R), 
Ї.Г(5)8(-5)А(х)4х-1|, 5,()097,(6)0)40(3) 
Proof. By (11) and (12), | 


dA 


(Ај 


1,5,()095,(6)0)46(3) _ J, Falf.)(4) Fa(e.)(A) 


р ал 
m A+ f АТ fo А. АТ 80 À 2 
METIR ORN 7 


=K +K}. 
By a Plancherel formula for the transform F a (see [14]). 
Ki = |. A (а), (х) А(х)ах 
Moreover, Бу using the identity 
F,(Ah)(A)=-(2* +p") л, (h(a), he S2(R), 
(see [14]), and by integrating by parts we get 


к, = |, A (ANAE (22,0) ^ 


| 10 
= |А) (а) (е) (х) A(x) ax 
= f Ал) ()2(,)()ах 


=—{_fo(x)go(x)A(x)dx. 
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E | — AE 


инээх ээг салж 


Hence 


K, + Ky = Ї [f (x) (x)- fo (x), (x) ]A(z) (z)dxf. f( f (x x)g( -x) A (x)dx . 
This completes the proof. 
Note 2.3. We denote by 


Зи? (В, A(x)dx),1< р <°, the class of measurable functions f on R for which 


Ifl < оо, where 


ІА, = (Ре 
апа ПА... d Ifl... 


1р B 
x)dz) ? if p<, 


І? (R,lol),1< p < e» be the class of measurable functions f on R for which 


Ifl plo) © °°, where 


2211 


and [flag - Il. 
Remark 2.6. By (8) there isa positive constant k>0 such that 


IZ (AON 580400, 


л) 0А)“, | if p<@, 


for all f eL (R, A(x)dx). | 
Lemma 2.3. For all f € L' (R, A(x)dx) and g eD(R). | 

Jf C) A G)ds -[, 5. (00095, (60949 0. | 
Proof. Fix g € D(R). For f € L'(R, A(x)dx) put | 


Р) = |Р) C) A (od 


and 
= AAAA (092502. 


In view of Theorem 2.4; l (f) zi (f) for each f € D(R). Moreover, 
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IG) < 11.181, 
апа 
I C) А, Jal (в): P del) 


by virtue of Remark 2.6. This shows that the linear functionals /, and /; are bounded 


on L (R, A (x)dx). Therefore 1,:41,, and the lemma is proved. 
An immediate consequence of the lemma above is 


Corollary 2.1. The generalized Fourier transform F, is injective on 


L (R, A(x)dx). 
For t>0, the Gaussian kernel E, associated with the differential-difference 
operator is A defined by 


E (x)= f е "5. (-x)do(A), xR (14) 
This kernel enjoys the following properties. 


Proposition 2.3. (i) E, € S^ (R) and 
E -(X p?) 
F,(E,)(A) =e ‚гай Aem. (15) 
(i) Е, is even, positive and Ї E, (x) A(x)dx =1. 


(iii) ^ The function u(x,t)=E,(x) is С" on R x (0,54| 


: and solves the partial 
differential equation 


burst) Lula), 


where A is given by (5). 
(iv) There are two positive constants С,0) and C,(t) such that 
x 2 


C, кте (x)<C, (t) 7 9 (16) 


(v) Let pefo, 


co[ . Then there exists a positive constant M(p,t) such that 
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f 


(E.G) < M(p,t)B,,,(x), an 


Proof. Assertion (i) follows directly from Theorems 2.1 and 2.3. A combination of 
(3), (12) and (14) yields 


E ()- ы =. 
2 j n|c(4.) (18) 
But according to [6], the right hand side of (18) satisfies (ii), (iii) and (iv). According 
to our assumptions on the function A, there is a constant k>0 such that B(x) zh 


for all x e R. The majorization (17) is then an easy consequence of (16). 


3. Hardy and Cowling-Price theorems. The following technical lemmas 
will greatly simplify the proofs of our main theorems. 


Lemma 3.1. [1] Let g be an entire function on C. Suppose that 


le (x) < м(1 + le y е(Вез for all ze C 


and 
lg(x)|S M for all xem, 
for some a, М»0 and m e N. Then g is constant on С. 


Lemma 3.2. [1] Let q € [1,°°| and g be an entire function on C. Suppose that 


f.le ax <= 


and 


|g(x)|< м(1+4))" e^ F9" for all ze C, 


for some a,M>0 and m e N. Then g=0 on C. 


Lemma 3.3. Let q є[1,°°[| and g bean entire function on C. Suppose that 


|81,6 < 
and 
le (e) мн) e for all zec 


for some a, M>0 and m € №: Then g=0on C. 
Proof. By (12) 
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20 id 
lel. s fale 990) 

E: a A -ip ах Ў 

ЫР aei 

a dÀ fe 

е 0) ale (nh le 


?**! for all AER 


-2 
According to (13), there is a constant k>0 such that lc (А) >k [A 
with n >1. Therefore 


g(a) Ap" ах» 23 


Ч | 
m he 22 


lel? 5= ANT 


which shows that I|. < со. The result is now a direct consequence of Lemma 8.2. 


Lemma 3.4. Let a,b>0, d 2 1,y € R and 


= -0(Х- Я 22 
g(y)= Ї ee) (14x) ё"фх, y>0. 
Then there is a positive constant C such that 


g(y) s С(1+ yy e» forall y20. | 


Proof. By the convexity of х we have 


g(y)=e™ Ї ge (1+2+Ьу) dz 
« 2) ои (1 +|2|+ Ьу) @5 
< em |7 ena Hil (1+|2|+by)" at 


= er ЇГ ew tie (1+2+ Ьу)“ аг 


< conste [" “хи (1 +2%+ (57) m | 


! 
= toy | [ауе 
const.e Ї e e La a )dz +(Ь eg 
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вэ 


< const. (1+ у“ jes 1 


< const. (14 y)" e” 
which achieves the proof. 


Lemma 3.5. Let 1 < q <œ and а>0. Then there is a,positive constant C such that 


for all А= £ ine -1Ё: 
Ф [Ee CHA); 


«сүү e ms 


ш (Е.Ф | 


Proof. Ав the function 1/./В(х) is bounded, it follows from (8) and (16) that 


| 
|E, м. (5) Ф. (x)| < const. (L p) (1-4 ёс “лэ 


"2 -а||х|-— 
= const. (1 + I + la) esse 


which proves (i). For q < »» we have 


2 1/q 


заса 
[УКОЛ = const.(1+|A|)e" ^" IE of 2) (1+х x)"e eU P dx 


2 (2-a)elni a)elni | 
< const. (1 + |2) (1 + шүр ча? 204 1 


1 
2 (2-8) | 
< const. (1- үе“ 8 


by virtue of (10) and Lemma 8.4. 


Lemma 3.6. Let 1<р,р S” such that 1/р+1/р'=1. Let f be a measurable 


<œ for some a7 0. Then the generalized Fourier 


function on R such that |Ёи./1,, 


transform of fis well defined and entire on С. Moreover, there is a positive constant 


C such that for all A-&*ime R &iR: 
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ЭР 
22 | 
0 1550) о Ma if p=]; | f 
P Qro | 
Gi) JF (NAS C( (1+|A)) e 22 2c» 7021 | 


Proof. The follows easily by using Lemma 3.5, Hólder's inequality and the | 


derivation theorem under the integral sign. 
We can now state our main results. 


Theorem 3.1. Let 1€ p,q<~. Let f bea measurable function on R such that 
Ej, f € L (R, А(х)48) (19) 
and 
e" 7, (f)e L (Е), (20) | 


for some positive constants a and b. Then | 
If ab 71/4, we have f=0 almost everyhwhere. 


If ab» 1/4, for all t € |b,1/(4a)| , E, satisfies (19)-(20). 


Proof. We divide the proof in two steps. 
Step 1. ab» 1/4 


Let t€ ] 1/(4a),d[ and g(A)=e” F (f(A), AEC. 

By Lemma 3.6, g is entire on C , and there is C>0 such that 
мО)ес(1-р| ee 

for all Ae C. Furthermore, 


а 04 


< оо 


aoj : GE (f) ajloj 


(0) — If q«», it follows from Lemma 3.3 that gA) 
5,(7)0)-0 for all AER, Therefore, f= 
(1)  Ifq2e, 


=0 for all лес. That is, 


0 ae. on R, by virtue of Corollary 2.1. 
then by Lemma 3.1 there is a constant K eC such that g(A)=K 
for all A e C. That is, 39 (f)(A)= 


But due to assumption (19), thi 
Step 2. ab<1/4 15 is impossible unless K=0. Thus f=0 a.e. on В. 


Let te Jb (4aj[. 
(a,£) such that 


Ке for all Aem. Hence, f = Ke' E, a.e. on R- 


[By (1), there are two positive constants C,(a,t) and Cz, 


Е E (x)E, (х) « C, (ege ee | | 
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15, 


23 


for all x e R. This shows that EE, e I? (ЕЁ, A (x)dx). Moreover, 

a э 
-0-5)2 ae 
aloj 


b 2 2752 
le 7, (E, 1, loj =e" lle 
» 


pune of (15) and the fact that о is tempered. This ends the proof of Theorem 
8:1: 


Theorem 3.2. Let 1< р<2 and 1<q<o. Let f be a measurable function on р 


satisfying (19) and (20) for some positive constants a and b. If ab=1/4 then f=0 
almost everywhere. 
Proof. Let 


g()- e" F (FA) VEC. 


Let р” be the conjugate exponent of p. As by hypothesis р'> 2, we deduce from 
Lemma 3.6 that g is entire on С, and there is C70 such that 


|g (A)/sC(1+ А) een 


for all 1 eC. The rest of the proof is now analogous to Step 1 in the proof of 
Theorem 3.1 
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ABSTRACT 

In this paper, an S-I-S epidemic model is investigated for carrier dependent 
infectious diseases such as measles, cholera etc. by considering explicitly the 
cumulative density of ecological factors, which are conducive to the growth of 
carrier population. It is assumed that the density of carrier population is governed 
by a generalized logistic model, the growth rate of which increases as the 
cumulative density of ecological factors increases. This density assumed to be 
human population density dependent. The model is analyzed by using qualitative 
theory of differential equations. It is shown that as the parameters governing the 
cumulative density of ecological factors increase, the number of infectives in the 
population increases. The model is discussed through the basic reproduction 
number and through numerical simulations also. 

2010 Mathematics Subject Classification : 37C75, 92805 
Keywords and Phrases : Immigration; Logistic Model; Stability Analysis; 
Simulation; Carrier Dependent Infectious Diseases, Ecological Factors. 

1. Introduction. There are many carrier dependent infectious diseases 
which afflict human populations mainly in the developing world. Various kinds of 
carriers such as flies, cockroaches, ticks, mites, etc., grow and survive in ecological 
habitats consisting of plants, bushes, parks in residential areas, water storage 
tanks, ponds etc. It may be pointed out that their population grow in the mounds 
formed by the roots of the bushes and plants as well as on the branches and the 
leaves and their populations increase further as the bushes and plants grow and 
become more dense. These carriers transport bacteria of cholera, dysentery, measles 
and gastroenteritis, etc. discharged by infectives of these diseases in the 
environment to susceptibles and thus spread these diseases in human population 
indirectly. 

Ecological interactions 
together or in sequence [1],[8], 
growth of suburban communities an 


have contributed to the rise in reporte 
spread of Lyme disease in United States [3],[6],[14]. It 


can be complex with several factors often working 
[11], [18]-[20],[22], [23], 1241,(261, 1281, [35] e.g., the 
d the reforestation of the eastern United State 
d tick bites and infections and thus to the 
is further noted that densities 
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of these ecological factors may increase due to various human population density 
related factors such as housing, agriculture and industrialization. 

The modelling and analysis of various infectious diseases have been 
conducted by many researchers in the past as Bailey [2], Brauer [4], Brauer and 
Vanden Driessche [5], Chavez and Yakubu [7], Diekmann and Heesterbeek [9], Dietz 
and Eichner [10], Gao and Hethcote [12], Ghosh et.al. [13], Hethcote ([15]-[17]), 
Kalajdzievska and Li [21], May and Anderson [25], Murray [27], бай et. al. [29], 
Singh et.al. (8301-1320), Zhang et.al. ((33-34]), However, it is noted that very little 
attention has been paid to the study of such diseases by considering effects of 
bacteria transporting carrier population, the density of which increases due to 
natural or human population density related ecological factors. Therefore, in this 
paper, an SIS model with constant immigration for a carrier dependent infectious 
diseases is proposed and analyzed by considering explicitly the effects of ecological 
factors which depends upon the density of human population. To be specific in the 
modelling process, it is assumed that the density of carrier population follows 
logistic model, the intrinsic growth rate and carrying capacity of which increase 
with the cumulative density of ecological factors. The cumulative density of 
ecological factors (such as biomass of leaves in bushes and plants etc.) is governed 
by a generalized logistic type equation, which depends upon human population 
density related factors. 

This paper supplements the research work done by Singh et.al. [32] 
таша the spread of carrier dependent infectious diseases with environmental 
effects. 


2. An SIS Model with Ecological Effects. Let the total human population 
density N (t) be divided into two sub classes: the susceptible density X (t) and the 


infective density Y(t), thus N=X+Y. Let B(t) be the cumulative density of 


ecological factors (such as biomass of leaves in bushes, plants, etc.) favourable to 
the growth of carrier population and is governed by a generalized logistic equation. 
Let C(t) be the carrier population density which is governed by a logistic model 


whose growth rate coefficient s(B) and modified carrying capacity L(B) increase 


as the cumulative density of ecological factors B(t) increases. It is implicitly 


assumed here that the number of agents of the disease from the environment 
transported by carriers is proportional to their density. 

Keeping in view the above co | 
action interaction, an SIS model is р 


ах 
(2.1) Gp A BAY -AXC+0Y - ax , 


nsiderations and by assuming simple mass 
roposed as follows: 

DY 

ap PAY +AXC-(u+atd)y, 


DN 
—, = A-dN-oY, 
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-s(B 
mu) jay o 
ОН 
qu^ Pol K)B' -7,B +r,BN, 


X +Y =N,X(0)>0,¥ (0) 2 0,N (0) » 0,C (0) > 0, B(0) > 0. 


For feasibility of the model (2.1), we must have s(B)- s, » 0 for B» 0 and л» 5. 


In model (2.1), A is constant immigration rate of human population, d is natural 
death rate constant, D and X are the transmission coefficients due to infective and 
carrier populations respectively, œ is the disease related deth rate constant, v is 
the recovery rate constant, s, is the deth rate coefficient of carriers due to natural 
factors as well as by control measures. Further, ry is the intrinsic growth rate 
coefficient of B(/), г; is the natural depletion control rate of В(1), r is the growth 
rate coefficent of B(t) due to human population densty related factors and K is the 
carrying capacity of B(/) which is assumed to be constant. 

In (2.1), s(B) is the growth rate per capita of the carrier population density 
such that s(B)-s, is its intrinsic growth rate as compared to the usual logistic 
model. In view of the assumption that the growth rate capita s(B) increases as the 
cumulative density of ecological factors B(£) increases, we have 


ds(B) 
(22) 8(0)-в,»0 and "ag ОШ 


where sọ is the value of s(B) when it is independent of B and s} >s,. 


Similarly, the fourth equation of model (2.1) governs the density of carrier 
population, when it is compared with usui logistic model, gives the carrying 


capacity as (1(B)/s,)(s(B)-5] where L(B) is the usual carrying capacity of 
carrier population density. It is assumed that this modified carrying capacity, L(B), 
increases as B(t) increases, so we have, 


dL (B) 
= 20 
(2.3) L(0)=L,>0 and 755 , 


where L, is the value of L(B) when is independent of B. ` E 
It is leas from the fourth equation of model (2.1), (2.2) and (2.3) that when B=0, 


the carrier population density increases.in its natural environment and it tends 
to L,(1- (s, /s, )) which may become zero as 5, > sa: : 

3. Equilibrium Analysis. For analyzing the model (2.1), we consider the 
following reduced system (using X+Y=N ,(ү«М,Х»0): 


вл) DY -p(N-Y)Y+A(N-Y)C-(v+a+d)Y, 
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nB” 


dB = аагар 


тор 
where n -n 570. 

To analyze the above model, we first need the following lemma. The lemma 
establishes a region of attraction for the system (3.1). 

Lemma 3.1: The set 


Q={(Y,N,C,B):0<Y<N<A/d,0<C<C,,0<B<B,} 


m? 


attracts all solutions initiating in the positive orthant, 
where B, =(K/7){r+(7,A/d)} and C, ={L(B,,)/s)}{s(B,,)-s,}. 
Proof. From the second equation of (3.1),we have 


ам 
п ^^ dN-aY sA-dN which implies 0 x Y<N<A/d. 


Further, from the last equation of model (3.1) we get 
2 
EE d 
which is a logistic type equation. Hence we have, 
0<B<B, where B, = Inr (4)/d) " 
Similarly, from the third equation of model (3.1), we get 
dC e 


сае e rer fe] 


which is the logistic type equation and gives: 


0<C<C, where C, = {L(B,,)/s)}{s(B,,)- sı}; provided e(05 ecw 


The region of. attraction, therefore 18, as given by Q defined avove 


It may be noted here that for the feasibility of the model s(B,,)>s 
Ы т i 9 


ng theorem: 
a: 


The results of equilibrium analysi i i 
i YSIS are given in the followi 
| Theorem 3.2: The system(3.1) has following six ius. 


| 0) E, =(0,A/d,0,0) 


(0 — E, =(0,A/d,0,8,,), where В, - (Inr (пАу4) 
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CO (D be 


| А 7777 ВА 
1 »=(Y,N,0 ase ee Л 
| (ii) 2, ( ,N, ,0) where exists if d(Us oma) 51 where KANA basis 


| reproductive number, Y ={BA-d(v+a+d)}/P(a+d) and 
| Ñ=[PA+a(v+a+d)}/B(a+d); 

(iv) — E, - (Y, N,0, B) which exists if R>1, 

where B=(K/n)(r+r,N),Y,N and К are defined above; 


ma (v) 19) = (Y,N,C,0) where 


Y =| B+ [B хад(0/4)00/4) | (6(o.2)1], 
B {(BA/d)-(v+a+d)-AC(1+(a/d))}, 


| C - (L,/s,)(s,-,) and Ñ 2 (A-oY J/d; 
(i) Е, -(Y*,N*,C*, B") the endemic equilibrium, which exists in the subregion 


О, ={(Y,N,C,B):0<Y < A/d,A/(at+d)<N<A/d,0<C<C,,,0< B<B,} of o, 


e 
" d(v+a+d) 
if the equilibrium point Е, exists then E; and Ёз will not exist. 

Proof. The proof of existence of Ey,E,,E2,E3 or E, is trivial. The endemic 
equilibrium point Ej is given by the solution of the following set of equations 
obtained from (3.1) by putting left hand sides to zero: 


| (32) QY*«Y (vea d)- BN € AC] - ANC - 0, 
(88) Y=(A-dN)/a, 
(34) C={s(B)-s,}{L(B)/s}, 


(3.5) B=(K/m)(r+7N). p 
Eliminating Y from (3.2) and (3.3), we get the following relalion: | 
(86) F(N)=(B/0)(A-dN) +{(A —dN)/a}{(v+a+d)-BN +AC}-ANC =0, 
where C can be obtained in terms of N by equations (3.4) and (3.5). 


provided the basic reproductive number, Ё <1. It is noted here that 


BA 


It is noted from (3.6) that F(A/(a+4)) is positive provided R= d(v+a+d) = 


and F(A/d) is negative, bacause 5(В„)—5,> 0. 


| Thus, it is clear that there exists a root N* of F(N)=0 in A/(a+d)<N<A/d 
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provided К «1. Further, this root N* will be unique if 


dF(N) o for A/(o-*d)« N «A/d. 
dN 


(3.7) 
To show this, we differentiate (3.6) to get 


(3.8) E = «(Ba - av) (£t oca) - t + AC} (2 -dN) 


АКР, s з N(a+d)-A 
-AC 25: s,}L'(B)+s ey" 
0'0 
dF(N .. RA 
It can be easily seen from (3.8) that d ) <0 provided R= d(v+a+d) <1 in the 


region Al(a+d)< N<A/d. 
| 

Now, knowing the value of N*, the values of Y*,C* and B* can be uniquely | 
determined from (3.3), (3.4) and (3.5). Hence, E; exists in Q, provided R<1. 
Remarks. 
(i) Ме note that 

dF(N) 

dN 


= —(d/a){BY *+(AN*C*/Y *)}-BY *-AC* 


E5 


-(AKr, Iss (B*)L(B*)+{s(B*)-s,}L(B *)|(N*-Y*)<0. 
Thus, it is clear that the above condition (3.7) is automatically satisfied at 


the equilibrium point E,. | 


»0 


E; 


= dY 
(ii) We note that mE 


e 20 
and dr " o 


On substituting the value of N from (3.3) into (3.2) and (3.5), we get 
py? +Y[(0+0+d)-B{(A-a¥)/d}+AC]-2{(A-a¥)/d}o=0 
and B=(Krin)+(r,KIm){(A-a¥)/d}. 


On differentiating the above equations wrt. 


25) Q а 1 
3 r' and on simplification, we 


dY 
(89) (ВУ «(o8/d)Y +(ANC/Y) + (20/4)|-3(М-ү)96 and | 
Р 


аВ | 
(3.10) pu In) - (or, K да). | 
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| Now from (3.4), we have | 


ас 
(3.11) dB a I(B)/s, , 


аге I(B)={s(B)-s,}L'(B)+L(B)s‘(B). 


Further, we note that 


dC dC dB 
(3.12) Л» dad- 
Using (8.9)-(3.12), we get 
he | 813) 27 - HAAK a мн 30 
BY *4(AN*C#/Y*)+(oB/d)¥" +(AC*a/d)+ ^ A I(B*) 


| Thus it is seen here that as the growth rate coefficient of cumulative density of 
ly | ecological factors increases, the infective human population density increases at 
the equilibrium £E;. 


from (3.13) into (3.9), we get 


dY 
Further ,on substituting the value of dr |; 


E Ed K[BY *«(AN *C*/Y *)+(aB/d)¥ *+(AC*a/d) |1(B*) »0 
A 807, СЕ СА нову «etapa Ч ко 


smd 


which shows that as the growth rate coefficient of cumulative density of ecological 
| factors increases ,the carrier population density increases at the equilibrium Ez. 


In the similar manner, we can show that _ 

A(N*-Y *) (KN */syr, (82) 2 
s (Му *)aKr (B*) ^ 

Ma BY *4(AN*C*/Y *)+ (aB/d)Y *+(cAC*/d)+ - eae 


Syd 
Thus, it is clear that as the growth rate coefficient of cumulative density of 1 
ecological factors due to human population density related factors increases, the 
infective human population density increases at the equilibrium #5. 
(iv) Also, we have 

(KN*/syn)[BY *+(АМ *С* гҮ *)+(оВ/ ду * 0*0) |B") P 


Ne 


dC 
| S * A(N *-Y *) Кол, E 
| dn Es СЕ =o JY *)+(aB/d)Y *+(AC ша) 4 snd he ) 


а eUam TENERE TNI meat menn gne 
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which shows that as the growth rate cofficient of cumulative density of ecologica] 
factors due to human population density related factors increases, the carrier 
population density increses at the equilibrium Е;. F 
4. Stability Analysis. In the following, the local stability results of 
equilibria Ep E ЕЕ з and E; are stated,the details of which are given in the 
Appendtx A. 
Theorem 4.1: the equilibira Ey, E} Ёс, Ез and E, are locally unstable and the 
if exists, is locally asymptutically stable provided the following 


eqilibrium E;, 
conditions are satisfied, 
(41) oC оар * /3, 


2 


* Í 
(4.2) 30" (N*-Y *) D (B*) (B+ s (В) К° < аВ?в2У sg | 
Remrak. It is pointed out here that the inequality (4.2) is automatically when 
s(B) and L(B) are independent of B i.e. when the density of carrier population is 
not affected by the cumulative density of ecological factors. This shows that | 
ecological factors have destabilizing effects the system. 
The global stability result of E; is given the following theorem: 
Theorem 4.3. In addition to assumptions (2.2) and (2.3), let S(B) and L(B) satisfy 


0<5'(В)< p and 0<L’(B)<q for some positive constants p and q in О,, then £s, 


if exists, is globally asymptotically stable in О, provided the following inequalities 
are staisfied: 


(43) OC? «2gp!y + 3, 


2 
ao RANTS Bo pes B къз capiens 
0 


Remark. It is noted here that if p=q=0, the inequality (4.4) is automatically 


Би which shows that the cumulative density of ecological factors have 
destabilizing effects on the system. 


Proof. See Appenix В. 
The above theorem impl 


: y that under certain conditi i itv of 
carrier population increases, th itions if the density 0 


en the number of infectives in human population 


| 
| | . 5. Simulation. In this section 


we anal e | 
vaules of parameters to show the glob ee QD for appropria ш 


al stability behavior of E.. 
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For numerical simulation we choose s(B) and L(B) as follows: 


(5.1) s(B)=s,+s,,B and L(B)=L,+L,,B. 3 


In the following ,we shall first show the existence equilibrium point and then its 


global asymptotical stability behavior by taking the following values of the vari 
ous parameters: 


В=5.1х10"',0=0.012,0= 0.0005, d = 0.0004,А.=2.1х10*, A = 10,5, = 0.9,s, = 0.6, 


su 22.0x10 ^, L, = 100000, L, = 4x 10% 7, =0.9,r = 0.7,r, = 0.225 x 10, К = 26000. 


The values of N* and Y* are determined as the point of intersection of the 
following equations (which are obtained by eliminating C and B from (3.2) to (3.5): 


(52)pY? +¥[(v+atd)-BN +A{(s, Б +8,,(K/7)(r+7,N))(L, +L, (K/r)(r+7N))/sy} | 


-АМ f(s — 5 + 5|| (К /т,)(т + r,N))(L, tL, (K/m)(r-+r,N))/s,} -0 
and 
53) Y -(A-dN)/a. 


Using the above numerical values of various parameters, these equations are plotted 
in Y-N plane by using MATLAB and are shown in following Fig.1: 


6010 


4.5 


ш 
an ~ 


[^] 


N 


Total Human Population Density (N) 
О 


- 


0 05 1 15 : 2 25 
Infective Human Population Densily (0) -. x 10 


Fig. 1 : Existence of (Y*,N*) 


| Here, the point of intersection in the first quadrant gives Y*=3411.91, 
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34 


N* =20735.10 and the values of C* and B* can be obtained from the equations 
= 5. 


anm 
| F 
үл 

(3.4) and (3.5) respectively. . . 

Further, the endemic equilibrium point E, can be obtained directly by | 

numerically solving the equations (3.2) to (3.5). Thus, solving these equations 

along with (5.1) (with MAPLE) and using the numerical values of various 

parameters given above, we get 

Y* = 3411.912520, N * = 20735.10935,C* = 38126.70052, B* = 21570.00433 . 

Hence the existence of the endemic equilibrium point is shown for the given set of 

parameters. The variational matrix corresponding to this equilibrium Æ., is 

obtained as, 


-0.006602 0.002532 0.000363 0 

-0.0005 -0.0004 0 0 

о Ят O 0.34200 0.07600 
0 0.048370 0 -0.74400 


The eigenvalues of this matrix аге:-0.0064, -0.0006, -0.7440, -0.8420 which are 
all negative. Hence [7, is locally stable. 

It is also pointed out here that for the above set of parameters, the conditions 
for local stability (4.1,4.2) and global stability (4.3,4.4) have been checked and they 
are satisfied. Also for the above set of parameters, a computer generated graph of 
Y vs N for the model (3.1) is shown below which indicates the global stability of 
the point (Y*,N *) in the Y-N plane : 


2000 4000 6000 800 10005 13000 14008 
b . Fege 2; GlobakStabitiegion Glos eGa 
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Further, variation of infective human population (Y) with time (0) has been sion 
for different values of parameters r and го in Figs. 3 and 4. 


c 

E Өс 27 
E r-2. 

a 

©. 

© 300 г-1.7 
5 г=0.7 
ES 

т 2: 200 

ge 

бүс 

= 100 


0:775 4 6 8 10 12 14 16 18 2 
Time (0 


Fig. 3. Variation of Y with t for different values of parameter 
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8 r, = 0.00008 
Ss = r, -0:00000225 
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I $5 200 
2$ 
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0 777777 в 8 10 12 14 16 18 20 
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Fig.4. Variation of Y with t for different values of parameter г» 


From above figures, it can be seen that Y increases as г ог Го increases. Thus, it is 
clear that the infective human population density increases as the growth rate 
coefficient of cumulative density of ecological factors and growth rate coefficient 
of the same due to human population density related factors, increase. 

6. Conclusions. In this paper, we have proposed ара analyzed as SIS model 
for a carrier dependent infectious disease by considering «оошо effects. We 
have assumed only those ecological factors which are conducive to the growth of 


density of the carrier population. It is assumed that the cumulative density of 


such ecological factors is governed by a generalized logistic вазада which is 
population density dependent. The equation governing the carrier population has 


been assumed to be a genralized logistic model, the intrinsic growth rate and 
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carrying capacity of which increase as the cumulative density of ecological factors 
increases. The model has been analyzed analytically as well as by computer 
simulation. The effects of parameters governing the ecological factors conductive 
to the growth of carrier population have been found to increase the density of 
carrier population, leading to fast spread of infectious diseases. It is noted that | o 
the ecological factors have destabilizing effects on the system. Further it has been | 
found that an infectious disease becomes more endemic due to immigration. | 
Appendix A. | 
Proof of the Theorem 4.1. In the following, the local stability behaviour of each | 
of the five equilibria Ej to E, is studied by computing variational matrix and the | 
endemic equilibrium point E; is studied by using Lyapunov's theory. 

The variational matrix M; corresponding to the equilibrium points Е, | 
1=1,2,3,4,5 is given by 


220228088 
~ 
= 


BN-2BY-AC-(v+a+d) BY«AC — A(N-Y) 0 | 
р», -d 0 0 | 
М,- 0 0 8(В)-8 228 s'(B)C о 
( ) leery ( B) ( ) 12 ( В) ( 
0 nB 0 Е 218 trN j 
K ( 
For equilibrium point Ер, we note that one of the eigenvalues of My is r+(r,A/d), 
which is positive. Thus Ey is unstable. d 


We note that one of the eigenvalues of M ; is 5(8,) 


р -5,, which is positive. Thus £} ) 
. is unstable. й ; | 


We find that one of the eigenvalues of M, is гэл 
2 


N , which is positive. Thus E, is б 
unstable. 


We find that one of the eigenvalues of E, is s(B) 


—$,, which is positive. Thus Е, is 
unstable. Е ^ 5 


Е СЕ 3 а | 
or equilibrium point E,, we find that one of the eigenvalues of M, is r+r,N; 3 
4 2 


which is positive. Thus, E, is unstable. 


Stability Behaviour of E, - (Y*, № сж B*) ; 1 
5 | st 
ince the nature of Е, cannot be seen easily from the variational matrix, | А 


we st it’ ili i 
study it’s local stability by using Lyapunov’s theory. For this, we use following Р 
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positive definite function 


| 

| 
ter | x2 i 
el AD У (№ /2)(Y -Y * +(k,/2)(N-N *)’ + (k,/2)(C-C * «(&,/2)(B- B*y. 
of | Thus, dV /dt along the lineraized system (3.1) can be written after rearrangement 


iat | oftermsas: 


1 аг" (BY *-ha)n- (RC NY) - (spy */3) 
ich -[(8Y */3) у? - (AC *) уп +(k,d/2)n? | 

he : 

E [or "Во ee Ne Es J d 


M е} ые 


-[(%В */2К)ь° - (kyr,B*)nb + (&d/2)n |. 


C h Ар=1апа k =(a/BY *), we get following inequalities for ү to be negative 
definite, 
(А.2) oA?C* < 2dp*y * /3, 


1 ) 2 БЭ жез х 2 B* K 
| (А.З) za c Ж 5 < k, SS 2 % d L 2 ky, 
50 2 2 ; So 
L(B*)C*4s +5 L'(B* Ü 
zl (*6* 107): 6 
в, (А4) ^ SESS E йг x. 
Now, by choosing k, -r,d/2Kr; B*, the inequality (A.4) stisfied automatically. 
is Further, we can choose hy satisfying (A.3) provided that | 
2 | 
7,  80d4(IV*-Y*)*17(B*)\s'(B*) +L (89) К 227 $B'Y *' rj is satisfied. | 
V, (N*-Y *) P(B*Js'( I(B*) 
_ It can be seen that ү is negative definite under the conditions (4.1) and (4.2) as 1 
| Ё 
| Stated in Theorem 4.1. Hence the result. | 


ix, | Appendix B | 
ag Proof of the Theorem 4.3. We prove the theorem by using the following positive 1 
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definite function 


(B) v-A(Y-Y*-Y* inY/Y*)« Aw- N*) +k (C-C*-C* inC/C*) 


+k, (B-B*-B* in B/B*). 
Differentiating (B.1) and using (3.1), we get 


(B.2) < = (&B-&a)(Y -Y *)(N -N *) -(kACN /YY *)(Y -Y *)’ - (kB/3)(Y -Y Y. | 


(8/8)(¥ -Y È -(kAC/Y *)(Y -Y "(N - N *)+ (d/2)(N - N *'] 


2L(B") 


al 
I У-Ү) (Ау *)(N*-Y*)(Y -Y )(C- C) 59 _ c-o] 


| 1,5, 


-[06/2) (мм) -kr (N - N *)(8- B9) (yn, /2K)(B—B 1 
where f(B) and g(B) аге defined as 


s(B)-s(B*) pap 
(B3) /(в)-| B-B* ' 

ds 

ав” DID 

. 4 

and 

SURGE) ae 

B-B* твір) 2*B* 
(ва) &(B) ES L(B)L(B*) 

L’(B) ав? B=B 


(B.5) 0< s'(B)< p and 0< L'(B)« 
Hence, from (B. 3), (B.4) (B.5) and us 


(86) |f(B)<p ana Ї (8) < ал12. 


q for some positive constants p and g. 


ing the mean value theorem, we have 
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2L (B9) € (В) +5,Cg(B)}(C-C*)(B- Be)+( 22 Jm 3 sy 


I 
| 


[ 


— 


Cl 


(В 


нэ... оо... 


Choosing №, =1 and k, =«/B, we have following i x 
8 inequalities for ү t 7 

| definite V to be negative 

(B.7) oA Ch «2dp^Y *° /8, 


| 3 od L(B*)(N *-Y *) 
25. 5-2 hb, < TS 
(888) 9 зүс <њ <= 
| Sof KL(B*){p+s,C,,q/ BY 


| B9) h«nd/(Kr;). 


| Ву choosing k, =r,d/2Kr, , the inequality (B.9) is automatically satisfied. So, we 
| сап choose £, satisfying (B.8) provided that 


2 
| ЗОЛ (B*)(N*-Y Урьд! K^rj «rj dB's;Y * is satisfied. 
| | Hence, E; is globally asymptotically stable if (4.3) and (4.4) are satisfied as stated 
| in the Theorem 4.3. Hence the result. 
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) ABSTRACT 

| In this paper we have obtained solution of certain new а integral equation 
in which Saxena's I-function occurs as a part of the kernel. Erdélyi-Kober fractional 
integral operators have been applied to transform the kernel in desired form, so 

| that the conditions of unsymmetric fourier kernels are satisfied. These results 
may be useful in many applications of integral equations. The special cases are 
obtained from the known results. First we shall transform the given equation in a 
desired form, by making use of fractional integral operators. Then the integral 
equation thus obtained, will be split into a pair of equations, with the help of 
Laplace transforms. The resulting integral equations will, then be solved by using 

| the formulas, already derived. 

| 2010 Mathematics Subject Classification : Primary 33C70;Secondary 45E05 
Key Words and Phrases : Saxena’s I-function, Erdélyi-Kober fractional integral 
operators, Integral equation, Laplace transforms. 

1. Introduction. Solution of various generalized sets of dual and triple | 
| integral equations have been given by Charles Fox ([1],[2]), В.К. Saxena [5] and | 
| VP Saxena [6]. In each case the responsibility of final solution goes to integral " 
| equation, having H-function as the kernel of the type corresponding to the set. | 

But due to non-availability of much inversion formula the, kernels аге more ог 


| less restricted to some standard forms. ; | , 

In this paper we shall give solution to the integral equation having 

considerably more generalized kernel. We shall use Laplace and Ее 

| Laplace transforms of Meijer [4] and operators of fractional integration. The 
! 


| generalization of Fox’s H-functions, called, I-function is given by Saxena [7], its 


| definition is given as 
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(a, aur 1 3, n;(a; Sj ji +l,pi 


(858,) (aBa) | al p 


(11) I(Z)eIm. | 


e -A-o +a,,6) 
УП 8-5, 58:01 (85-82) 


with the conditions of existence given in (2.17). 
(i) If we put r=1 in (1.1), it reduces to Fox's H-function 


122122 201000 
J1 Ју? Л? Ji [nay pi Оа п] 
I" n 2 Dp “НТ: л Ш п+1,р 
р.а 

? 


o(6)= 


"o, Bj), ini), "(6 B; hc (b. B;) 


where p, = p,q, = 9,9; =4;,0;, =b; =0.,, and 8,-08,. 
(ii) Ifwe put a=b =a;= 


m*l,q 


= п) а, 


2ni 
where 


үхэн, 


S res ree) 


t=1 | ј=т+1 nd 


(iii) If we put r=1a;=B;=a,=8 


function, i.e. 
Inn f (в, d), (аа), s d es Ы (а), (8), 
SOR 6 nc), цоо (бл), (5,1)... : 


2. Problem and Solution. We consider the integra] equation 
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В,-1ш (1.1), it reduces to a function J defined as | 


ji Pj =O in (2.1), it reduces to a Meijer's G- | 


| 


| 

| 

| (2.1) g(x)- f 1 ux 
| 


D rae 
| 45 3 


(1 = ае) ses), a] 
(or), fa e Ба) „| 


Here g(x) is known and f(x) is to be determined. The notation of the [- function, 
used here in different from (1.1). 


(u)du 


Mellin transform of J-function is given by 


| (2.2) МІ ке р СО) 2 Пг, -Aer t-a +8) 
| ? (851), hla- bif); 123 SITI re-s- ГЭ г(а,- : e| 


jzq*1 


The Mellin transform of a function f(x) is defined as 


Q3 F(5)-M[f(u)]- f w"f(u)du 


! The Parseval's theorem of the transform is 


(2.4) f° &(ux)f (и) аи - lim gf (I-E) 


Тә=» От” ov-iT 


| From the equation (2.4) and (2.2), we obtain 


 їгф,-8 0-2 


| @ 5) =“ х*Е(1-&)4&= g(x), 
lu 27145 2, | s 
| sire EU 


| Where F(5) = M[f(u)]. 
Now we shall use following definition of Beta function 

| 5 -p- ..Г(В+сѕ) 

| (9.6) [| (v^ -29) р 1 jile-ble-s-ldy 2 cl (b-B) x” ECT 

(v^ ERU үес, USES 


Replace x by v in (2.6), multiply by 
[ x to co where x »1 to 


| 8ndintegrate under the integral sign with respect to v form 
| obtain A 
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S j hcm Oh = 
Foie март унаанд, EG 

Accordingly from eq. (2.5), we arrive at | 
| 2 S[re-e ед) | 
| ШОШ лд ГЭЛЭГ +) L xt (L-E) 
| j=2 25 
| enai : р l'(b + AE) | 
E oe 
| i=l (/=р+1 j=q+1 
| 1 vh (ий ХУЛ n p! ich (yd, 
3 = fT (b - bh) > 


Let us denote right hand side of above equations Бу g,(x). Then it is evident that | 


(2.8) а, (x) = 9[5, =H /fi :1/f,]g(x), 
where 9; is second operator of fractional integration defined in (2.6). 
We apply the operator again, so that we get 


1 [г г(ь, +Г®г( H * fiE)T (m +AT- a-e) , | 
з 7 ян =F (1-&)dé = g, (x) 
УП г(1-Ь i Fé) TI Plantet) Pler te) 


=1 | j=p+l j=q+l 
where | 


(2.10) в,(х)=%[%—-н„ц„//,:1/}Ь]а\ (а). 


In a similar manner, we have 


ainw &с оо га. 


Again making use of Perseval’s th om 
8 D 
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{ and so on. | 
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| respectively | 
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| (2.18) 8:(3)-1, I| ux ddr 
| | (AeA HOA) l-ba) i: 


| 1-aj,e;) joe), ‚(ооё 
hat | (2.14) &(x)» f, I| ux ( hal ( 1 Je 2 f (u)du 


(eA ea HOA Etpa — ' 


where the functions g,(x),g,(x) and g,(x) are given in (2.12), (2.13) and (2.14). 


| Now, we are in a position to split either of the form (2.13), (2.14), (2.15), 
whichever is suitable, into the pair of known integral equations. We have already 
indicated that the choice of the form depends on the type and number of param- 
(x) | eters in (2.1). The criteria for the selection will be clear when we arrive at the final 
solution. 

First we assure that any one of /,(7=1,2,... р) is equal to unity or reducible 


to unity say f, 21, where p, =0. 
In this case we consider the integral equation (2.11) which may be written in the 


form 


| c (a- 97 е), T : 02 5s f dx i 
их (u)du = g, (x). | 
(x), (2.15) |, 1 T г), а-в), 2 


Now we use the following known result 


gout 
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The I-function converges absolutely in the Ё -plane if 
(2.17) G) X» 0,агө | “лА/2, 


(ii) A > 0,argz| > 14./2, Re(u + 1)«0, 


18197 | . 


Q Р Si Ti 
LIES ЭЕ 8 | 
ј=1 


3 (4,-26)-ms| 2 (А,-ЕД)- 2 (8-8 )8.2| 


74041 Ј=Р+1 


If €=0 then the set of above expressions (2.16) and (2.17) reduce to 


Р Q S, Т, 
БЭ | УЕ, + У 5, 
j=l 741 Шин 17:22! Ј=9+1 


DOPOR Y. (4)- X (@)-8 2.2012] 


J=Q+1 Ј=Р+1 


But by using the formulas 


(0-4), 4683... 


(8,5), E ВР), .) 


(8.5), a ВЕ) 
(1-4 Е) (4585), 


ae 
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(2.19) x?I|x 
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22227 
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Applying equation (2.20) in (2.15), we get 


(0-4,8) 6 53,.л| 


dyd 
(вв), „}, {a- ВБ, эч 


(2.21) (ШЕТ ху 
ae 


1 
and С, (х) = P (1/x). 


Reversing the order of integration in (2.21) we arrive at the following set of inte- 


gral equations. 


(2.22) [rese (u)du = h, (x), 


| 122277:2-2) 


“glu h (u)du=G, (3) - 
(2.23) [| Цих (8), fa- 224555 


| Clearly, is Laplace transform of f(u). Therefore iheunverson tonne 
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where 6, = Re(z) 


Hence, if the integral equation (2.22) is assumed to have known inversion, 


then h,(x) will be obtained in terms of prescribed function С, (x). Consequently 


the uknown function f(x) will be determined with the help of (2.24). Actually the 


selection of the form (2.11) implies that the resulting integral equation (2.22) is 
standard one. Accordingly we have the following 


Theorem. If f(x) and A, (x) are solutions of the integral equations 


{(t-a,,¢,),, Ё COA 
(6:456), -bifita 


2.25) 809-1, |“ f(u)du- & (x) 


and 


COTON 
(B.E) hi- В,, Е 


(2.26) (х) = |, Цих h, (u)du 


respectively, then 


? 1 гонт xx 
(227) f(u)=lim f eh (z)dz, 
where 
(2.28) G,(x (3) -R[b,0; 1]e( 1/x), 


under the existence conditions of I-function. 
Special Cases. Now, we proceed to the transform of (2.14) in the desired form. We 


assume that any two of HO =1,2,...,p) are equal to 1/2 or reducible to 1/2. For 


convenience we take f, = Б 21/2. Since the 
: constants ц апа Е 
chosen, we put By u, are arbitrarily 


By =1/4+0/2,y, =1/4-v/2., 
Hence the relation (2.13) takes the form 
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Using Bessel’s formula 


(2.30) |, K, (Br)t" dt = 2? виг (u/2+v/2)xT(u/2-v/2), 


we get the following integral 
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provided Re(s)>0,Re(3/2+v+minB,/F,)>0(j=1,2....P). 
Making use of this integral (2.28), we have 


ENTM 


ауди = б, (x), 
(ВЕ), „  ! 5 Bo Fs), s. 


(2.32) [7 f (u) шу)” K, (2uy)1 xy 


where G, (x) 1/x). 


epi 

= Wax”? ( 
Finally the reversal of integration in (2.33) gives the following pair of integral 
equations 


(2.33) In (2ux)? K, (2ux)f (u)du = h, (х), 


(1 E АЕ), 2) : (27:79 (aue s. (Ox 
(85), | Ja та ВЕ; ЕН 


function transform of f(x). Therefore applying 


(2.34) Їл их 


Obviously 1500) is Meijer's Bessel 
the inversion formula of Laplace transform, we get 
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where 6, = Re(z). 
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ABSTRACT 

In this paper we consider certain elementary models of finite but new 
population models of non-interacting species. This population is divided in three 
age groups, namely juvenile, reproductive and non-reproductive. It is assumed 
that the three classes maintain a fixed ratio throughout the period under 
consideration. The length of each regular subinterval is called a generation and 
depends on the particular animal or insect. The formulation is in terms of a system 
of linear finite difference equations incorporating birth rate, death rate and 
migration. The solutions are obtained in terms of finite polynomials having variable 
dependent on the important parameters. Graphs have been plotted for various 
sets of migration rate for different generations. 
2010 Mathematics Subject Classification : Primary 92D25, 33C47 
Key Words and Phrases : Linear Finite Difference Equations, Finite Polynomials, | 
Discrete Population, Migration 3 

1. Introduction. Population dynamics is an important area of 
biomathematics. The population studies are associated with many processes In 
life sciences, such as natural growth and culture of microorganisms, тош апа 
decay of viruses, growth, decay and migration in aquatic and wild life and 
demography. These studies differ from one population species to ше ад 
оп the type of the growth and environment involved. Some populations are isolate 


i i ++ dividing itself or grow out of heterogeneous | 
ks: SUM independently yth асам е Ап eGangotri Initiative " 
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54 
pair contributing to reproduction. While in рет cases two or more species of 
population live together and may be dependent or independent on each other. In all 
cases the environment plays a key role along with other factors such as availability 
of nutrients, presence of toxicants and physic-chemical condition. Thus, population 
dynamics is a subject with many applications. It also plays role m other studies | 
like epidemiology, pollution analysis and ecology. In particular, it plays a major | 
role in the management and manipulation of protected wild life. | 
Any biological population has the following important factors: | 
* Growth/Reproduction 
* Decay/Death 
The Mathematical behaviour has the following types with respect to time 
* Discrete 
5 Continuous 


x Piece-wise continuous | 


Transitional 

The mathematical modeling of each may differ with others in seveal respects. 
This difference may be due to various factors like natality, mortality, fertility, 
community behavior, the impact of environment and genetic effect. | 
Populations are groups of individuals belonging to the same species that | 


live in the same region at the same time. Populations, like individual organisms, 
have unique attributes such as: 


5 growth rate 
age structure 
sex ratio 
mortality rate / 


* 


* 
* 


+ 


. а | 

| if a population A of the same species has finite number of individuals | 
(possibility of ош the sexes and ай the normal age groups) is called a finite | 
ШЕ In this case А is discrete and is governed by difference equation Gotteli | 


A, = А, + f (AnA); 


where A, denotes population after n generations with A, as its initial value is an | 
appropriate function depending on the type of species and environment, 2, is control | 


parameter. If f is linear, say f «aA, | +À, then 


n- ( 

Ac 0 Aj АД ахай +. eor) e o^ AL ХУ oF | 
i=0 | 

шиг Ут мг 


! 
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of approach for quantitative modeling. In many cases we have to consider a given 
all population size as discrete number and its growth as a step function of time. In 
ity such situations the population size is relatively small and the growth rate is 
on generally fixed during a given time interval. This interval we call a step or a 
ies | generation. For all practical purposes we take this interval as a time unit. The 
jor | population under consideration is self-contained in healthy environment and has 


enough amount of nutrition available. There are few population studies in which 
| applications of special functions are classical growth of finite populations Saxena 
[3]. In This paper we make an attempt to express structured finite population in 
terms of polynomials of two variables. These polynomials are close to Hermite 
polynomials of two variables Khan, M.A. et al. 121, Xu 151. 

Formulation of the Problem. We assume that A,, denotes the population size 
after n/^ steps or generations with initial population Ay. Medium and large animal 
species with longer generation of life the population has to be divided in at least 
three groups; we can take Vaishya [4]. 


ts, | A, =A% + А? + AO 40) 
ty, where 

hat | AD — Population of infants and juveniles (pre fertile), 

118, 


AD = Adult population in the fertile age group, 


A = Population of aged non-fertile. 
The change in population of first age group will take due to births, deaths and 
) migration. Therefore. the change in the population A” is given by the equation 


als | ЛАО = AQ — AU. 
2 | The governing equation can be written as 
| MA = B, AC. — D, AD - AP - a QAO, + М,. 
| Accordingly, we get the following difference equation 
an | А9-(1-0,-1,-40)39,-8,А2,4М,. 49) 
ло! | Similarly, we have 
| AA? = AD - AB. 
& The governing difference equation 18 ze 
| A9 = (1 - D, - 7 -a,9) AD, n AD + Ms : 
the | Similarly, we сац х бад Kangri Collection, Haridwar. An eGangotri Initiative rA 
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A9 =(1- D, -а,@) АЎ, +T,A®, +M, (4) 
where B,,D,,T,,M, and q; are respectively birth rate, death rate, transition rate, 
migration rate and predation rate in ith generation. Here we have assumed that 
B, = В, =0 as there is no contribution of birth in second and third age groups. 


Also the entire population is exposed to the predators having by and large a slow 
growing number of predators represented by Q. 
Solution of the Problem. Let 


1-D,-T, -qQ-x. 
Then, equation (2) becomes 

AU =хА® +B AO. +M,- (5) 
Further, let 1- D; -T,-9¢,Q=y. 


Then, equation (3) becomes 


АЁ = yAD, «TL AD. + M,, (6) 
and for 1- D, - q,Q =z, equation (4) becomes 

АР =2A®) VT, AD, +М,. 2 
We find the solution of the equations (5), (6) and (7) by induction method, i.e. 

AD 2 xAD * BAD +M,» ...(8) 
AP = yAD «T AD +М,, .44(9) 
АР = (x? +B,T,)A® + B (x+ y) AP +(1+x)M, +B,M,,, ...(10) 


AP = yAP e T AD + M,. 
Putting the values of A and А, we get 
(2) _ (,,2 
AP co^ * B) AP +T (y) AP (1) M, +T,M,, 
Aj) = xAD BL AP + M,. 


Let us define 


(QD 


р„(х)=1+х+х?+...+х" 
апа 


q(x)e y* y! 4. y^. 
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) Then, putting the values of AU and A® , we get 

3 2 

: AP = 2° AP + BIT, (2х+ y) AP + BAP Y xy + BT AD + M, p, (x) 

г-0 

+B,M, (p, (x) +9, (y))+T,B,M, . (12) 


Similarly, in the equation 
АР = yA) «T AP +M, ; 


putting the values of A® and A® , we get 


) | АР = y AP Bye) AP AD Y y 7 + BT; AD + M, (1+, (y)) 


107 (5), (y))T.M, +B,T,M,. ...(13) 


Likewise 


AD 2 xAD + B A +M,. 


) 


Putting the values of AU and A, we get 


) 
AQ =x! (1) +B 1, А Ў\(3- r)x*"y + BT? A” +В, 0000 TENE 14 ВТ, (2x+2y) AD 
r=0 
) 
y | Map ()* BM, (12x y) * BM; (p (5) * xy +q (y)) + BM, 4 
Again, 
) 
AD = yAP -т,А0-М,. 
Putting the values of AU and А?, we get 
| : m": peu 2 A(2) (1) ху r-1 +2B,T? (x y) AD 
| AD = 54424 BT ADV (4—r)x y + BIT? Ay’ 23, » 11 
) | 4 y А; 1 148 2 | ? 
BTM, ...(1 
| +М,[1+а()]+ BTM, (р (3) 2a (9)) +TM, (р,(х)+ ху + а, (0) + ie 
From 


AP = хА0 + B, AD +M,- 
Putting the values of А! and Af, we get 


(1) B АРУ ху r-l Р 
) xt ty) + ВУ (3х+2у) Ар + / 
1 А ) = зас p +B má Seuk angri Totection. о An eGangotri Initiative=! Ч 


2 
+B?T, (3x7 4xy + зу?) AD + BT? AD + Мур, (х)+ В.Т, (= + Pa(x)+ 2x +q, 07 


Ы 2 
мр) "y +9 (y ум, + В,Т? rez Ty Jus e eret, 


a г=1 


«В| eem ул“ 437116) |м, +В mee m yu, Bere, 


From 
AQ = УА + TA”) + M, , 


putting the values of AU and A”), we get 


3 
AD = y AP + BL ADV (4-r)y x «В(Т, (2x + 3y) AD +T, ух» лур 


г=0 


+B,T? (8x? +4xy +33”) AD + BPT AD + M,[1+4,(y)] 


2 
*B,T, P +1+q,(y)+2xy+ р, (9-1, 
ral 


г=1 rzl 


From 
АО =xA +B AD +М,, 


putting the values of AU and А, we get | 


4 
AP = x AD +B, TAY (5—r)at7 У + ВИТ (6x? + 6xy + 3? АР 


г-0 
5 
t: BT? AD BA) yey ror + BT. {б DH r 27: 
*BjTy (8x + 3y) A + Мур, (x) + BIT? (1+ 3x 2у)М, 


3 1 
+B,T, (È= + py(x)* X (3-r)xy Y" +259 Уу Ju Jm, 
r=) r=0 


r=2 


r=0 


2 
5 )+ ау г Y" Уу хужа, (y )]м,+вутгм, 


2 2 
*B?T, p +р,(х)+4ху+ V ry’ +q, 0))u, 
г=1 A 


r=) 
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-4(16) 


(Т) 


...(18) 


REE a 


М 


) 


8) 


eo 


From 


1 


АР = AQ T AD +M, 


putting the values of А) and А, we get 


4 
A® = y AP + ВАУ ( (5- r) )y* "х" + 3B? 218 ар вэ 2-r Do =| 
r=0 


From 


r=0 


+B°T2 AD + T AC Эрэ Кузгу нэвтр 002+ у) 
+В,Т. 2 AQ |o y cds i y | ма 0) 


BTM.) Ууу + 2y p, (x) )-1] 3x7" + p3(x)- | 


r=0 


+T,M, [ Ps (x)+ х?у+ х? (y) + ха; 57547 (у)]+ BM, (1+2x+3y) 


2 
+B,T? M, P 1)x" + 4xy + 20; +» entre, : 


r=0 


AY = хА® + B, AD + M, 


putting the values of A‘) and AB, we get 


& г,,5-г 2,,3 
A = 2749+ BTA Ў) пету «ey 5 най | 
г-0 


3 A(1) 
+B?T, 12 Б э +BY y «à Baty |+ BIT (4x+3y) 


r=0 


4 2 Gau OP? 
: 720000 0) ry 4-3 x^ y 224 
+B, ADV x*7y' + BT Ао bx y уу 


г=0 


+В?Т2 А? р xi y + 2 + Mp. (x) 


r=0 


: 5 -r A 
+вдм|[ ук+)” edat rtt y +у{рз( 


r=0 
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(19) 


дарш | 


Р, 
! 


D = ss esse 


EE юн a 
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+B,M, | ps) [0s 9) 1) 0)» {Po (х) -1] ху(х° +9") +95 (9) ] 
2270 +1)" +6x7y + 4xy + Oxy? + Уе + D| 
= = 
+B°T?M, [1+ 3x +3y]+ BAT? AD + В?ТЁМ,. ...(20) 
From 


АЎ ур ERAD +, 
putting the values of AU and A”), we get 


(2) _ 7 AQ) (2) ч = x -r r,3-r r, 3- 2 r S 
АГ-уАГВДА, » (6-7) y + BIT? A; р УР i “уу йн "og =i) »"| 


HET AP (edo) AD Ў у Sr 4 Ox 22 ТА UD у +3 "y +х?у i 


r=0 


2 
Bi? AD D» х" у" dt 2 + M, [1 ar 96 (»)] 


r=0 


4 
*BTM, b (r+1)y + Зху?ру (х) + 4xy? + 2xyp, (x)- p, 2 


r=] 


УТЬМ, [as (2) хур, (x) * ху?р, (х)+ y* (р, (9-1 + xy" р,(5)| 


1 тм | х 
| +в? Гг! ву 1+2х+2 2 ; 
| = ( y) |+ BTM, a (y) 42y6xy (x) e 25? + Уух” 
rzl 
© HBITIM, (1+ 3x+3y)+ BT AD + BiT M, . 
| Now, from equation (7) we get E 
| AP -z40 т AD +М,. 
| Therefore, A9 = 240) n AD + M,. 
5 From d 
AQ = „2 A(9) 
p= 2 Al Tie y) AP ANTAN + (1+ 2)M, +т,м, 22) 


— 


AD = 2A + ТАЎ +M, 
Let, . 


T, (z)= le. +2, 
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Then, putting the values of AP and. AU, we get 
2 
Ape 24D LADY ату + BAD VIT (yz) AD +, (М, 
20) +T,M,| r, (z)+ y ]+T,T,M,. ...(24) 
From 
AP = zAP +7,A + M; , 
putting the values of AU and А), we get 
3 
3 AQ = АО T; ADV z^" y' +B, TT, (x 2y  z) AD 
1 r=0 
s| +TT, (5 gi Uy Хул ge sy a? ar ВЕТА + Му» (2) 
г=0 
iT, (a. (9) (2) o2) M, T (m (2) +942) M, + BT TMs. (25) 
From 
АР = zA® +7, AD + М,, 
putting the values of A® and AP, we get 
4 ? 210223 z 
i] {AMA RAMS ery BTA ажан? +249) | 
r=0 
21) 3-r or Ў ху "у" 
| RITE ДӘ nap Ур "x +ye(z+y+z)+ У, 
2 А 
+B,T?T, (z+ 2x 25) AD + Mar (2) € EM [а (0) + 220) e] 
2) | TT. (r.(z) 4 x -2y)M, 
| ST (p. (o) xy + x2 +92+1(2) 1+9 (0). +В, inae | 
| ...(26) 
S» 1 +B, T; T,M 1 
» From 
AD = zAD v T AD + M, 
| 0) and АХ, we get 3 4 
! Putting the values, feda And Kangri Collection, Haridwar. An eGangotri Initiative 7 
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AD = А +Т, ТУУ gy | 
r=0 
калт As iat Hey гу" “УР ry + y?z+ Qxyzt By? + 3y x + Qyx? | 
| +В ТАУ TTAN р ЕН ро 7 ap Door "+ xyz(x y+ 2 
BTA |S zx exzex + (x+y) +2у(х+у+ 2 
r=0 


+B?T T, AD + Mar (2)+Т,М, | ч, (у)+т, (z)+z{1+9, (9) +г°@, (у) +yz’ | 


УЛААН, )+е[у+(@-1+а(у+г)ру) х)+ угт, (z 2+2 ry" + у: d 


ҮВДТ,М,|т,(2)-р,(4)-14(2у +z) p, (x) +2y(2+3y) | 


+B,T;'T,M, (1+ 2x+2y+z)+ В2Т2Т,М, . 221) 
Finally, from 


AP) 2 zAP +T, AD + M,, 
putting the values of A® and A®), we get 


AD = Z1A0 47, AY zy -r,r 


r=0 


реа 


ea O ne 


r=0 


(2) гг 
+B,T,T,As ba x +У( 5-r)y* "Z +xyz(2z+3y+2x De ха ул | 
г=1 


2 
+B TPT, АО р zy +х(22+3у)+ 2) yan зу xy? | 


r=0 
r=0 


эи 


5 -rar ч -r : 
MAD Y a xe У 25у" > у" +x? Y oa 5 yr (y^ + 2) 
= г=0 г=0 г=1 


y *xy 


*B, gem (1) 3-ror 3-r or -r | 
ЗАХ ре х Èz y D х” ES y Hye yay 1222 | 


г=1 


+В At їе +3x+3y) 


+2,M, [n (2) - 122 (у) + (2) a, (y) 2а (0) +? (y) yz (y? 42? ) | 
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n 


da | 


| enm, eo) os em) pi) o) o) Уг 
(+) аг (у) + xem (x) + y2* n. (x x) xy? z |+ B,T,T,M, [z (2) )* z' p. (x) + p, (x) 


2 3 
-1+ хәр, (x) + 2yr (2) + 3^ (z) + 2хуг, (2)+ Y (т +1) у^х” “+B. TM, 


r=1 


Ї г,(2) + 2(x * 9) (z) + зуу" ХЗ 21 B?T?T,M, |r, (z)+2x+3y | 
+B°T°T, AD + BIT; T,M, - .4(28) 
Numerical Examples. Numerical calculations done with the help of Mat lab 7.5 
(2007) programming and also all the graphs drawn with same. 

If g,=0.001,q, =0,D, =20/100,D, =15/100, T,-T,- 


Casel : Examplel. 
T =1/100,B, = 25/100, A® =1000, AS? = 2000, 
21) 260,310,Q = 300, 400,500,600 


--0-300-к0-400--0-500- Q=600 


M, = M, = M =10,60,110,160,210, 


— 
[09] 
e 
e 


210 260 310 
—_> 
AU for Q= 300,400,500,600 


| 50 110 160 
M 


Figure 1: Graph between M and 


1 СС-0. Gurukul Kangri Collection, Haridwar. An eGangotri Initiative 


! 


NT 7 - - "Зэр 


-e Q-300 = Q-400 -+ Q-500 +- Q-600 


64 


10 60 110 160 210 260 310 


Figure-2: Graph between M and AQ for 0-300,400,500,600 


= -8-0-300 = Q-400 =*-Q=500 -$- Q-600 


1 
60 210 260 310 


sc 


| Fi 
| gure-3: Graph between M and A for Q= 
: 4 for Q=300, Q=400, 0-500, Q=600 
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-6-0-300-8-0-400-9-0-500-9-0-600 


10 60 110 160 210 260 310 


M — 


Figure-4: Graph between M and AS” for Q-300, Q-400, Q=500, Q-600 


Тө-0-300 = Q=400 -*- 0-500 -ф- 0=600 


M RRÜÓM M má 4i 


1400 -- I iuc menm 


pM 


en Р 


msc 


i: © d 
iem titan sommes rnnt 
ieee 


Яаг Со eem mn : 


i i 310 
2045 60 110 
—— 
M 
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-@ Q=300 =- =400 -ө--500 -0--600 


Figure-6: Graph between M and AS"? for 0-300, 0-400, Q=500, Q=600 


Extended Hermite Polynomial of two vairable with application in Large 
Size of Animal Population. Now, let z, = Н, (x, y), where 


12 - 
Їл/2| yx 2r j 


Н, (x,y) =n! 2 rin 27) Hermite polynomial. 


n-lz -H,(x,y)-1 
| e y 
| n=2,z, = H, (x, y) - 2! "aca +27 


| 1 
| n-3,2, = H,(x, y) - 3! sala rex 


| È S yx 
| n=4,z,=H,(x,y)=4! Уау testy +1294 


2 
n=5,z, =H,(x,y)= Dara =x" +20x°y + 60xy? 


3 

n=6,2,=H,(x,y)=6!—* __ в 

«7 (x,y) D +30х*у +180x7y? 4.120? 
3 г..17-2г 


= = 9 = Ул : 
п=1,г =H, (xy) тту +4229 + 202 y 4840y, 
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Example Let А; 1000, A = 2000, B, = OD TO 2005 0.01, 
M,- M, = M =10,30,50 


Bliss) BIB) и: Aa ao A) at as 


| 1200 
1000 
-  800- 
AU 600 
400 -|- 

Hs(x,y) 200 
0 

0.8 


Figure-7:Graph between х, y, ACD and Hs(x, y) at М-10, 30, 50 


Example 3. Let А( 21000, A? = 2000,8, = 0.1,7; = T, = T = 0.001, 


M, = M, = M =10,40,70 
` HE Hyxy ШШ А0 atM=10 ШИАО at M=40 BH А at M=70 


(а) =10, 40, 70 
Figure-8:Graph between x, у, 47 and Hz(x, y) at M 


| 
| 
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Conclusion. The.figures drawn depict three significant age groups and 
are self-explanatory. The pupulation under discussion is also under the influence 
of predators whose population does not change so rapidly in comparison to the 
main population of medium size animals. The migration rate is counted as the 
| main factor and governs the growth pattern. 
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In this paper, we present a computational method used oh Cheb 
wavelet for solving a multi order linear and non linear fractiona 
equation. Chebyshev wavelet operational matrix of the fractional integration is 
derived and used to transform the differential equation to a system of algebraic 
equations. The validity and applicability of the technique is demonstrated by some 
examples and comparing the methodology with the known technique shows that 
the present approach is more efficient and more accurate. | 
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Key Words and Phrases : Chebyshev wavelet, fraction 
matirx, multi order linear and non-linear differential equation. 
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ABSTRACT 


AÍ ё б d 1 | D 
GE es 


ifferential 


al calculus, operational 


1. Introduction. In recent years, the use of fractional calculus deals with 
er, in mathematical models, has become 


calculus, as generalizations of classical 
e increasingly used to model some problems 
mechanics [20], dynamics 
s and subtrates [7], non-linear oscillation of 
3], mechanics [27], signal processing [23], control 
ed noise[21]. There has been significant interest 
s for the solution of fractional differential equation. 
ds for solving multi order fractional differential 
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equations as a special kind of fractional differential equations [1,8,9,10,12, 13,14,18, 
19,201. However, few papers have reported applications of wavelets in solving | 
fractional differential equations [26,28]. ! 
Therefore the purpose of this paper to introduce the method of Chebyshey | 
| wavelet to solve multi orders fractional differential equation which include linear 
| and non-linear differential equations. In this paper, the operational matrices of 
| fractional order integrations are derived and a general procedure for forming these 
| matrices is presented. In the proposed procedure, the Chebyshev wavelet expansions P 
| along with operational matrices of fractional order integrations are employed to ( 
reduce the multi order fractional differential equation to systems of non-linear 
algebraic equations. Illustrative examples are given to demonstrate the efficiency 
and applicability of the proposed method. 


2. Basic Definitions and Notations. \ 
2.1. Fractional Calculus. In this section, some necessary definitions and р 
mathematical preliminaries of fractional calculus theory are given, which will be f 


used further in this paper. 


The Riemann-Liodville fractional integral operator I* of order «>0 on the usual 
Lebesgue space L [a,b] is given by [24] as І 


(TA) 0 06-070 211 


: 1 
(г) (8) = (2), (9.1.9) 
where f €C,, 2 -1. 


Its fractional derivative of order «»0) is given by 


| 123285: (12255107 for n-l<a<n, 


| (2.1.8) : 
where n is an integer and feG. 

1 Now Бу the Riemann-Liouvilles definition | 

| Tet ={(Г(ь+1))/(г(а+ь+1))}+ | 

| ITF (t) =I" F(t) ( 

| I'Pf(t)- ГІ“) : сис 


where o, 2 0, >0 and [> 
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18, There MR MR. a Riemann-Liouvilles derivatives when trying to 
ng model rea. AE p pau gih fractional differential equations. Therefore 
we need to introduce a modified fractional differential operator proposed by Caputo 
lev [25] 
"ar a far 1 n-a-l „л 
of D IUe rcd) [(t-2) f” (z)dz,(n-1<a <n) (2.1.5) 
25е 
ns where n is an integer />0 and fe Ce 
to Caputo integral operator has a useful property, 
заг n-1 w t 
Су 1-1) Fe)=F@)- dF (0 ),(@-1<«<л), (2.1.6) 
where n is an integer £20 and ѓєС". 
nd . 2.2. Chebyshev Wavelet. Wavelets constitute a family of functions constructed 
be from dilations and translations of a single function called the mother wavelet 
al w(t). We have the following family of continuous wavelets in case the dilation 
1 : 
parameter a and the translation parameter b vary continuously, 
a (t-b 
wal) 5422) aber, ono. (2.2.1) 
.1) a 
If we suppose parameter a and b to take discrete values as ага, 0-1, ‚а, » L, 
2 b, >0 forn апа positive integers, we have the following family of discrete wavelets 
(2.2.2) 
Vaa (£) = w (ait нь) 
Where w,,(¢) forms a wavelet basis for 12(В). In particular, Wy, (t) forms an 
3) orthonormal basis, if 00:22, 5971, that is 


(Van (t) QUA (9) = быбл- 


The Chebyshev wavelets ү, „ 
п=1,2,...,9-1 and ñ=2n-—1. m is the order of the Chebys 


(£) 2 w (E fim,t) have four arguments; keN 
hev polynomials of the 


first kind and £ is the normalized time. 


A) Chebyshev wavelets w,,,,(¢) are defined on th 


e interval [0,1] as 


Wee BCE ep op (2.2.3) 
= otherwise ? 
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1s, m=0 | 
| where 7, Т„(@) = 6, 750 
and m-01,.,M-1n-12,.,2 where k is any positive integer and T, (t) are | 
Chebyshev polynomial of the first kind of degree т, which are orthogonal with 


; 1 
respect to the weighted function e(t)- ШОШ z n the interval [-1,1] and Т, (t) can 


be determined by the following recurrence formulae, 


T, (t)=1,7,(t)=4,T,,,, (2) - 27, (1) - T. (2), т=1,2,3,... 


Now in dealing with the Chebyshev polynomials the weighted function @ = o(2t - 1) 
has to be dialated and translated as 

х: o, (t)=0(2'¢-A), (2.2.4) 
to get orthogonal wavelets. 
А function f(/) defined over [0,1) may be expanded as follows 


= Y Y C, ws (t), T (2.2.5) 


n=0 mez 


by the Chebyshev wavelets, where 


25101710) = 210121012077 (2.2.6) 


| 
In which (5,5) denote the inner product in Lo, [0,1]. If the ifininte series in equation 


(2.2.5) is truncated, then it can be written as 


f(t)s SSc amam (t) = C™¥(t), 


n=l m=0 


(2.2.7) 
where C and (t) are ока M x1 matrices given by | 


C= sese. „С 


| 
| шан (2.2.8) | 
| апа E 
| 
| КОСЕ (2.2.9) | 
Taking the collocation points as | 7 
"m T TY 
i ar? 1-12,.,2- M, (2.2.10) | 


CC-0. Gurukul Kangri Collection, Haridwar. An eGangotri Initiative 


Нил : 


vith 


we define Chebyshev wavelet matrix 9... 


as, 


1 3 2m-1 
= ур —— —— 
Omm | (85) (2: Е | (2.2.11) 


where m= 2 М- 


For example, when M=3 and K=2, the Chebyshev wavelet is expressed as, 


2.2568 2.2568 2.2568 0 0 0 
1.0688 9.5746 18.0854 0 0 0 
-2.4828 54.2562 201.7761 0 0 0 

Цэн 0 0 2.2568 2.2568 2.2568 
0 0 0 1.0688 9.5746 18.0854 
0 0 0 -2.4823 54.2562 201.7761 


A function k(x,t)e L, ([0,1][0,1]) may be approximated as 

k(x,t) = (х) КУ(0), (2.2.12) 
where K is 2! M x2** M matrix with 
2.3 Operational matrix of the fractional integration. The integration of the 
vector ¥(ż) defined in (2.2.9) can be obtained as [vat = P'¥(t) where P is the 


2^ М х9 ү operational matrix for integration [17]. 

Our purpose is to derive the Chebyshev wavelet operational matrix of the 
fractional integration. We rewrite Riemann-Liouville fractional integration as 
follows . 


(7710) = zs 16-97 farm entro зл 


: : ti d 
Where œe R is the order of the integration, T (a) is the Gamma function ап 


t" * f(t) denotes the convolution product of ;** and f). 


Now if is expanded in the Chebyshev wavelet, as shown in the equation (11) 


the Riemann-Liouville fractional integration becoems, 
e а-1 [3 x P(t). (2.3.2) 
(rf f)(t)= E Fay’ «f(t)*= = CT T(a ig 4 ) 


lets, the 
Thus if у=, f () can be integrated, then expanded in Chebyshev wavelets 


hev wavelets. 
Riemann-Liouville fractional integration is solved via the Chebyshev 


Next, we define a $6-8eufcBloelkqvalee; нена (RE PIS Initiative 
л ГҮКҮ 


$ 


14 | чау 
31, i/ms<t<(i+1)/m 
| (Ð= 
50) js otherwise й (2.3.3) 
where 1=0),1,2,...,(т-1). 
The ВРЕ has two useful properties which will be used further, 
1 
| 
| 
| 
| 


(0) Disjointness. | 
0, 121 à 
ИО is E оза | 
(ii) . Orthogonality. Y 
0, EA | 
[nhe T (2.8.5) | 
Similarly Chebyshev wavelet may be expanded into an m-term BPF аз. 
Wm (2) = Omm B, (t), (2.3.6) 


where B, (7) [5 (!),^ (t)... (01. Ў | 


Kilicman and Al Zhour [18] have given the Block Pulse operational matrix 
of the fractional integration Е° as Ps 


(IB, (t) « F*B, (t), | | 


where 


NS ЧАЙ (2.3.8) 


With. & A E | 
зогсоо the Chebyshev wavelet operational matrix of the fractional | 
Let (I°Y,, )(t)= Ре ле (t), (2.3.9) 
where matrix Рэліз called the Chebyshev wavelet operational matrix of the | 
fractional integration. Now on using equations (2.3.6) and (2.3.7), we have | 
(^v, (t) = (178... B, (t) =, (I*B, (t) = bnn Pp. (0). (2.3.10) | 


From equation (2.3.9) and (2.3.10), we get 
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3.3) 


3.4) 


3.5) 


rix 


PETS (t) = Bose В, (£) = QUIE (4). (2.3.11) 
Then, the Chebyshev wavelet operational matrix of the fractional integration p* 
is given by бэр 
D 2: po AES 


msm“ (2.3.12) 
Therefore, we have found the operational matrix of fractional integration for the 


Chebyshev wavelet Р“ 


The fractional integration of the function / was selecteds to verify the correctness 
of matrices. : 
That is because the fractional integration of the function f()-t is easily obtain as 


a Г2 А 
follows, (1 f)(t) = TE ! which is easily used to compare the result obtained 


by the proposed method. 

3. Applications and Results. In this section, we will use the Chebyshev 
wavelet operational wavelet operational matrices of the fractional integration to 
solve linear and non-linear fractional order differential equation. These examples 
are considered because closed form solutions are available for them of they have 
also been solved using other numerical schemes. This allows us to compare the 
results obtained using this scheme with the analytical solution or the solutions 
obtained using other schemes. 

Example 1. Our first example deals with the homogeneous linear fractional 
differential equation 


D*y(t) = -y(t),y(0) 21,0 «2 «1 where 0<a<1 (3.1.1) 
Let, D*y(t) = KEY, (9, 
y(t) = Кера (4) ur [ШЖ B, () 
Since E (0) = Omxm Bm (2), 
y(t) m КО pu DA (2) + [11,...1В, (t). 
Therefore by equation (3.1.1) 
(3.1.2) 


каз (EIAS 05-85 (¢)+ [1,1,...,1]В„ (¢)=0. | 
This is a system of linear algebraic equation. ув оҳе MATLAB to solve equation 
(3.1.2). The exact solution of equation (3.1.1) is given by 


y(t) = E, (-t*), 
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С“ -_[.’__ҸМ и 


у z! i 1 р 1 f order 
б = y ———— js the Mittag-Leffler function о a. 
where E, (2) LTR) 


| | 4 
Numerical results for example 1 when о = 0.5,0.75 and 1 with comparison to ref. 
| [22] are given in the table (1). 

| Table 1 

| 


25 
y(Appro.) y(Appro.) y(Appro.) 
1.000000 1.000000 | 1.000000 1.000000 1.000000 
0.723854 0.723578 | 0.828250 0.828251 0.904837 
0.644794 0.643788 | 0.732486 0.732585 0.818732 


y(exact) 
1.000000 
0.904837 
0.818731 


0.592179 0.592018 | 0.660340 0.660337 0.740816 | 0.740818 
0.553728 0.553606 | 0.602023 0.602121 0.670319 | 0.670320 
0.523256 0.523157 | 0.553606 0.553603 0.606529 | 0.606531 
0.498126 0.498025 | 0.512289 0.512285 0.548811 | 0.548812 
0.476605 0.476703 | 0.476557 0.476555 0.496586 | 0.496585 
0.458149 0.458246 | 0.445294 0.445292 0.449328 | 0.449329 
0.441954 0.442021 | 0.417684 0.417682 0.406570 | 0.406570 


0.427611 0.427584 


0.393106 0.393108 0.367079 | 0.367079 
From table 1, it is clear that the approximate solutions of the equation (3.1.1) 
in high agreement with the exact solutions for different values of a. 


Example 2. The second example deals with t 
equation 


D*y(t)=[y(t)} - F(2) 


are 


he nonlinear fractional differential 


(3.2.1) 
where 


2 
(1) “үгээр”? ын ~2,t > 0,0 «ас 1. 
Let 


D*y(t)= KT, (1) 


then y(t) - КТРЗ nY, ()-2[1,1,...,1]B,, (2). 
Since Y, (2) = Ф „В, (t), 

y(t) = КЕРЕН. (t) = 2[1,1,...,1]B,, (t) 
Let 
KT p 


m Tue =[а,,а„,...,а„] апа using equation (2.3.3), we ауе, 
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Ч шиг - 


NE. =| a,b, (0) сад, (t) a,b, (От 18В, (t) 
fue н ‚в ]B, (t) )-2[1,1,...,1]B,, (1) 
Now f(t) = 7+, (£) and +, (г) =, „В, (t) 


тхт т 


So f (t) = fn mmB, (t). Hence by equation (3.2. 1), we have 


KaYa (0) [het ] B, (0) + 2[11....1]B, (t) лт, „В, (1)-0 (3.2.2) 
This is nonlinear system of algebraic equations, here we use MATLAB to solve 
equation (3.2.2). The exact solution of equation (3.2.1) in case of a=1 is given by 


y(t) =-2/(t +1) (3.2.3) 

Table 2 shows the approximate solutions of the equation (3.2.1) for different 

values of « using the method of Chebshev wavelet operational matrix. The value 

of 0-1 is the only case for which we know that the exact solution is y(t)=-2/(t+1) 

and our approximate solution of the equation (3.2.1) in this case is in high agreement 
with the exact solution. 


Table 2. Numerical results for example 2 when a =0.2,0.4,0.6,0.8 and 1.0 


-2.00000 


-2.00000 -2.00000 -2.00000 | -2.00000 


-1.57189 -1.58883 -1.56450 | -1.73540 -1.81817 
-1.48649 -1.48992 -1.53360 | -1.59188 -1.66666 
-1.42058 -1.41462 -1.44370 | -1.48260 -1.53845 
-1.36657 -1.35255 -1.86907 | -1.39286 -1.42855 
-1.32105 -1.29978 -1.30710 | -1.31635 -1.33333 
-1.28225 -1.25414 -1.25245 | -1.24985 -1.25000 
-1.24859 -1.21416 -1.20425 | -1.19121 -1.17649 
-1.21925 -1.17881 -1.16135 | -1.13896 -1.11111 
-1.19340 -1.14720 -1.12280 | -1.09200 -1.05264 
-1.170478 -1.11891 -1.08801 | -1.04955 -0.99999 


Example 3. Our third example covers the inhomogeneous linear equation 


2 Yo E 3.3.1) 
pre 1275-51-44, ( 
901» (0* «5-5 res) 

with у(0)-0,/»0 and 0 «a «1- 


2 ze рер 
"mee. e) 
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Шилэн 


EEUU aü—. 


m 


Then y(t)= KARA ES |: ши )= mx m B, (t ). 


Hence y(t) = K p Pim mm Bn (0) S 
Now f (t) = 77%, (t) and V, (t) = mm Bn (0) 
Therefore f(t)- f74,,.,, B, (t). (3.8.4) 


where f7 is a known constant vector. 


Now making an appeal to (3.3.2),(3.3.3) and (3.3.4) from drive (3.3.1), we derive 


Ki m(t )+ Kp. „Ф. „В„(ї )- ihe ол B, (t )= 0. (3. 3. 5) 
This is a nonlinear system of algebraic equations, here we use the MATLAB 
function to solve the equation (3.3.5) 
Th exact solution of equation (3.3.1) is given by y(t E 2-1 

The linear equation(3.3.1) is solved by K. Diethelm using the fractional 
Adams-Bashforth-Moulton method [11]. Table 3 shows the exact solution by [11] 


and the approximate solution by our proposed method of equation (3.3.1) for 


different values of «. The results compare well with those obtained by Adams- 
Bashforth-Moulton method. 


Table 3. Numerical results for example 3 when «=0.5,0.75 and 1.0. 


a 
Es tS y (exact) 

0 | 0.000000 | 0.000000 | 0. 000000 .0.000000 | 0.000000 | 0.000000 
-0.089820 | -0.090000| -0.089894 | -0.090000 -0.089925 | -0.090000 
-0.159791 | -0.160000| -0.159856 | 0. 160000 | -0.159956 | -0.160000 
-0.209818 | -0.210000| -0.209876 | -0, 210000 | -0.209965 | -0.210000 
-0.239851 | -0.240000| -0.239913 | 0 240000 | -0.239976 | -0.240000 
“0.249894 | -0.250000| -0.249917 | 0 250000 | -0.249974 | -0.250000 
-0.239963 | -0.240000| -0.239972 | 9 240000 | -0.289986 | -0.240000 
0.210026 | -0.210000| -0.210015 | о. 210000 | .0.210012 | -0.210000 
-0.160089 | -0.160000| -0.160066 | о 160000 | -0.160022 | -0.160000 
-0.090175 | -0:090000| -0.090102 


-0.090000 


th 
сора он Chebyshev Wavelets expanaioas &ogethapri 


| 
T 


р roposed to obtain the numerical solutions of mul 
equations. The approximate solutions for some fr 
are obtained by our method. The proposed method i 
as it is computer oriented. Several examples a 


ti order fractional differential 
actional differential equations 
5 very simple in implementation 
re given to demonsdtrate the 
powerfulness of the proposed method and numerical results show that the method 
is very efficient and correct for the numerical solutions of linear and nonlinear 
fractional differential equations. 
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ABSTRACT 

Almost hundred years ago research in Summability Theory and Applications 
started in India (precisely, in South India). Slowly but surely, it started to attract 
the scholars to work in this field. For more than sixty years, this fascinating field 
of research remained popular among the scholars. The first aim of this paper is to 
introduce the Summability Schools in India, mostly in Central and North India, 
scholars and research topics on which researches were carried out their work at 
thesecentres of research. Second aim is to let the readers know the migration of 
the scholars from one centre to another mostly to establish a new centre for 
research, 
2010 Mathematics Subject Classification : 01-02. 
Key Words and Phrases : Survey article, Summability Schools in India. 

Introduction. The present paper (survey article) is based on my two 
lectures. The first lecture, which I gave in the “National Conference on History of 


ns 


А, Mathematics” organized by Department of Mathematics, Kumaun University, 

Nainital(Uttarakhand) during October 13-16, 2000 and the second “U.N.Singh 
al Memorial Lecture” was delivered in the “International Colloquim on Historyof 
2) Mathematical Sciences and Symposium on nonlinear analysis” organized by Indian 


Society for History of Mathematics at SSJ Campus, Almora, [шап University 

Nainital, during May 16-19, 2011. The whole paper is divided into the following 

sections: к ditv Th : 

1. General Information regarding early researches in Summability Theory in 
India. 

2. The first Summability Sc j Е 

3. Other Summability Schools and Scholars in the Nor 


me i igrations. 
India in Summability Theory and their migra ; d 
l. General Informations regarding earlyresearches in Summability 


Theory in India. By Summability School, we mean a Doe Bi € m 
of a University or Institute or even College where the p es 2 rd d 
theory and its applications have dominated over the researc 5 ? цэн Бэ 
branches of Mathematics. Through publications, Т could locate the fo g 
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hool in the Central and North part of India. 
th and Central part of 


или 0077 


where the researches in summability theory were carried out. 


PLACES IN INDIA 

1. Aligarh 2. Allahabad 

3. Annamalainagar 4. Baroda(Vadodara) 
5. Berhampur 6. Bhubaneswar 
7. Bombay(Mumbai) 8. Cuttack 

9. Gorakhpur 10. Jabalpur 
11.Kanpur 12.Kurukshetra 
13.Madras (Chennai) 14. Madurai 
15.Rohtak 16.Roorkee 
17.Sagar 18.Sambalpur 
19. Ujjain 20.Varanasi 
21.Waltair. 


From publications, it appears that researches in summabilitytheory and 
applications in India particularly in South India are being done since long, precisely, 
not later than 1918, which is clear from the following publication of a senior most 
worker of this field (late) K. Ananda Rau, Professor of Mathematics, Madras 
Presidency College, Madras: 

Proc. London Math. Soc.(2), 17(1918), 334-336. 

T. Vijayraghavan andC.T.Rajagopal were among those pupils of K. Ananda 
Rau who were very much inspired by him. T. Vijayraghavan, who 
influencedRajagopal in his later career, was a senior college-mate in the Honours 
classes of Rajagopal. Rajagopal began his teaching career in the Annamalai 
University in 1930-31. Those days, a very few universities were having their own 
teaching departments especially in Mathematics, and therefore it was very much 
difficult to have a team work. It was 1950, when Ramanujan Institute of 
Mathematics was founded by AlagappaChettiar, a well-known educationist of that 
time, in the nem Madras University, Madras. In 1951, T. Vijayraghavan took the 
charge of headship and later C.T.Rajagopal also joined the Ramanujan Institute of 
Mathematica] After that other people also got associated with Ramanujan Institute 
л Mathematics(RIM). Before we give a few names who worked at Madras, Waltair, 
sine ee Oba in ТЕЕ 

е d que contribution in *RIESZ SUMMABILITY 
which is due to K. Chandrasekharan and S. Minaksh 
published for the Tata Institute of FunndamentalRe 
Oxford University Press in 1952. : 
NAMES OF SOME SCHOLARS 
1. K.Ananda Rau, Madras 


2.  T.Vijayraghavan, RIM, Madras (1918) 
СС-0. Gurukul Kangri Collection, Haridwar. An eGangotri Initiative (1926) 


isundaram, a monograph 


earch(TIFR), Bombay, bY | 


г —ии = : 
5) 2 


V. Ganapathylyer, Madurai 
V. Ramaswami, Madurai 
C.T. Rajagopal, RIM, Madras 


(1935) 
(1935) 


3 
4 

: i i (1935) 
6, S. Minakshisundaram, Andhra University, Waltair (1936) 
7 

8 

9 


К. Chandrasekharan, TIFR, Bombay (1942) 
K.S.K. Iyengar, RIM, Madras (1943) 
M.S.Ramanujan, RIM, Madras (1956) 
10. M.R.Parameswaram, Madurai (1957) 
11. M.S.Rangachari, RIM, Madras (1961) 
Note: Against each name, number indicates one of the years of his publications, 
not necessarily the first one. 
2. The First Summability School in the Central and North part of India 
(Allahabad). It was July 1932, when B.N.Prasad (12 Jan 1899-18 Jan 1966) returned 
back to Allahabad from Europe after earning his degrees of Ph.D. and D.Sc. 
ly respectively from Liverpool(England) and Paris(France). He started building up 
at Allahabad a strong school of Research on Summability Theory in which he 
succeeded. In the North and Central part of India, it was the first Summability 
School. 

Earlier, B.N.Prasad completed his M.Sc. in 1921 from Banaras Hindu 
University. On the advice of Ganesh Prasad, he joined the Department of 
Mathematics of Banaras Hindu University as a research scholar in July 1921 
with Ganesh Prasad and later, after one year, he joined the same Department as 
Assistant Professor of Mathematics. But after the departure of Ganesh Prasad 
from Banaras Hindu Universityto Calcutta he could not continue at Banaras Hindu 

Universityfor more time and left it to join the Department of Mathematics of 


Allahabad Universityin 1924. And in 1929, he left Allahabad for England for 
heory of trigonometric series by 


intensive study of the works done in the t 


W.H.Young, G.H.Hardy, J.E.Littlewood, E.C.Titchmarsh, H.Lebesgue, etc.. He 


Worked with E.C.Titchmarshfor Ph.D. degree at Liverpool and later with celebrated 


French Mathematician Arnaud Denjoy for D.Sc.degree at Paris, which was awarded 
оп the 4th June 1932 by the University of Paris. 
After establishing a School of Summabi 


founded the Allahabad Mathematical Society in 1958 an 


Journal of Mathematics” to facilitate the workers for the research opportunities 
ned an international status very soon. 


lity at Allahabad, B.N.Prasad 
d began publishing “Indian 


at this center of learning. This journal attai 


2.1. Scholars of the first generation of t 
f some scholars o 


Prasad, that is, pupils of B.N.Prasad: 


his School (Pupils of B.N.Prasad) 


ration of 
The following are the names о f the first gene 


theSummability School headed by B-N- 
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M.L.Mishra (retired from Sagar, M.P) 

PL. Bhatnagar (retired from Shimla, Н.Р) | 

U.N.Singh (1917-1989) (retired from University of Delhi) 

J.A.Siddiqi (Canada) 

.M.Mazhar (Kuwait) 

idm ma from Allahabad Universityand expired in 2008 July 02) 

S.R.Sinha (1931-1985) (expired while in the service of Allahabad University) 

D.PGupta (retired from M. L. N. Engg. College (now NIT), Allahabad; currently 

President of Allahabad Mathematical Society and is settled in Allahabad) 

9. S.N.Bhatt (retired from Allahabad University; he is settled in Allahabad) 

10. Tarkeshwar Singh (retired from M. L. N. Engg. College (now NIT; he is settled 
in Allahabad), Allahabad) 

11, N.D.Mehrotra (retired, Allahabad) 

12. Pramila Srivastava (retired from Allahabad University; she expired in 2008 
April 12) 

18. L.M.Tripathi (retired from Banaras Hindu University) 

14. SulaxanaKumari Gupta (Gorakhpur) 

15. G.D.Dixit (retired, Auckland, New Zealand). 

Most of the first generation workers(pupils of B.N.Prasad) of this school 
migrated either to take active part in establishing the learning centres in various 
universities of India or abroad as indicated in the brackets of the above list. I will 
give some brief commentaries about some of these scholars. Presently, I give some 
observations about PL.Bhatnagar and U.N.Singh. 

PL.Bhatnagar who turned out to be an applied Mathematician was primarily 


associated with this School of Summability established by B.N.Prasad. Under 

Prasad's supervision, he published the followin 

1. On the convergence of the conjugate series 
Nat. Acad. Sci. India, 7(1937), 116-122. 

2. On the summability of the conju 
Sci. India, 7(1938), 85-90, 


8 papers on summability: 
of the derived Fourier series, Proc. 


gate series of Fourier series, Proc. Ind. Acad: 


The absolute summability of the conjugate series of the derived Fourier series, 
Bull.Calcutta Math. Sco., 30(1938), 17-26. 


4. Onthe Abel summability of the co 
Math. Soc., 81(1989), 31-44, 
However, PL.Bhatnagar found Pl 
SummabilityTheory so he switched over to Ast i 
roph ear 
worked with A.C.Banerji for his И cem 


(i) D.Phil. (Allahabad University, 1940) 


njugate series of Fourier series, Bull. Calcuti? 
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(ii) D.Sc. (Allahabad University, 1974). 
UNSingh, after completing Post-Graduation from Allahabad University 
in 1944, he joined ш Summability School and completed his thesis entitled “Strong 
Summability of Trigonometric Series” and was awarded D.Phil. degree in 1949, 
Later in 1951, he went to Paris to work with famous Franch Mathematician S. 
Mandelbrojt. U.N. Singh was awarded D.Sc. degree in February 1954. After coming 
back to India, he changed his interest to Modern Mathematics. 
2.2. Scholars of the Second and the Third Generationof this school: 
1. §.N.Lal (retired from Banaras Hindu University): worked with S.N.Bhatt. 
He is settled in Varanasi. 
2. Nand Kishore (retired from Berhampur University, Odisha): worked with 
S.R.Sinha. 
39.  ZU.Ahmad (retired from Aligarh Muslim University): worked with T.Pati 
008 . for (i) D.Phil. at Allahabad and for (ii) D.Sc. at Jabalpur. 
4. НР Dikshit (retired from Jabalpur University): worked with T.Pati for 
(i) D.Phil. at Allahabad and for (ii) D.Sc. at Jabalpur. 
5.  J.BSharma: worked with D.PGupta, MLNR Engg.College(now NIT), 
Allahabad. 
6. Lal Bahadur Singh: worked with S.N.Bhatt. 
7. Vikramaditya Singh (retired from HRI, Allahabad): worked with S.R.Sinha 
for D.Phil. &D.Sc. at Allahabad. He is now settled in Allahabad. 


1001 
ous 
will 


me | 8  UshaKakkar: worked with S.N.Bhatt. 

, 9. 5.5. Tripathi: worked with S.N.Bhatt. 

у | 10 Anand Prakash: worked with Pramila Srivastava. 
der 


11. Ashok Saxena: worked with S.R.Sinha. 
12. K. Prasad: worked with S.R.Sinha. 
13. Shyam Lal: worked with L.M.Tripat 
University in 2001 and stayed from 
period he supervised (i)Prem Narayan Sin 
Kumar Tripathi (iv) Anjali Shrivastava (v) Ramaa 
Yadav. >: 
2.3. Topics of Research on which researches were carried ou 
(Allahabad) are: 
(A)Dirichlet’s Series (D.S.) 
о Multiplication of Dirichlet series. ds a 
ti) Multiplication of absolutely summable 455. : 
(ili) The РЕ. of the Dirichlet’s product of ee E -— 
(V) Тһе abscissa of (a) Strong summability of D.S. аар 
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hi at B.H.U. and joined Allahabad 

28-04-2001 to 17-05-2004. During this 
gh (ii) Gopal Krishan Singh (iii) Vijay 
kant Yadav (vi)Ramashrya Singh 
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(v) 


Strong RieszSummability of D.S. 


(B) Infinite Series 


q) 


(1) 


(iii) 


(iv) 


(v) 


Limitation Theorems for Absolute RieszSummability Method. 
Tauberian Theorems for 
(a) Absolute NórlundSummability Method 
(b) Absolute Matrix Summability Method 
(c) Absolute RieszSummability Method 
Consistency Theorems: 
(a) Consistency Theorems for RieszSummability 
(b) The Second Theorem of consistency for the absolute RieszSummability 
(c) Unified Theorem for consistency for the absolute Rieszsummability 
(d) Consistency theorems for absolute Rieszsummability with indices. 
Summability of Infinite series by 
(a) Absolute Abel Summability 
(b) Absolute CesàroSummability 
(c) Absolute NorlundSummability 
(d) Absolute RieszSummability 
(e) Absolute RieszSummability with indices 
(f) Strong RieszSummability 
(g) Discontinuous Riesz means. 


Absolute Summability factors for Cesaro, Riesz and Nérlundsummability 
methods. 


(C) Fourier series and allied series 


(i) 


(ii) 


(iii) 


Summability of Fourier series by 


1 
(a) (n=) -and (b) (N, p,)— mean. 


Absolute Summability of Fourier series b 
RieszSummability methods. 


Localization Problems for Absolute Sum 
factored Fourier series for 


(a) |A|-summabilitymethod 


y Abel, Cesàro, Nórlund and 


mability of Fourier series and 


(b) IC, k|(k > 0) summabilitymethod 
(с) IN, p,|- method 


(d — |R,X,,o|(0.» 0) method. 


(D) Some Other Areas 


G) 


Gibbs Phenomen б 
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() ^ Convergence of Ultrasphericallacunary series 
(iii) Local property of Absolute Summability of Laplace series and rth derived 
series of Fourier series 
(iy) | Summability and strong RieszSummability of orthogonal series 
() Absolute NórlundSummability of Jacobi series 
(vi) Strong Riesz means of orthogonal series 
(vii) Bounded variation of power series. 
3. Other Summability Schools and Scholars in the North and Central part 
of India and their migrations : 
3.1. ALIGARH 
J.A.Siddigi, U.N.Singh and S.M.Mazhar were among the main workers in 
the theory of summability and applications in early days in the Department of 
Mathematics, Aligarh Muslim University (A.M.U.), Aligarh. They were migrated 
from the Summability School at Allahabad. These scholars could not stay more 
time at Aligarh but during their stay they trained a number of persons in this 
field. In 1958, U.N.Singh left Aligarh for Baroda as Professor of Mathematics to 
establish a new Department of Mathematics in M.S. Rao University of Baroda. 
J.A.Siddiqi left for Canada and settled there, he is now late. S.M. Mazhar also left 
Aligarh and settled in Kuwait, he is also late. In the meantime 7. U.Ahmad, a senior 
most pupil of T.Pati, joined Department of Mathematics of A.M.U., Aligarh. 
There were a number of persons who worked in Summability Theory and 
Applications at this centre of research . Names of some of the scholars who worked 


and still working are: 
J.A.Siddiqi(retired & expired in Canada) 
U.N.Singh(retired from Univrsity of Delhi) 
S.M.Mazhar(retired in Kuwait) 

Sarfaraz Umar(retired from A.M.U., Aligarh) 
Z.U.Ahmad (retired fromA.M.U. , Aligarh) 
Niranjan Singh (retired from Kurukshetra) 
А.Н. Siddiqi (retired from A M.U, Aligarh and 
Abrar Ahmad Khan 

EM. Khan (retired from A.M.U., Aligarh) 

10.  PM.Abdul Rahiman 

П. МА. Khan 

12. Q.A.Khan 

13. АН.Ѕаій 

М. K.C.Varshney(retired from A. 


15. HH Khan | 
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17. RS Khan 
18. H.TJalal 
19. —Mursaleen. 


3.2. VARANASI ; > 
Pramila Srivastava and L.M. Tripathi were pupils of B.N. Prasad. They 


left Allahabad to join the Department of Mathematics of Banaras Hindu University, 
Varanasi. Later, S.N. Lal, who worked with S.N. Bhatt, alse left шиш to Join 
the Department of Mathematics, Banaras Hindu University, VERE They trained 
a number of students in Summability Theory. However, Pramila Srivastava came 
back to Allahabad in 1970 and retired from the Allahabad University, Allahabad, 
while L.M. Tripathi and S.N.Lal, continued to be in B.H.U., Varanasi, and retired 
from there as Professor of Mathematics. We now present a list of scholars (may be 
incomplete) in Summability Theory who worked in the Department of Mathematics 


of Banaras Hindu University, Varanasi. 

(a) The pupils of L.M.Tripathi are: 

1. G.C.Dvivedi 2. Kanhiya Prasad 3. A.PSingh 

4. K.N.Singh 5. Vaagesh Narayan Tripathi 6. VP. Yadav 

1. Shyam Lal 8. Shyam Narayan. 

(b) The pupils of S.N.Lal are: 

9. Sheo Ram Singh 10. Raj Deo Ram 

11. Syaram 12. Kailash Nath Singh. 

1.2.1. There had been another centre for research in Summability Theory located 
in VARANASI “Harishchandra PG. College, Varanasi”now affiliated with 
Mahatma Gandhi KashiVidyapith, Varanasi. 

| In 1993, Shyam Lal of B.H.U. joined the Department of Mathematics of this 


college and supervised the following pupils in Summability Theory: 
1. Sunita Verma 


3. K.N.Singh Yadav 


In November, 2007, Shyam Lal rejoined the Department of Mathematics of 
B.H.U. as Professor of Mathematics and supervised the following : 

1. Dr. Binod Prasad Dhakkal 2. Dr. H : Si 
S icu a - Harendra Prasad Singh. 


2. Ajay Pratap 
4. Surendra Singh Yadav. 


Т ВЭБ S 
5 
Un 
= 
3 
Б 
er 
- 
гэж 
Ф 
o 

d 
О 
3 
Ф 
o 
к 
2 
Ф 
м 
E 
Ч 
Е 


upta after her marriage to a Doct 


resigned from the учеа нр њ B RM teles вав gane from England 


— ұы ӘӘ 


| ce «ue University During his period, B.P Mishra, who expired in 2009 June 


01, supervised the following students in Summability Theory at Gorakhpur: ' 


1. Dinesh Singh 7. Anand Prakash Shrivastava 
2, Sanjay Kumar Singh 8. Eklakh Ahmad Khan 
oy 3. Anupama Shrivastava 9. Parthiweshwarji Pandey 
y, 4. Mamta Rani Jaiswal 10. Suyash Narayan Misra 
In 5. Sachchidanand Pathak 11. Harish Chandra Singh 
əd 6. Gunjeshwar Shukla 12. Dig Vijay Nath Pandey. 
ne 3.4. KANPUR 
d, Kanpur is not much known for the workers in the Theory of Summability 
ed and its Applications. However, R.S.L. Shrivastava, Department of Mathematics, 
be LI.T, Kanpur contributed here to the Theory of Summability and its Application. 
cs He was associated with B.N.Prasad of Allahabad. 
The other place in Kanpur, which I could locate is H.B.T.I. In fact, Shyam 
Lal, who was from Harishchandra PG.College, Varanasi, joined H.B.T.I., Kanpur 
in August 1999 and stayed there till April 2001. He supervised the following 
students for their Ph.D. degree of Chhatrapati Sahu University, Kanpur: 
1. Hare Krishan Nigam 2.Virendra Nath Tripathi 
3. Purnima Yadav. 
3.5. ROORKEE 
In the Department of Mathematics of the University of Roorkee(now L1.T.), 
ed the work in Summability Theory and Applications were initiated by O.PVarshney. 
th He worked in Summability Theory with M.L.Mishra (perhaps the first pupil of 
B.N.Prasad of Allahabad), Department of Mathematics, Saugar University(now 
is Dr. H.S.GourUniversity), Sagar and left this university quite early to join the 
University of Roorkee. In September, 2001, this university was declared 
LI.T.Roorkee because of the institution of National importance. O.P. 5 
worked there andeducated some scholars in Summability Theory and Applications 
of Those who have carried out or are still carrying their research work Gu this 
institution) in Summability Theory and Applications are the following: 
O.P Varshney (Roorkee University) : 
M.L. Mittal (worked with Prof. О.Р Varshney, Roorkee 209 
С. Prasad (worked with Prof. O.P Varshney, Roorkee Час ern 
Neeraj Bharadwaj (worked with Prof. O.P. Varshney, Roorkee 2 


Uaday Singh (.1.Т,, Roorkee) 
Vishnu Narayan Mishra (worked with M. 
Smita Sonker (I.I.T., Roorkee) 

Shailesh Kumar Srivastava (LLT, Roorkee) 
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3.6. KURUKSHETRA 3. 
Niranjan Singh, who worked with S.M.Mazhar in the Department of 
Mathematics, Aligarh Muslim University, left Aligarh to join the Department of ў U 
Mathematics, Kurukshetra University, Kurukshetra. He was a grand pupil U 
ofB.N.Prasad. He along with his research students worked there in the Theory of re 
Summability and its Applications to sequences and series. Niranjan Singh retired Т! 
from KurukshetraUniversity as a Professor of Mathematics. Some of his pupils, 1 
who worked with him are: 2. 
1. Babu Ram 2. Rajiv Sinha 3. Pushpa Suri 3. 
4. Vinod Bhardwaj 5. K.M.Sharma. 4 
3.7. КОНТАК 5 


Babu Ram, who worked with Niranjan Singh in the Department of 
Mathematics, KurukshetraUniversity ,Kurukshetra, migrated to Rohtak to join Si 
the Department of Mathematics, M.D.University, Rohtak. He worked in re 


Summability Theory and its Applications and supervised the following research m 
students: 6. 
1. Renu Chug 2. Suresh Kumari 3. Satvinder Singh Bhatia T. 
4. Navneet Hooda 5. Kulwinder Kaur 6. Jitendra Deep Kaur. 8. 

Babu Ram retired from M.D.University, Rohtok, as Professor of 9. 
Mathematics, and in the meantime his first research pupil, Renu Chug became 10 
Professor of Mathematics in the M.D. University. 11 
3.8. BARODA(Varodara) е 


It was 1958, when U.N.Singh joined The M.S.University of Baroda as 1 
Professor of Mathematics. B.S. Yadav, who was a student of U.N.Singh at Aligarh 
was so impressed that he came to Baroda to work with him. He got hisPh.D. degree E 
on a thesis entitled “Operator Theory" possibly the first in India. Some of his l 
works on Fourier series are referred to in a book: R.E.Edwards;Fourier Series(A 
Modern Introduction), Volume 1,SpringerVerlag, 1979. 8, 


In 1969, U.N.Singh left Baroda to join the Department of Mathematics of 


Delhi University as Professor of Pure Mathematics. After a few years, B.S. Yadav 5 
also left The M.S.University of Baroda to join the University of Delhi from where I 
he retired. л 

| VM. Shah, who retired as a Professor of Mathematics from The M.S. Ха 
University of Baroda,also worked with U.N.Singh on Lacunary Fourier series: al 
There ce been some other workers also at Baroda in Summability Theory and H 
Applications, in particular, J.R.Patadia, a pupil of VM.Shah. Thefollowing are 4 * 
few names who worked at this centre of research: 
1. U.N. Singh 2. B.S. Yadav 3. VM. Si E 

. BS. . VM. Shah 1 

5 O.PGoyal 4. C.M. Pate 8. 


сс: 4.Б. ац oection, Hare Agtewabi Initiate, N.V Patel. 
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3.9. SAGAR 


f sages University(now Dr. H.S.Gaur (Central) University) is the oldest 

f ү University of М.Р The Department of Mathematics and Statistics of Saugar 

1 University was headed by M.L.Mishra, themost senior pupil of B.N.Prasad. Before 

Ї retiring from this university, he educated a number of scholars in Summability 

d Theory and Applications to Fourier series and allied series including the following: 

Basudeo Singh (later migrated to Vikram University, Ujjain ) 

PL. Sharma (Saugar University) 

О.Р Varshney (later migrated toRoorkee University) 

Sushil Kumar Shrivastava (Saugar University) 

B.N. Sahney (migrated from Sagar to Canada). 

f After the retirement of M.L.Mishra, it was the turn of his pupilsBasudeo 
Singh and PL.Sharma and their generation to educate scholars in this field. It is 
really difficult to name all the workers at this centre of research but we could 

h manage to report the following: 
6. B.Deokinandan (worked with Basudeo Singh and migrated to NCERT, Delhi) 
7. S.B. Dubey (worked with Basudeo Singh retired from a college at Gwalior) 
8. Mradula Shah (worked with Basudeo Singh and migrated to Jodhpur) 

yf 9. B. P Katariya (worked with Basudeo Singh,Sagar) 

e 10. PD. Kathal (worked with PL.Sharma, retired as Govt. College Principal) 
11. R.K. Jain (worked with PL.Sharma, Sagar University) 
12. O.P Rai (worked with PL. Sharma, retired from M.A.C.T., Bhopal) 

8 13. S.C. Kori (worked with PL. Sharma, retired from Madhav Science College, 

h Ujjain) 

e | 14. Ashok Ganguli (worked with R.K.Jain, retired from G.S.TL., Indore) 

iS 15. R.S.Varma (worked with PL.Sharma, retired from a college at Bilaspur, 

А Chhattisgarh). 


3.10. JABALPUR | 
After Allahabad, perhaps, Jabalpur had been one of the greatest centres of 


Summability Theory and Applications in the North and Central part of шин 
T:Pati, who laid down the foundation of the Department of Mathematics in the 
newly started Jabalpur University (later on renamed asR.D.University), Jabalpur, 
5 came from Allahabad University(Summability School). He stayed here for about 
5. eleven years and established а good centre of Summability Theory I 
d He educated a good number of scholars in this field including the fo g 

a | Most of them migrated to different places in India and abroad: 

Ll 711 Ahmad (A.M.U,, Aligarh) 

2. НР Dikshit (R.D.University, Jabalpur) 

Ч, G. Das (Utkal University, Bhubaneswar) 
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20. 
21. 


Ram Narayan Mohapatra (migrated to U.S.A.) 
B.D. Malviya (migrated to Canada) 
Р Panwar 


` VP Srivastava (Delhi University) 


Prem Chandra (Vikram University, Ujjain) 

D. Rath (Berhampur University, Odisha) 

G. C. N. Kulshrestha (Govt. College Principal) 
Satnam Narang (Tiwari)(Govt. College Principal) 
E.C. Daniel (a college in Jabalpur) 

Indu Karmohapatra (Sukla)(Sambalpur University) 
M.L. Chandak (Agriculture University, Jabalpur) 
Maya Mudiraj (Jabalpur University, Jabalpur) 

S. Varma (a college in Nagpur) 

S.V. Lokre (Nagpur) 

K.N. Tripathi (a college in Mirzapur, U.P) 

VS. Sai (Agriculture College, Rewa, М.Р) 

P. Ray (Migrated from Jabalpur) 

J. Banerji. 

After the departure of T. Pati from Jabalpur back to Allahabad University, 


it was H.PDikshit who took over the charge of the Department of Mathematics. 
He also educated many workers in the Summability Theory and Applications and 
in some allied topics. Some of his pupils are: 


22. 
28. 
24. 
25. 


Arun Kumar (R.D. University, Jabalpur) 
B.L.Mishra (R.D. University, Jabalpur) 
Namita Srivastava 

Anil Dixit. 


Arun Kumar, who worked with H.P Dikshit, also supervised a number of 


scholars in the Summability Theory and Applications andin some allied topics. 
Some of the scholars who worked with him are: 


26. D.P Sahu (Raigarh, Chhattisgarh) 
27. L.K. Govil (Govt. College, М.Р) 
28. S.N. Singh 

29. — S.N. Dubey. 

3.11. UJJAIN 


Department of Mathematics and Statistics 
University), Sagar, 
he worked with Laurent Schwartz and obtaine 
India, he could not stay inSagar for a lon 


BasudeoSingh who worked in Summability Theory and Applications in the 


| » Saugar University (Dr. H.S. Саш 
eft the Saugar University to Paris for further studies whe!? 
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| the meanwhile, G.S.Pandey, 


Ujjain, as Professor of Mathematics to establis 


h a new Department(School of 
Studies in Mathematics and Statistics) and wor 


i ked there till his retirement. In 
| a pupil of D.PGupta, and Prem Chandra, who worked 
with T.Pati at Jabalpur, also joined the School of Studies in Mathematics and 


Statistics and educated several scholars in Summability Theory and Applications. 

Some of the names of the workers who were associated with this research centre 

are: 

1. Basudeo Singh (Vikram University, Ujjain) 

2. G.S. Pandey (Ph.D. & D.Sc. both from Vikram University, Ujjain) 

3. N.V. Deshpandey (worked with G.S. Pandey and retired from D.A. 
Vishwavidyalaya, Indore) 

4, Prem Chandra (Ph.D. from Jabalpur and D.Sc. from Vikram University, Ujjain) 

5. Ashutosh Pathak (worked with Basudeo Singh, Ujjain) 

6. Anurag Jaiswal (worked with Ashutosh Pathak, Neemuch) 

7, Archana Beohar (worked with Ashutosh Pathak, Sonkachha, M.P) 

8. B.P Khoti (worked with D.PGupta, Science College, Ujjain) 

9. В.К. Beohar (worked with G.S.Pandey, Vikram University, Ujjain) 

10. B.L. Jadia (worked with B.K. Beohar, retired from anadministrative post of 
M.P Govt.) 

11, В.Р Raghuwanshi (worked with Prem Chandra, Govt. Polytechnic, Ujjain) 

12. H.C. Jain (worked with Prem Chandra, Shahdol, M.P) 

18. H.C. Sharma (worked with Prem Chandra , retired from R.B.S. College, Agra) 

14. J.K. Gupta (worked with Ashutosh Pathak, Barwani, М.Р) 

15. Jyoti Chaurasia (worked with Ashutosh Pathak, Itarsi, M.P). 

16. L.N.Jatwa (worked with Basudeo Singh and retired from Nani-Daman) 

17. Mahesh Dalke (worked with Ashutosh Pathak, Mhow, М.Р) 

18. Manisha Chaure (worked with Ashutosh Pathak, Khandwa, M.P) = 

19. M.M. Sharma (worked with B.K.Beohar,Govt. Madhav ScienceCollege, Ujjain) 

20. Neha Beohar (worked with Ashutosh Pathak, Ujjain) : 

21. Rajkumar (worked with Basudeo Singh and E us ва 

22. R.C. i (worked with D.PGupta. Govt.College, Dewas 

23. a with Prem Chandra , retired from эы i 

24. Sunita Gupta (worked with G.S.Pandey, Govt. College, E Ms 

25. Saroj Joshi (orked with Prem Chandra ,Govt. Sci. College, Ratlam, М.Р 

26. Suresh Jain (worked with Ashutosh iren зоны М.Р) 

27. S.K. worked with B.K.Beohar, Ujjain 

28. 5р ХЭЛЭ Цай with G.S.Pandey, Govt. Sci. College, Rewa) 


i lege, Rewa) 
29. T.PSingh (worked with Prem Chandra, Govt. ws Uni) 
30. Tikam Singh (worked wikb, D-EGupte, Елве: puerum Initiative 


94 


81. ао Gupta (worked with Ashutosh Pathak, Ujjain) | 
82. VarshaKaranjgaokar(worked with Prem Chandra,Govt. Sci.Gollege, Raipur) at 


33. VK.Gupta(worked with Ashutosh Pathak, Ujjain) i th 
34. VS. Yadav (worked with Prem Chandra, Agra) u 
35. Yusuf JalaliBafti (worked with Ashutosh Pathak, IRAN). al 
3.12. CUTTACK U 


The oldest centre of Summability Theory and Applications in Odisha ce 
(Orissa) was the Department of Mathematics, Revenshaw College, Cuttack, which 1 
was established by RamanathMohanty, the first and the most outstanding 5. 
Mathematician of Odisha and the first President of the Orissa Math. Society, after 8) 
coming back from London. Earlier, after completing Honours course in Mathematics 3 
at the graduate level, he went to Patna to complete Post Graduation in Mathematics 
from Patna University and served there as lecturer for a period of one year. Later th 
he joined the Department of Mathematics, Revenshaw College, Cuttack. In 1945, w. 
he got the opportunity of proceeding to U.K. to work with J.E. Littlewood, and of 
came in close contact with Mathematicians of world rank like G.H. Hardy, J.E. St 

Littlewood, L.S. Bosanquet and others of the British School and within a couple of U 


years he submitted his Ph.D. thesis to London University, under the supervision Pi 
of L.S. Bosanquet, one of foremost analysts of the Hardy School. R. Mohanty spent 
the best period of his active academic life as teacher and researcher in Revenshaw ki 
College, Cuttack. His residence was a haven for scholars from different age groups jo 
includingBrajabehariPatanaik , T. Pati, G. Das and DiptyRath. of 
Mohanty's researches in London and thereafter weremostly concerned with Se 
Fourier Series. Several of his results led to many later investigations. For example, | le 
his criterion for the absolute convergence of Fourier series at a point is the first Fe 


theorem of its kind. His investi 
Lebesgue Fourier series at a poi 
Pati, N. Matsuyama (J apan), B. 
with several workers like S.Iz 


gations into the absolute Rieszsummability of the 
nt paved the way to further studies like those of T 12 
D. Malviya, Prem Chandra, etc.: He also collaborated Р 


umi(Japan), B.Mishra, M. Nanda, J. Manmadha Као, C 
PC. Rath, B.K. Ray etc. Some of those i 


, who wereeducated by him forresearch are: Е 

1. B.K.Mishra 2. B. Mishra 3. B.Patra 4. G. Samal of 

5. J.N. Patnaik 6. М.К. Nayak 7. M. Nanda 8. N. Tripathi G 

9. S.Mohapatra 10. S.P Misra . U 
Further details may be found in Jour Ori | 12- s 

dos . Orissa MENS Soc., Volumes 

3.13. BERHAMPUR | t 
Perhaps, Department of Mathematics of Berhampur University was th? J 


CE EE_@_—_ S S 


21 
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second generation of B.N. Prasad’ 


8 School of Summabilit 
educated a number of persons in Su lüty at Allahabad. He 


mmability Theory and retired from 
| there. 1 
the meanwhile, D. Rath, a pupil of T. Pati, also joined the Department of мана е 


under the dessen of Nand Kishore. Rath, who was heading the Department 
and carrying out researches in Summability Theory, also retired from the 
University. The following are a few names of those who were associated with this 
centre of research: 
1. Nand Kishore 2. D. Rath 3. U.N. Samantray 4. G.C.Hota 
5. H.C.Jena 6. U.K.Mishra 7.Md. Abid 8. PPadhey 
9, Sakuntala Swain. 
3.14. SAMBALPUR 
When the University teaching Department of Mathematics was opened by 
the SambalpurUniversity atSambalpur, R.N.Mohapatra and InduKarmohapatra , 
who worked with T.Pati for their Ph.D. degrees at Jabalpur, joined the Department 
of Mathematics and began teaching in Mathematics and carried research in 
Summability Theory. However, Mohapatra could not continue inSambalpur 
University for more time and left Sambalpurin 1970 for abroad. Presently, he is a 
Professor of Mathematics in the University of Central Florida, Orlando, U.S.A. 
After some time, Sambalpur University was shifted to BURLA, a place 10 
kms. away from Sambalpur and GokulanandaDas, who was in Delhi University, 
joined theSambalpurUniversityas Professor and Head,Department 
ofMathematicthes.Earlier, G.Das started his teaching career from G.M.College, 
Sambalpur in 1960 and started meeting with T Pati(at Jabalpur). He proceeded on 
leave from G.M.College ,Sambalpur to Jabalpur on C.S.LR. Senior Research 
Fellowship for the period of two years to work with T. Pati for research work. 
He completed Ph.D. work in 1965 and came back to Sambalpur. P 
awarded Ph.D. degree in 1966 by UtkalUniversity, where he was РА ts 
Ph.D. degree. From September, 1967 to July, 1969, G.Das went to ee is 
Commonwealth Scholar to work with L. S. Bosanquet at the University маны, 
London, for Post-Doctoral Research work. He maintained his cae 3 Өө: ; 
of Jabalpur who was his Ph.D. supervisor. After coming back m n хө s БН 
G.M.College, Sambalpur but he could not stay there for pis ao Е а xm 
University as reader in Mathematics. In the meanwhile, he w 
degree by the University of Jabalpur, ш 1905 hed the Department of Mathematics 
In the Sambalpur University, he establis Functional Analysis and Sequence 
to carry out the researches in the fields of 7 ility Analysis and Fourier series. 
Spaces (11) ErgodicTheory and (10) Бие this university also and migrated to 
owever, G. Das did not stay for more time п 
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Utkal University, Bhubaneswar as Professor and Head, Department of | 
Mathematics. As a result, only Indu Karmohapatra (now Mrs. Shukla) was left to 
continue the work in the Theory of Summability, from where she retired ag ў | 
Professor of Mathematics. Some of the names of the scholars who worked at ] 
Sambalpur (Colleges and University Department) are: 
1. M.M. Nanda(G.M.College,Sambalpur) 
2. G. Das (left forUtkal University,Bhubneswar) 
3. R.N. Mohapatra (migrated to U.S.A.) 
4. Indu Shukla (retired from this University). 
5. S. Nanda (worked with G.Das, 1971, retired from IIT, Kharagpur) 
6. S.L. Mishra (worked with G.Das, 1972) 
7. S.K. Mishra (worked with G.Das, 1981) 
8. S.K. Sahu (worked with G.Das, 1983) 
9. B.K. Patel (worked with G.Das in Ergodic Theory, 1981) 
10. P Mohapatra (worked with G.Das, 1972) 
11. RC. Mohapatra (worked with G.Das, 1975) 
12. N. Guru (worked with G.Das, 1985). 
3.15. BHUBANESWAR 


When G. Das joined the Department of Mathematics of Utkal University, 
Bhabaneswar, perhaps, there was no one who was working in the Summability 
Theory and Applications in the Department but soon after, a number of scholars 
joined this research centre which was developed by G.Das for research in (i) 
Functional Analysis and Sequence Space (ii) Summability Analysis and Fourier 
Series (iii) Fixed Point Theorems (iv) Probability Theory (v) Approximation Theory. 


However, somepersons who joined this research centre were from the different 
colleges of Utkal University, 


inresearch, collaborated with ( 
M.M. Nanda (4) Ram Naraya 
R.C.Dash (8) J.PDebata (9) Elli 


yo pe es, SD) em СШ I TL (2-1 гч 


У 


ГО) (3 А (255 (Үе) (о2) р су сл ыы сэ NO к Ср = 


and many more. 


+ 
=> 
tS 


ы. кш a remarkable and valuable contributions, I wish to mention 2 
3 E ve ae London Math. Soc.(3), 19(1969), 357-384], he developed 4 
е technique for a Tauberian Theorem which was adopted by a number 
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of workers іп their research works includin 
and very important contribution on Sum 


handle factor problems quite effectively. S 
the determination of new sublinear functio 


 С" 


ome very important contributions on 
nal that generate Banach Limits have 


been published in three papers that appeared in Jour. London Math. Soc. in 1973. 
After retirement from the Department, he took over the charge of Vice-Chancellor 
of Utkal University and continued his research work and supervising a number of 
research scholars. Later, he took over the charge of the Chairman, Governing 
Council, Institute of Mathematics and Applications, Andharua, Bhubaneswar, and 
kept alive “Summability Theory and Applications". He is still attached with the 
research centre through his collaborators and research pupils and continuing 
research work. B.K. Ray, who worked with R. Mohanty and retired , is also associated 
with this research centre . I am giving below most of the names of the scholars 
who worked with G.Das in various fields in Utkal University till 2012: 


1. J.PDebata 

2. M.N.Mishra 

3. M.Mohapatra 

. Tulika Ghosh 

. LaxmidharGiri 

. BijayMohanty 

S. Mohanty 

. Sangita Mishra 

. Minati Mohapatra 

10. Premlata Mohapatra 
11. Arati Mishra 

12. Anusaya Nath 

25. Krishna Manjari Sahu 


13. A.K.Ojha 

14. G.K.Panda 

15. K.C.Panda 

16. Sumitra Patel 

17. A.K.S.S. VK. Rao 

18. B.N.B. Ray 

19. Sudhakar Sahoo 
20.ManojaManjari Swain 
21.PratimaSarangi 
22.M.Samal 

28. E. Savas 

24. Dr. Nandita Tripathi 


For further details, one may refer to Jour. Orissa Math. Soc., Oe eae 
2005). Finally we give a Road-map of the scholars who migrated from a ua 
centre to another in India. To read the map, we indicate the meaning of notations 


used in the map. 


D XE 
А-—›В °С A>B 


indicate respectively migrations of pupils and grand pupils 


of B.N.Brasad from the research centre A to B. 
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(ii) I+II indicates pupils and grand pupils of B.N.Prasad. 
ROAD-MAP 


Baroda 


A = Allahabad 


Bh = Врати, ёс 
C = Cuttack | ын” 
та Tabula. 

k= K wuakeh a 


Sam = Sambakbur, 
S = Sagor 

| 1 | У = VARANAST 
222565 Pref: BN. радд 
т = тч вм et Prot. B.N. (28777) 
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ABSTRACT 

In this paper, a mathematical model is formulated in order to study the 
effects of infected immigrants and partial failure of controlling strategies for HIV 
positives in heterosexual population of uniform promiscuous behavior. The model 
incorporates a newly launched controlling strategy, named as Detection, 
Treatmentation and Rehabilitation Program (DTRP) and its partial failureness 
due poor socio-economical conditions, illiteracy and myths of targeted population. 
We also introduce the immigration of infective males and infective females as 
functions, which are proportional to their corresponding density of sub-population. 
We show that there exists a positively invariant region where the model E 
epidemiologically and mathematically well-posed. The model does not шоо 
infection-free state in general so we have investigated the criteria for T io 
endemic Steady State along with positively invariant region of P ad Be vm 
compute reproductive rate and doubling time of infection using in 5 19 bed 
state which is obtained in a specific case when inflow a шэг t set 
the population are restricted and the whole population is йг Ногоо 
immigrants. Finally, we determine local and global rus and make conclusion. 
of the system using stability theory and computer та 05 
а Mathematica ае омар Population, 

ey Words and Phrases : r : 


p ^5 1 Ш roduction rates. zt 
rede Human Immunodeficiency Virus im bean 6 

. Intro б ae 
of the maj ic health problems in the World. HIV has | Ў г 
major publio рн ce in 1981, making it one of the mos 


S : 5 еп 
million people globally since its emere 
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destructive epidemics in recorded history. The infection continues to inflict 4 
significant morbidity, mortality and social-economic amc public health burden, 
Since vertical transmission of HIV-virus can occur during р ерпапсу, delivery o; 
breast feeding. It is influenced by many factors, including maternal viral 1044 and 
type of delivery. It is documentary evidenced that, about 20010 of the children 
infected with HIV virus develop AIDS symptoms in the first year of their lives апі 
more of them die by the age of 4 years, more over up to 80 years of infected children 
develop HIV/AIDS symptoms at school age to adolescence age, CDCP[6]. About 
40% of all transmission cases of HIV virus result from mother to child 
transmission. So, HIV infection is a serious threat to future development of 
developing countries as India etc., where the level of literacy is low, the poverty 
level is very high and quality of health services is generally very poor. 
Various preventative and therapeutic measures have been embarked upon, 
alarming at combating one of the greatest pandemics of modern times. As the 
current treatment therapy has been proven beyond reasonable doubt to reduce 
transmission, it is imperative to detect those who are infected and put them on 
treatment when they eligible, so they will be recovered from their infectiousness, 
(if it would be possible). Therefore, it is desirable to encourage voluntary testing 
that will increase case detection and also promote the researches toward desirable 
mediations or therapies. So it would be helpful in reducing the number of secondary 
infections in males and females by providing treatment to primary infective 


generations. Currently antiretroviral drugs (ARV) are known to be effective in 
lowering viral loads. 


et.al.,[16]). 


In recent years, numerous mat 


hematical ively 
developed, to explore theoretical as M. cen exten® 
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e | that also incorporate vertical transmission. Brauer et.al.,[3], considered models 


en. for disease ош Mentical transmission with non-linear population dynamics and 
17 finite carrying capacity, and analyzed the stability of equilibriums in the special 
and case in which the overall birth rate does not depend on infective population size 
ren Li M. Y. et.al.,[12], proposed a model for an infectious disease that spreads in ne 
and host population through both horizontal and vertical transmission. 
ren A little attention has been paid to study the role of vertical transmission 
out in HIV/AIDS models. Agarwala et.al.[1] developed a density dependent HIV 
ild transmission model for Canadian population by taking into account the vertical 
| of transmission and by using simple mass action type interaction. Also Tripathi 
rty et.al.,[17], examined a similar model to Agarwala et.al.,[1] with vertical 
transmission but did not consider infected immigrants coming into the community. 
ОП, The major difference between past frame work done and our study given in this 
the chapter is that they do not consider HIV immigrants and also do not visualize & 
on targeted population. And we considered not only HIV positive immigrants but 
OSs, also their inflow, in infection circulation, is proportional to the density of their 
ing corresponding sub population, which is become an important aspect of HIV/AIDS 
ble in developing countries. 


In this chapter, we formulated a deterministic model in the light of above 
considerations. The model also incorporating visualization and evaluation of a 
newly launched controlling strategy, named as Detection, Treatmentation and 
Rehabilitation Program (DTRP) and its partial failureness due poor socio- 
economical conditions, illiteracy and myths of targeted population. Thus, this 
study is proposed to assess the effect of infected immigrants and partial failure of 
controlling strategies for HIV positives in the distribution dynamics of HIV /AIDS 


ase 

inn 8t heterosexual setting. | 

ous 2. The Mathematical Model. Let at time /, PU) is die targeted 

ave heterosexual population for our study. We divided this population into five sub- 

thi classes, which are denoted as S,(£),1, (t), S (£), I, (!), are known as susceptible males, 
infective males, susceptible females and infective females respectively The E 

2! sub-class denoted as D(t) and visualized, for those males and females in a- 

с 1 iti tment, when they 

i ositive and then got trea : 
< oM EET es cations or both, which would able 


uce evaluate the best strategies that would able to suppressed HIV/AIDS in any given 


Were eligible by providing such therapies/medi 
to recover their HIV positiveness and finally 
according to a newly launched controlling strategy, 


del, 
els Treatmentation and Rehabilitation Program (DTRP). In the model, g, 


they availed rehabilitation policy 
named as Detection, 


D(t) is the 
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proportion of those individuals, who were successfully detected and got treatment, 


when they eligible. But they could not received effective rehabilitation, which may 
be happened due to the scarcity of proper education, employment най also of many 
mythological issues appeared in their life as happened in developing countries, | 
like India. So, they are supposed again as susceptible to infection circulation. Ang | 
is the small proportion of those individuals, who were detected, got treatment and 
successfully availed rehabilitation. So they are emigrated from infection circulation, | 
We assumed that susceptible male and female become HIV positive following 
contacts with the HIV infected female and male are at rate k. At instant of time 1, 

| a constant inflow of susceptible males and females in targeted population are at 

| 

1 


rates 0, and 0, respectively. We also introduced the immigration of infective males 
and infective females as functions, which are proportional to their corresponding 
density of sub-population and defined as f,(J,)=5,2/, and /,(1,)-56,1/,. Where л 
is proportionality constant. In the model we ignored birth process in susceptible 
individuals because it is a sexually transmitted infection. Since, the individuals of 


D(t) population are exposed, so they supposed to be incapable for producing their 


children. In the model, pb’, gb’, гапа p’ are rates of same sense that we have 
defined in previous modeling. These above assumption make our model different 
to other model have studied on the field. 


The dynamics of infection' is assumed to be governed by the following system 
of non-linear ordinary differential equations under the above assumptions: 
S, - 0, + pb' BI, € I,)+hg,D-kS,I, - r$, 


| 
| 1, = бул, + qb' II, - L) € kS I, (го), 


| S6 «p-BU, +1,)+(1—h)g,D - kS,I, - rS, 

I, = бул, + qU'Q - DI, * L)* kS,I, =(r'+0)I, , se» 
D - c1, +1,)-(r+g,+4,)D 
with initial conditions 


S,(0) = Sip > 0,1,(0) = Т 20,D(0) = D, >0 
S,(0) = $, >0,1,(0) = L, > 0 


T 2) 
p+q=1,P@)=S,@)+1,(t)+S,(¢)+ 1„(@)+ DE) a p 
The total population changed is 


Р-(0,-0,)-(т--57-г9(1, *L)-r(S, 48-08 *g)D 
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Here c 1з the rate at which infective males and females are acknowledge into D(t) 
population due to preventive strategies or/and treatment therapies. 


§,,/,8,.-B),.—-2),6, are the probabilities of being male and female respectively, 
so that $, +ô, =1. The flow diagram of model (2.1) is givenin Fig. (1). 


3. Positive and Invariant Region. Let r, = min{r,(r +g, +g,)} then the 


equation (2.3) becomes as 


P (0, 462) & (€ b), + 1) - P - rn), +1,)- (81) 


By using differential inequality, above equation become as 
Р«(0,-0,)4(л--59(1,41,)-1Р» 


(л+Ь')(Т, -1,)4(0, 0) 


n 


which implies that P- 


= КЕ m+qb’—(r’ +e 
"emm 
S? = S? +S? =Initial susceptible males and females at time /—0. 
From (2.1), we get 
D=c(I,+1,)/r+&, +8. 


Therefore from (8.1), we derive 


LI b'-( c) 
P+r P = (Ө, +0,) + (7+0) - (r^ Зиг +52). ... (8.2) 


Hence complete solution of (3.1) is 


: 1: qb'- (r* c) l-e" 
PO < Poet + 6, «02 [ine b0- wwe 7 е). 


1 


... (8.8) 


Неге Р(0) is the total number of population at initial level of study. 
Thus it is suffices to consider the solution in the region 


о-1(Р,(, +1,),D:0< PS P08 (+h) (2 +1,),0<D<D)} 


under the condition (k++ 9) > (т'+ с). 


i iti invariant regio 
Henceforth, О 1$ positively in 
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where all solutions remain in Q, 
the dynamics of model system (2.1) is cons 


mathematically well posed in the region Q. 
4. Steady State Analysis. The endemic steady state of the model can be 


derived from the following set of equations 
Ө, + po'BU, + 1,)+ hg, D – К;1, -r$,-0 
5,11, + qb' II, + I5) * kS,T, —(т'+с)1, = 0 
Ө, + pb'(0 - B)U, + 1,))+(1- h)g,D-kS,I, - rS, = 0 
5,1, + qb (1 — D, + LZ)  kS,I, —(r'+c)I, -0 
e, + 1,)-(r+ g; + 8,)D =0 


Endemic Steady State E'(S;,I;,S;,15, D"): 


for all time. So it is sufficient to consider that 
idered to be epidemiologically and 


wees (4.1) 


. | 0, * (a, £n - e)I; +а;1; 
б = — conc |, 


S; »04 I; 20,75 » 0 and 6,x» e, 


Where a, ТЕЕ 
rtg tg 


1 2 


ip 


F- Ё +41, +a, +(8,0-e)I; | 
Diam ——————— 


S; >0,iff I; 20,7; >0 and 8,1» e. 


Where % = иас 


r+8, +8» 
г.) “1 +I) 
(rtg, +g) p 


D' >0, iff I 20,7; » 0. 
Now, we investigate criteria for positive steady state Е" to exist. We again 


use the set of equations (4.1). Firstly, 


а isocline in J; and I; from set of equation 
(4.1), we obtain RN 


a[a, —(1—8,)л]/; +[а, -(1-8)n]I = (6, +6,), (42) 


| which implies that 


I; > I, when I; 5 0and I; > 1, when I; 40. 
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| í =e А ld 
be 


where Ө = шах(0,,0,):0, | нас 
r+8& t8, 


dí __|a,-(1-8 
We can also obtain —z- Gz M. 0, 
1) аг, з 


а,-0-6,)/л 
dl; : 
am Ot a, > (1-6,)n and a, > (1- 8,)n. 
2 


Hence, we find that // is decreasing function of /;. Again, second Isocline in 
I; and Г; from set of equation (4.1), we get 


b eb, +, (I; 5) + gl} +51; =0. (4.3) 


Here 


НЫ а, +ӧл-е е О ш а+а,+ёл-е 
r r r 


80 
dete etd 
r 


| which implies that 
(nn)? +b, I) =0 when J; =0, 


such that = I;j[b1; +b, |=0, 
80 that I; >Í » 0 when J; 2007 
in | Here] =(b,/b,)>0, iff b, >0 and & 20, 
on Where 
b, оаа e and б, > Oif (a, +б,л)> 6 
r 
and 


Ay +67 EX when Ij = 0, 
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such that I;[b,I; +b,]=0 so that I, > I, when I; > 0. 
Here 1, =(b,/b,)>0 iff b, » 0 and 5, 20, 


where b,»0 if (n+qb’+h0,/r)>e and b, >0 if a,» 0 and r» 0. 
We can also obtain 

dI _ (n+qb' 4-0, Ir- e)+ (2a, /r)I, + bli MO 

аг, 25,1, + n+ qb’ +k0,/r—e 


where л+ Б + А, /г » e,6,1» e and n qb' & k0, /r» e,6,n» e. 


АЕ АЛА. 
Гал o ~ = 
= Ss е. <=. 


Hence, we find that J; is decreasing function of J; . у 
From above we see that two isoclines (4.2) and (4.3) will be interested 
provided that Г, > Î and J, > 5: Hence with the above considerations, equation 
‚ (4.2) represents Г; as decreasing from i, and equation (4.3) I, decreasing from T 
b Therefore two isoclines must intersect provided that intersection value and 
then positive steady state р" exists in I; -I, plane, see Fig. (2). Т 


5. Computation of Basic Reproductive Rate. In order to achieve | $ 
infection free steady state, it is expected that the whole population is to be occupied f 
by susceptible males and susceptible females and other sub-populations go t0 | 
extinction after some time. So, we assume that at initial stage, when л=0, the 


population P « P*-[0,/r 6, /r] - 0/r, here 0- max{0,,0,}, which is maximum 


of the recruitment of susceptible males and females. 
From model (2.1) we get 


(5.1) 


1-1 ехр[{т + qb k/r — (r4- c))1]. 23 8 


If [n+ qu'a 9 (p peg th 5 А 
5 en the number of infective die out exponentially with Ї 


time, then there is по infection at infinite time, which may implie that R, «1 | 


T qb' РӨ | 
H - + — — = Basi 
Suo in (r*c) (rc) 3 r(r'+c) Basic Reproductive rate of model (2.1.1); 1 
so that (5.1) can be rewritten ав J = [° exp((R, ~1)/t) (56.2 


Here t=1/(r'+c)= Mean time during which males and females being infectioU* 


and Г = max [ I? GRA. бишка Kangri Collection, Haridwar. An eGangotri Initiative 


p. 28 


;2 | 


ous 


Let t= t4 and J= 27° then (5.2) AE 


| 91° = Л exp((R, - 1)/1) 4, . 


which implies that 
ЭР тїа2/(В,- 1) =Doubling time of infection of model (2.1.1). 

6. Local and Global Stability Analysis. To determine the local stability 
of endemic steady state E, the following variation matrix М. is computed around 


the steady state E* : 


-(kI, +r) pb 0 pbp- kS; hg, 
kl,-r 5,2+qbl 0. 91+ 87 0 
М,. = 0 pb(1-B)-kS, (Е +r)  pb(1-B) (1- А) 
0 4801-1488) А-г  8,n+qb(1- 0 o 
0 с 0 с -(r- g, t gj) 
The characteristics equation of M, is 
НОРГАГЛАГЛТАГ ТАГ ЛЭГ ЛЫА?! +» (2) 


Here, the values of A; (for every i=1, 2, 3, 4, 5) are given in Appendix. 
Thus, by Routh-Hurwitz criteria, E*is locally asymptotically stable as it can be 
seen for 

A, > 0,(Vi =1,2,...,5); А,А,Аз > A,’ + A,’ A, and 

(A,A, = A;)(A, A, Ag E Аз = АСА 2 A;(A,A, p AY + А 45 


along with [7+ qb 20/7] > (7+0). 


п the conditions for global stability of endemic steady state of 


Now we have give 
model (2.1) in the form of following 


Theorem. If the endemic steady state E* exists then it is globally asymptotically 


stable provided (r+ g,+82)>¢ in region Q. 

По) as v, (D, v, (D, vs (D), v.) and 
ТА tv and D =D" +0;(0) in set 
pulations the corresponding non 


Proof. Take small perturbation in E (Si, 
S, = 5; +000, = 


v(t) respectively and putting 
After some mant 


of equation of our model (2.1). 
linear system is given by 


= = (5( + о) то «(рэрВ-88,)г + Pie; 


RUE Ua + U,) + г + pb -Pv +{pb’(1-B)- kS, v, x (1— h)gyUs; 
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E = R(T, n Up), + {qb'(1 -1)+ kS, |, +{б,л+ qb'( - 1) - (r' + c)? U4, 

du, : 

qi estem -(т+ в +8»); ЖО 
Let us take the following positive definite function of system (2.1) is 

У-( + Ba? + Bv? + By; + B; |/2. .. (1.2) 


Here В, аге arbitrary constant (i=1,2,3,4,5). 
It time derivative is given by 
dV dv dv dv, dv, 
pt By, + By, a + Biv, — dt t Buy — 2 ЧҮ 25: (018) 
We have positively invariant region Q defined ав 
(S;5,S,1,D) 0<S,,<S,<P, 0<I, <I, <Р, 
Q= 0«S, «S, <Р, 0<I,<I,<P}. 
0<D<D 


Here S,,,,1,,S,,,,J, are positive constants and 
P-[(n«5)(1, +1,)+(0, *8;)]/ › 
where 7, = min{r,(r tg 8:) 


© sut (S, *S, [1e (n qp'- (r’+c)}/k]; and D- с(ї, +1,)/ 


Here (S, +5,) is the number of susceptible males and females at initial level of 


(r+g,+8,). 


study under the condition (Е+л+ qb’) > (r' + c) 


А ар, 
Using the val S05 
Б the values of di ( 1,2,3,4,5) from (1) and using positively invariant region 


О given above in (1.3). we get 

dv T 

— --([RI. 2 [fr = , T 

2; ( ЙГ ru; 10253) (бүл +4Ь I) Biv? - (1, + r) By? -((r'*c)- (8r - qb' I} 


2 2 р 
бийг rte s) Ви (АРВ + pvo, рыр sr] +hg, wu, + 
1*4 1 7175 


(88: + а") в, (qb’(1-1) + ЕЗ; )В, сл, + pb' (1— 8) 


+12РВ, -(рь (1- 
Ч (p (-8)-3 5,)В зоо, +@- -h)g B vv; + cB 4U; + CB,U,u; woo (le?) 
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ee 


Applying inequality ab < (a? + b) 


dv 
рео 2 2 2 
di [ uUi + Bou; + Во» + Виш, + Bsy; | 


/2 in (1.4) and after some manipulations, we derive 


n (1.5) 
Неге 
и =| (RI, *r)- (2pb'8- kPB, - kS; +hg,}/2], 
5) 
a т +0) (8,7+71)) B, +- {рыр + APB, +(qb/(1- 1) +hS;)B, +eB,}/2], 
Bee (аг ı +r) В, -{(2pb’(1-B)+hg, - kS; 2)B, + RPB,} /2], 
3 | Ba-|((r*9)-(&x«av(- 1))} B, - [pt/8— S; + (qb/L--ks;) B, 
*(gb (1-1) kS;) B, « (pb/(1-B)- kS;) B, + RB, «cB, )/2], 
Bs, =[(r + & + 8,)B, - (hg, « (1- h)g,B, +2еВ, }/2 |. 
From (1.5) it can be shown that 2 is negative definite in the region о under ће 
following conditions: 
(21, *2r +kS; ) > (2рьВ+ АРВ, +hg,), ... (1.6) 
[lr +e) - (8, qb’ 1))В, >= [pb B+ (RP +qb'l kS;) B, + pb’ (1-В)В, 
(qb'(1 — 2) kS7) B, +ев,), 2009 
of | +) В, > z (epa -8) «ha, -®5;)В, +kPB,}, qu 


(er +e)- (85+ qb/(1-1))} B, >= (рыр - kS; + (qb/1+ hS; )B, «(рь(1-8)-86:18, 


mm | (qb/(1-1)+hS;+hP)B, *cB,], (1.9) 


... (1.10) 
(27:17) В, >а +( +(1-h)g,B, + 2cB,} . 


Р 1) The above conditions given by (1.6) to (1.10) may be further manipulated to get 
following simplified conditions: 
From (1.7) choosing B, 
p. [Pe B+ pb/(1-B) B; » (qb/ (1 es Be Je 
“1/2 
4) A (rc) - (RP +S, +28,7 + 34b 2 
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the condition (1.7) reduces to 
(rc) > (kP +kS; 4 28,14 3qb/1)/2. 2 111) 
From (1.8) choosing as i 
[ЕРВ, |/2 
B, Pl үүл PS 
[(er, +г)-{2рЬ'(1—В)+ һе, -88:/2| 

the condition (1.8) reduces to 

[2(kI, +r)+kS; ]>(2pd(1-B)+hg,). (149) 

From (1.9) choosing 


| pb’B + (qb'1-- kS; )B, 4(рь (1-8)-28,)В, «cB, |/2 


3 > E » 
(r'*c)- (kP +kS, + 28,n+ 345 (1 - 1))/2 | 


along with рБ(1-0)» А5; ,the condition (1.9) reduces to 


; 1/5 5 
(r +0)>=(RP+kS, * 28,1. 3gb' (1-1). оо (Л 
From the condition (1.10) choosing 


m. [^g, +(1-h) ев, ]/2 
(r +8, t8; ) gor 
the condition (1.10) reduces to 
(r+g,+8,)>c. (1.14) 
Thus, endemic steady state E* i 
Q under the condition (1.6), (1.11) to 
7. Numerical Analysis. To 
(2.1) numerically, the model is inte 
following set of parameter values: 
0, = 0.840,0, = 0.893, p = 0.700, b’ 
б, = 0.500,5, 


8 globally asymptotically stable in the region 
(1.14). Hence, it proves the theorem. 

study the dynamical behavior of the model 
grated by Runge-Kutta 4th method using the 


di он = 0.8018 = 0.200, g, = 0.200, k = 0.578,r = 0.400, 
MS EM is 0-900,q = 0,300,1 = 0.200,r’ = 1.000,¢ = 0.600, g, = 0.136 
with initial conditions §,(0) = 0.5, 1,0) = 0.3, 8, (0) = 0.5,7,(0) = 0.3, D(0) = 0.6 the 


endemic steady state values are com i 
8 puted as shown i in Fi ze) 1025 
seen from the Fig. 3(a-e) that for any initial start: meme шш АШ 


(1) 5,(0)= 0.500, J, (0) = 0.300, S, (0) = 0.500, T,(0) = 0.300, D(0) = 0 600 
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(2) S,(0) = 0.745, I, (0) = 0.289,S,(0) = 0.845, I, (0) = 0.624, D(0) = 0.414, 
(3) S,(0) = 0.389, I, (0) = 0.654, S, (0) = 0.452, I, (0) = 0.498, D(0) = 0.547, 
(4) S,(0) = 0.687, J, (0) = 0.395, S, (0) = 0.712, I, (0) = 0.264, D(0) = 0.870, 
the solution curve tend to steady state р". Hence we infer that the system (2.1) 


may be globally asymptotically stable about this endemic steady state pr. The 
endemic steady statevalues of different variables are computed as follows: 


S; = 1.05014931169431,/; = 0.663578072430592, 


S; = 2.45490438471322, I, = 1.11629241558393, D' = 1.45098138082658 
We also compute numerically following variation matrix of system (2.1) around 


endemic steady state jr” 


-1.0452 0.1121 0 —0.4948 0.0400 
0.6452  —1.1019 0 0.6550 0 
Ар. = 0 -0.9704 -0.7835 0.4486 01600 


0 1.6112 0.3835 -0.9578 0 
0 0.6000 0 0.6000 -0.7360 


The characteristic equation of variation matrix A,. is given as 


25 +4.626БЛ* + 1.20822? + 5.09264? +1.5789A + 0.1442 = 0 
and the characteristic roots of variation matrix А„. are 


М = —2.892518235036159788256730120598, 

dog = -0.7199584759080117792625460275619 
40.3290294260549985883494208903371i , 

М = -0.14971697550915996051421011426455, 


А = _0.64231297971265656387726903425918. 
Routh-Hurwitz’s array system (2.1) around endemic 


We also compute 
Steady state p* as 
1.000 7.2032 1.5789 
4.6945 5.0926 0.1442 
6.1020 1.5477 0 


M.z 

E (3.9196 0.1442 0 
1.3233 0 0 
0.442 0 0 
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Hence from characteristic roots of A,- and Routh-Hurwitz’s array, we infer | ` 
that the model (2.1) may be locally asymptotically stable around endemic steady la 
state E'. | E. t 
In Fig. (4), the distribution of all sub-populations of ус: (2) with time js |r 
shown. It is found that susceptible (male and female) populations increase with is 
the time due to the recruitment of susceptibles (male and female) and of new born Ч 
susceptible offsprings. It is also seen that male susceptible population decreases g 
rapidly in comparative to female susceptible population as population is not closed. 5 
Whereas, male and female infective populations, starts with small decrement and | 
then increase rapidly with time due to initiation of Detection-program then it le 
again decrease as all infective males and females have been effectively detected as t 
HIV positive and got treatment. Its impact that the detected and treatmented a 
population increases with time and then decrease slowly as they moved from the 8 
| infection circulation due to either of their successful rehabilitation or dead from fi 
| infection. Also some of them again stand as susceptible to infection circulation с 
due to their improper rehabilitation. Finally the variations in all population are le 
being stabilized due to endemicness of virus. 
In Fig. 5(a-e), shows the variation of all sub-populations of system (2.1) with р 
time at different level of horizontal transmission rates as k=0.578 to 0.800.It is t 
found that all susceptible (male and female) populations decrease as transmission | s 
rate increase and all infective (male and female) populations along with deducted | i: 
treatmented population firstly increase as transmission rate increase and then i 
reach on its steady state position. The impact of recruitment (i.e. vertical b 
transmission) of new born offsprings of infected parents directly into their f 
corresponding classes, are explicitly shown in figure 6(a-e). It is seen in Fig. 6(2) e 
and Fig.6(b) that at all increasing values of vertical transmission rate (b’), all 
susceptible populations firstly increase for some time, since population is not t 
closed, then reach to its steady state position as earlier as transmission rate being 
increased. On the other hand in Fig. 6(c) and Fig. 6(d) as transmission rate increases 2 
all infected and deducted - treatmented populations increase and reach to its steady t 
state position as slowly as transmission rate increased with time. Thus if the 8 
birth of infected new born offsprings is controlled by medication, the overall 0 


infective population might be remain 
in reducing the density of infective 
population. 


The Fig. 7 (a-e) shows the effect of recruitment of infected immigrants t0 © 
susceptible, infective and deducted-t : : of ü 
CC-0. Gurukul Kangri cacdtmented populations at.different values 


ed under control. That 15, 16 would be helpful 
Populations which will transform into AIDS ý t 
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7 n= (0.912,0.947,0.984). It is seen that as recruitment rate of infected immigrants 
y | шоог eee ( шин are! nel) populations reached earlier to its steady 
state position .1 his Increase in rate also results that all infective and deducted- 
treatmented populations increase very quickly and then stabilized as slowly as 
is rate increased. Hence, it may predict that the recruitment of infected immigrants 
Л is not only accelerating transmission efficiency but also increased viral density, 
"n which can create more complications to implement the best controlling strategies. 
os Such satiation should be controlled because it may damage the controlling 
d. strategies more than vertical transmission. 
id The Fig. 8(a-e) predict the variation of all sub-population with time at different 
it level (0.3 to 0.7) of rate (partial failureness of D.T.R. program). It is seen that with 
AS the increase in the value of rate, the deducted and treatmented population decrease 
d and reach to stabilize position earlier. Which results internal increase and late 
le stabilize position to susceptible and infective populations. Hence, it would be 
m fruitful to plane such controlling strategy that will focus not only on infection 
n controlling but also on proper rehabilitation their education and socio-economic 
re levels, so that failure chances of plane may be reduced. 
8. Conclusion. In this paper, a deterministic mathematical model, have 
h proposed and analyzed to understand the effect of infected immigrants on 
is transmission dynamics of HIV infection in a heterosexual population, a controlling 
n | strategy named as Detection, Treatmentation and Rehabilitation Program (DTRP) 
d | is running. In this model, the growth of infection of infection is increased by sexual 
n interaction between susceptible and infectives. The infection density is increased 
al by new born offsprings of infected parents. The model also incorporated a partial 
ir failure rate of D.T.R. program to investigate the impact of socio-mythological, 
a) economical facts of targeted population, where D.T.R. program 1s running. А 
Ш In the study of model, we have derived а threshold parameter Rg and note 
ot that, if [n+ qb k0/r]>(r'+c), the infection always remains in population. The 
: model has a non-negative endemic steady state Е*(8:,1,9,1,0 *). It is found 


which exist only when? Бо? 1, is always locally 


y | i tate E* Rozem Sa улс: 
va s stable under the positively invariant region; 


16 stable and globally asymptotically 


Ш obtai 1: with certain conditions. T 
ul 23 AE e an increase in the rate of vertical oS ET uae 
: 1 ive; u e 
5 (| toincrease the populations of male infective and female mr х coat ie 
turned into AIDS individuals. Thus the vertical spread of е gr. 

: i ndom use that will be Ke 
f aey eee ME х be predicted that, by controlling the 


of infective density under control. Also, 
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recruitment of infected immigrants and by implementing an appropriate Отв. 
program according the situation of targeted population, the spread of infection 
can be reduced significantly. 
It is found that an increase in the rate of recruitment of infective males and 
females leads an increase in the population of infectives as well as transmission 
rate of infection. It is shown that when viral load rapidly increase, controlling 
program has required to increase process accordingly, which may be reduced the 
effectiveness of D.T.R. program. Consequently the steady state values of infective 
males, infective females and detected-treatmented individuals can be maintained 
at desired levels. Hence the analysis of model, speculated that the HIV prevalence 
can be reduced and/or controlled by an effective strategy, which must include 
detection, treatment and rehabilitation process according targeted population, so 
( that the progression to AIDS can Бе slowed stopped down. 
| Appendix 
A, = [a-g+k-s+u)]; 
A, = [k(a-g)-b(e-/)-ag-h(n +p)-l(q-r)+u(a-g+k-s)-s(a-g+h)]; 
Аз = [(kI-h(v-w))(q-r)-(h(n+-p)+1(q-r))(a-g+k)-w (ag +b(e-/)-k(a-g) +h(n-tp)+ 
Uq-r)+s(a-g+k))-k(ag+b(e-f))+s(ag+b(e-f)-k(a-g))-(n+-p)\(gh + (b-c)(e-f))- 
dt(e-f)-mt(q-r)]; 
авг оло +(b-c)(e-f)))(a-g+k)-((v-w)(gh + (b-c) (e-f)) +h(k 
ga u(lh(n+p)+Lq-r))la-g+k)-(kl-h(v-w))iq-r)+klag 4-0(6-0) 
p o )+(n+p)(gh+ (b-c)+(e- f)))-U(dg(e-f)-ad(e-f)+hm(q-r)) 
| ы к ОН ин бон» e 
| ks(ag--b(e-f)]; pue debe 
хэс aale Dt hm (4-р). (q=r) (mk? + Im(q-r) + d(e-f)(0- 
gag r)+d(e-f\(n+p))) +u(((kl-h(v-w))(q-r)-(n-+p)(gh + b- 
OC P)Ma-g+h)-v-w)(gh+ (b-c)(e-f))+k(hI-h(v-w)))( I Sta 
T)Yag-tb(e-f)-h(a-8))-(n--p)(e-f)(bh-a(b-c))4 ane se (np EC 
f)-(Gdg(e-f)-ad(e-f) - hm(q-r))4- urs gigh - (b-c)(e-f))) - ks(ag 5 
S\a-r)-km(q-r)+ d(e-f\(n+p)))(a-gth 


А,- 


ad C= RS, d = ghe S b =r, в =[бүл qu'1] 


214 ә) вав. (1- A) 
? T3625 cT 


n-qb' 1-1), =Ё$,, = Pf (T а 
0-0 cto. G oe p b ast Нат. АМ авандорива ра, 8, £;) 
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Fig. (3) : Phase diagram of model (2.1). 
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Fig. (5) : Variation of all sub population of model (2.1) with time at different values 
of transmission rate К as k =0.578, 0.689, 0.800. 
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| ABSTRACT 
In the present paper, we define a degeneration formula for a quadruple 


hypergeometric function K,,(.) due to Exton [3,4] as introducing a Dirichlet type 


integral involving that quadruple hypergeometric function [3,4] and show that it 
is degenerated into bilinear product of double Appell functions Ё; () (cf. 1,11], and 


ive fractional derivative representations, bilinear 


again, make its applications to der 
we present some 


generating relations and transformation formulae. Finally, | 
particular cases, interesting applications, conclusions and discussions. 

2010 Mathematics Subject Classification : 05A15, 26A33, 44А20, 11B37, 
38070, 33C99, 30635 E 
Key Words and Phrases : Weighted Dirichlet type integral formula, Exton’s 


quadruple hypergeometric function К,,(.), fractional derivatives, bilinear 


lusions and discussions. 


5 i i :mteresting applications, сопс 
| generating relations, inter Барр hlet integral introduced due 


1. Introduction. From the normalized Diric ao 
to Mathai and Haubold [6], particularly, in the four dimensional space R* c C^; 


We define following Dirichlet integral density 
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(1.1) Цр )-1,4,| k | pre 
mi ! Ї ЇЇ yh ade P RET щ-1 (12 Х)-Х,-Х,-Х Ae dx \dx,dx,dx, Pai 
where, x= NC R‘ such that 0< x, *x,*x,*x,Sl, u= (Hi Hols Hap ) 
Г(р,)Г(вь)Г (нь) (ы, )F (us) 
s B(u)-» ———3— 3 а | j-1,2 
ecs, ВШ T (u, +My +H +My +0) шиш, Vi- 123,45. 
Exton [3,4] defined the following complete quadruple hypergeometric function 
among twenty one quadruple hypergeometric functions : 
(12) Kj, (a,2,2,a,5,,5,,b,,5,;c,,C,, 0; x, y, z,t) | 
= Y а) (b, ), (5, ), (b, ), (b, ) e yt e iv 
m,n,p,q=0 (с, = (с), m! n! р! 9!” f 
where, the parameters a,b 15055550 6561 ‚со are real or complex and x,y,2,¢, are complex 
M» [T 
variables. du 
Now, for a>0,8>0,y>0,8> ORe(c)»Re(-uj),Re(u)»0, Vi=1,2,3,4 and | Pr 
R ) E" 
ea- Dh, а) 0,а,5 = (b,,b,,b,,,),¢ = (coc), a = (His 15,15, 14); A, ħa, hg and | 7 
h, EC we define a function due to weighted Dirichlet type integral formula 
(1.3) Еа (^о, hB, hyy,h,5) = 
Г(а)(о) =) (py 6» (yy te) (8) 0) 
T'(u, +0, Jr (m, *o;)T(u, *o;)r( T (u, +04) c(a- У( (шс, ) | 
1-1 1 
a eB ey 6d Е к e. о 
Jf, fef, roe! әв" -ayt cag e Мани тэл мон | f 
EE 157 1 


y z с; u XE c 
рн) 2-2 con 
до Еа? 97 1-х01- 3121-8. 


= рг 
Hy-1, jua 71 н, H4-1 
x gh- 7 4 К 
ý „(йа,а,а, bib, bb; 16,,6,,6,,0,; h, Jihaz h t)dxdydzdt 
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ЭЭГ... 


19) 


tion 


| 
] 


1 
| 125 
| provided that 05 xoc! + y^ + гү! +457 «1, 


Particularly, when 6, =0, from (1.8), we have 
(4) FHM“ (ho, рв hay, hyd) = 


Г (а) (о) 99 (py) (yy 9 (sy 9 


| P (A, *c,)T (u, * o,)T (us *o;)T (p ee r[a- X. 42) 


S f rae эв -eyt ey mnm ОГЧ 


y 92 2 95 Ї СЯ 
[ -xa - Xp -zy! -t87 | Ё -xa — X9! -zy! -to | f -xa — yp" - zy! -t6^ | 
gl yel pj tp po (a, a, a,b, ,b;,b,3¢, ,65, 05; ho y hisz, hat) dxdydzdt : 


provided that 0 < xoc! + yB” 4 zy! +15 <1. 
| 
| When, b, =0, and ù, =0, (1.3) becomes 


(LB) р) (ду, RB. hy, h45) 


T (u,  o,)T (Hy * o;)F (us + 05)F (п, +о,)Г eS 


2 


27 


-zy|-16! | 


x 


n ; f Ї (1 -xa — yp -zy -tT Ed Ї -xa -yB 


9, 


92 


оз 9 Ё 
2 um le c. 
III pA шшс == = | Se ee. =, 
в) Ї xa- yB" -zy 1521) (1-хо7"-10 -zy 


4 умт отрар (а, by b4; Ca; haz, ut) dxdydadt 2 


Provided that 0 € хог! * 58^ +2Y + 17 «1. 
For, b =0, and 6, =0, (1. 3) becomes 
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(1.6) Ее ӨӨ (ро hB, h.v, hò) = | 
: ШООК (gym aye eye 
х T (p, +0, )P (u, *o;)T (u; +оз) (и, +54) P| a- 24 + 2 


a 5 E e 5 -1 M71 7124371471 x Ч! 
[н op cert amm 


y Mme y 
суг - yp! -zy - t6"! | Р -yp!-zy'-t6* | 
t р 
ej gU REC ЛИЙ (a, 5,,b, 36,, 5; h, y, h4t)dxdydedt : 
provided that 0 € хог + yp! «zy! «15? <1. 
For, 6, =0,6, =0, and b, =0, (1.3) becomes 
(1.7) p? 093009 (f o, haB, hv. 8) 


д P(a)(a) 9) (By даа (yy 992 (oen | 


шилээ 


сү ey eb e x 2 -119-11-2-13-14- | 
Бөгж temm ste 


y үз 2 Ra 
——À áo E 
- - - - ТОО СШЛЭГГ н г 
БЭ *— 38^ — zy — 48 | тээ <=) 


t j^ 
Шоо ат 
provided that 0 < хог! + yg?! + zy +161 < 1 


Here, in our investigations, we show that when max (|/o] |.) «1 end 


1 

max | чл| 8,8) <1 the function Fens (^а haB, hay, h. 8) defined in (1.3) 
? > 2° %4 

degenerates into bilinear product of double Appell functions F (.) (see Appell and | 

ae , 


Kampé de Fériet [1] and Srivastava and M Me 
рё c anoch i 2 
applications to derive fractional derivativ: Ln двар меша 


е герге 1 ЯГ rating 
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we present some particular cases, 
ussions. 


relations and transformation formulae. Finally, 
interesting applications, conclusions, and ihe 


2. Degeneration Formula and Fractional Derivative 


Representations. In this section, first we derive that the function 
а,5,,,д 
р" (ha, hB, hay, h45) defined in (1.3) degenerates into bilinear product of double 


Appell functions F,(.). (see, Appell and Kampé de Fériet [1] and Srivastava and 


Manocha [11,p.53]). Again, we make its applications to evaluate some fractional 
derivative representations. 


Theorem A. When, max {|h a,l} <1 and шах, 


у ! <1, then, the function 
paremo (hia, h B, hy, һӧ) defined in (1.3) degenerates into bilinear product of double 
Appell functions Р, (.), that is given by 

(2.1) Ре (ho. №, B, Asy,/,5) = 

F,[u, +0, ,by + 05,5, 5,0, ho, AaB] F; [bs + 05,4 +5 43051043605 №] 


provided that all conditions given in (1.3) are satisfied. 
Proof. In the right hand side of (1.3), for the given conditions 0,0,ү,6 > 0, we set 


xo eu, yp -u,,2y uy and t§'=u, and then define the quadruple 


hypergeometric function due to (1.2) and thus change the order of integration and 


summation, we find that 
(22) pates (ha, h B, hy, №) = 


D(a) 


T (u, +6, )T (m, +с„)Г (p +95 )Г( u, +0,4) 18-31 Уи + 2) 
hy 


1 ps0: m-l 
(и) 


L2 a b, |= 27 „ (Bs), (5.), (ha а)" ( 
$ (н 00.0) 


| ! a q! 
m,n,p,q=0 (с, ae (с, his 24 i 5 


а- ) (u,+0,) F1 н (2.2) 
(u, уйлан (us yerum (и, унн (1 =p, = и, =й = и, | У И } du,du,du,du, 


provided that 0 <u, tu; + +4, 51: De 
| Now. under the given conditions, Ке(с;)> Re(-u;) Re(u.)7 0, Vi = 1,2,8, 


Re P у + make an appeal to the for mula 1.1 Inr side 
( ) ight hand 


per ukul Kangri Collection, Haridwar. An eGangotri Initiative 


4 


ЕГ г" 


128 
of (2.2), and again, in it define double Appell functions F; () [cf.1, 11], for the 
conditions, max (od pnl «1,and max (Lal) <1, we derive the result of (2.1), 


| 
Theorem B. | 
(0 If Re(e, -max{b, (u: +0,)}) > 0,8е(с,-шах(8, (u: +04)}) Оно < 1,98] «1, | 


|hyy|<1 and |h,8|<1, then the fractional derivative representation of 
po^»? (ро, hB, hy, hò) defined in (1.3), is given by 
(2.3) Е" (ро, hB, hy, А.) = 


Г(с, -5,)Г(с, ны б,)Г (c, - b,)r (c; =u, -с,) (h ue +l ( -с+1 Di*oitesm-a 
T (b,)T (u, o, )T (b,)T (us +03) " 


| 
JD SEU ER (h E (h; jars sere DRS е (@ у (he jio (1 2 КОЙ Те | 
| 


*6,,5, 


(1 apt дч +0,,0 а +hp- 2217 Л hay +h,d- 20 


25 
d 

h D ———— med 

where D, dh’ and D, Zn 

(ii) If Re(c, - max (b, (і, +0,)})>0,Re(c, - max [b,, (u, + оз)})> 0,0] 51,158 « 1 

|hgy|<1, and |h,5|<1, then, the second fractional derivative representation of 


| 
| poene (ра, h,B,hgy,h,5) defined in (1.3), is given by 
| (2.4) Ее (hc, If, hY, hò) = 


(e, -b)T (e, = = 02) (e; -5,)T (c; - 20 ) 
3 3 һ2+о2-сү+1 4305-05 40, 1540-01 
F(b)r(u, *e,)r(b)r(u, +6,) (8) (a, Jerre" prot 


Hs 40g + +04- Hy +0; +b) —c, 4+04+b;- Л M 3 2 
Di: % (h) (A, BD. 5, ‘Dy b,-c; (o. (ye пој с 
13303 *55-c; 3755 
Gh) mes бий анд 221 Ie У 48 — hh Jd 


| 
5 8 
4 4 
р = = 
where, ^, dh and р, dk 


4 
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Proof. Make ап appeal to the fractional derivative formulae for double Appell 


е у 

function F,(.) due to Srivastava and Manocha (11, р. 306], we get 
). - 

буй Г(у-о)г(ү-В) 
(2.5) Е, |0,0,В,3;ү;рх,су|= Y «-ү+1 pyara’ oft BB — 
| 1 ү;рх, oy] rore) ^ De "appe 
| 
2 222 o p^ 

7 Б "(17 px) un х ре+оу-реву],р, = ©. 

апа 


T'(y - a) (y - B) 


(6) Foo B, v; px,cy]- r(o^)r() 


VAM EDS pe s 


Та (1- oy) Эд Р px + Oy — 212 = n 
Y; dy 
Then, use (2.5) in (2.1), we derive the result of (2.3) in the theorem B(i). 
| Again, use (2.5) and (2.6) in (2.1), we obtain the result of (2.4) in the theorem B(ii). 
| | 3. Generating Relations and Transformation Formulae. In this 
' | section, making some changes into the parameters of (2.1), we derive various results 
on generating relations and transformation formulae due to following theorems: 


Bl] |ВҺТ 


L <r, such that rj +r nr, 
14 IT-1 


|, | Theorem 1. If [ho] 558 55 | <1, 


<r, such that r,*r, =, then, for all 


DA | 
Е KP 


and [hay] « |р сир 


conditions given in (1.3) and À be any complex number, following bilinear generating 


relation holds 


(3.1) уо genie (^а,љ,8,№у,№8)1" = 
-n=0 B 


(А), (uz +6; ): (н, +6; 1 h,h,poT. 


- anys Ber ena 0-7) | 
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ыг! 


: +n, c +n; hey, ha 0, —— 
Ba) ha ын, 56, HM Ha +O, bb nies л A шиг 2114) |p 


where, Е,|.| is the triple series defined due to Saran [8]. 
Ч | 


in both sides of 4 


Proof. Replace b, Бу б-п, b, by b,-n and multiply by 


(2.1), and then sum n from 0 to со for [T| «1, we get l: 
[К (3:9) хэв Fubbenboberens (р o. hB, hy, 45) T” |i 
l n=0 | 


у tt) (6), 2875 (и, +03) ), (6 s); бы) у 3), 


| р-0 ee r=0 (с, ) г! п=0 n! 


= (tp = (b, - = (1, c x: (1,8) л 


Now, in the right hand side of (3.2), use following formula due to Srivastava (10, 
p.71, eqn. (3.4)] 


(3.3) S (BO Stas ту 281 | m | ї 


: 705549” | 
| = (0), (A +n) 1( Т Y & (о) (A+n) TY 
TR u U u y 
2, l!m! Үп E x u'u! 8). (2 3 || «1. 


Then, make some series rearrangements and define the triple series Е, | 


introduced due to Saran [8] (Also, see, Srivastava and Manocha [11,p.68], we obtain 
the result of (3.1). 4 Ч 


In the similar manner, we have Н | 


ABT h,ôT 
Шешә ши лах |Т c апа |7|<1, then, for all 


conditions given in (1.3) and 
relation holds 


à. be any complex, following bilinear generating 


S (4), a 45 ,-n,b,-n),c,u,c | 

ВО Ye Phere 77737718 
= (2 Ey? yo (He +0), шоо al 

= ale а), 2 852, Gangotri Initiative 
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Е ТОС b nie +n; ha, т 


2 EUM CE ML 
-1 . 
freonem 3 Lt хо < „В| «вө «1, and |hs'y| < fh, 6|«c, then, бог all 


conditions given in (1.3) and } be any complex, following bilinear generating 
relation holds 


< À -n +л),с,и,с 
| (3.5) $ C. реб 7n b; ,6,*n)c p, (ha; hB, hay, h ò) T" ү, 


туа ха g(t +в), (и, +о,), | aBer)" 
(1 Т) 2, п\(с), (е), | = 


3 (C): (ius +0, +7); G ): (u, * 6,5, ), (6) (+n): h,BT 
pe [ OOOO. адолат 


p (s (—);(и, ДЫ +n), 4) (Us + 05,b,), 3(b, = А), (А+ n), : 
Е! | (с +п), Ми (-)3(-)s(-):(-)3(-)s(-): hyy,h,, То 21 


| where, por. ..] is the general triple hypergeometric function due to Srivastava [9]. 


Theorem 4. If |h,o|« =,|^,8|<=,|Т|<1, and || <= ,|h,3| « =, such that then, for 
all conditions given in (1.3) and À be any complex number, following bilinear 
generating relation holds 


(3.6) 50). кезүн (ро, hB, hy, A 8)T"^ = 


n=0 n! 


(A), (Ha +02), (Ha 54), заа 


Mo e 1127 


3| V4 =); y t0 +7),5( Ne (u, #155, ны -A) о: ha Вт Ee 
zi )::( еше ХООЛ a 
а) =) (ity tonb) ib 20) 3A р mu | 
хил Ши 


i i Srivastava [9]. 
Where, F®[...]is the general triple hypergeometric бариад dua юу 
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AY | eo 
Theorem 5. If max Tio], In; pl «1 JU sib ren ТАА « 4В,9| <1, then, for ай conditions 


given in (1.3) and д be any complex number, following bilinear generating relation | Р 


holds 18 
Е 
(3. 1) хэ. A peh (то, hf, hg, Mg §)T" | 
п=0 | 7 
| 
=, (ho) g В), 0) (uto), | ABT | 
Чахол) элр (Љо) (0), (Ha +о), (5)... 09), s * 64), eae | 
-11-8,5) (1-1) Ус) 5 m! һо(1-Т) | i 
Fo) (~):()3(c -.),:(7) : (He +52 +т,,), (А + тра +05), (84 +04 Js уй hy hô =] E 
Ө): +) C): 1-7-1 | 
A) T" | 
| Proof. In (2.1), replace b; Бу Б,-п and b; by А+ папа then multiply by ( л in | | 
both sides to it and thus sum п fro 0 to со for | «1, we may write 
< (А), a (b b -n nb.) o " 
G9 York (Iso, B, hy, hið) T" = 
ж h, 5 (4*4) (17 +С 18 (b 1), (ha)" ы л), b, =N, u, + Ga: 
=( ) 2 т(с,), 2, ni oF сү +m; hp 1 
xen -b: quse, AB hd ла | 
C2: C= b, h; 5- 1 | 


Now, for |T| «1, in the right hand side of (3.8), making an appeal to the formula | | 
due to Srivastava and Manocha [11,p.155] | 


= (А), p-n,(g,); rape (е a): X n.a, (c 1 Я 
вэ» Жы | (hn); ds (6):(6)s(.) her 2 


5 0), 16), 
(-7/:У--8-- = 
=D п] [(^;), 


ja 


Th oe NNI 
OEA: (nm) #7), s(6,)) (8, кел 
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where, Fi | isthe 1 2 
ar] two variables Kampé de Fériet function due to Srivastava and 


da [12 (8) ! 
Вап P ] and F [+] is the general triple hypergeometric function due to 
Srivastava [9]. Then using series rearrangement technique > 
(8.7). ques, we obtain the result 


hy 
TT 


Theorem 6. If max [[/o| |^; [) « 1,[T| « 1, 


< |[^48| « 1, then, for all conditions 


| given in (1.3) and д, be any complex number, following bilinear generating relation 


holds 


| = (А), а -njÀtn. сн 
(10) Уу Efe mie (ho a hy hy)" = 


| given in (1.3) and } be any complex number, following transformation formula 


a-ra" a-p $ BA у Creed... ео par. T 


та ЇЕ m=0 m! 


ha(1-T) 
=H, 501) ;(—): (He +0, +m,b,) ;(A +7, ц, +03), :(6,), : 


hp L M 
Theorem 7. If max ЇЇ 251 «1, | «1,шах л lhl} <1, then, for ай conditions 


holds 


e (A 
(8.11) yg pne (ha, hap, hy, h,d)T" 
n=0 о 


4540) (He +02), (1*9), [ёт |" 
= 1 = Т А ( n 2 2/n n 4 
eMe | 
hô 
pet A*n,bs Хаа) (eene tn) ua nolis Mun) (ro Af ; 157 1 А a 1 
<1, then, for all conditions 


Theorem 8. If ах о « Ir «mox АА) 
given in (1.3) and à be any complex number, following transformation quam 


holds 
(3.12) У O), pot Ens „А+л),сс (ha, hb, hY, hj5)T" 
n=0 n: 
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130) (to +r), (44 +50), | AI B9T | 
= (1-7) P» те | d 


ҺВТ , M | I 


groben bin) mazim) n tainn uin) | кл" abl T 


4. Special Cases and Interesting Applications. Setting À--ha, | 
| 


: 1 

hg = —hgy,b, =p, + с, =} * 65,5, = +0, and b, = t0, M (2.1), we find that 
(4.1) рбетот) со (Ао, -ha hy, - hy) E ( 
2 H, +O, H +с,,(щ +0; +} +0,)/2,(0, venari нот DUE. i 
wes р, +0, + Hp +оз,(с,)/2,(6)/2+1/2 | 


Us 03,1, C, (I 03 +, 0,)/2, (I5 +05 +H, +0, +1)/2; TOY 
Fi ! (hey) |. 
Is +0, +, *6,,(c,)/2,(c)/ 21/2; 
| Again, setting 6,=1,b,=1,b,=1, and 6,=1,c,=y,+0,+",+6, and 


€; = H, +0; +p, +0, in (2.1), we find that 


(4.9) pou so ates datos tuo, ) о (ha, 1,8, hay, h,5) = (ha + hp - hh) 


ООЛ go havran Bie | 


Hi tO, +H, +0,3; Hı +O, +H +03; 
Из +03,1; H4 +04,1; 
hy 8| +h,5,F ANC 
| 21 о шо 2 Tdi ожа ^ F 
Setting с,-28,,с,- 2b,, 1, = b, —1/2—-вс»„,ц, =b, -1/2 — 6, in (2.6), then, we find 
a b (2b, 25;)(u 5, -1/2-0, 5,-1/2-03 в,)а 
(43) Р : ^9? (ho hB, hy, hò) 


= | 


sa T (b, +1/2)Г(Ь,)г(Ь, -1/2)r (b 8 
= 4 ун EU 3172) Г(5,Г(Б, 4172)Г(5,) 1/2-8, 1/2-6, 1,10, 75, 1/2 1,504,751? | 
К Tee IPG Л, (А) "е Dt 


(о) ав) ^ (1 TED | 


provided that Ве(5,)»-1/2,Не()»-1/2 | 51 «1 «1, and |һ,8]<1- 
Now, from (4.8), we obtain ап interesting transformation fo 


: rmula 
(4.4) preter hae 


(а, 6, hay, h,8) = 


e b,+b;-1 M* R 8,1/2-5/2,1, +6 5 
( ) 3^2 b, b 41/2: "hy (1,8) 
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b ,1/2- 9/2 to, 
М” E| 4 2:7 4? 
| | 5,,6, +1/2; м |ы) о (4.4) 


provided that Ве(5,)»-1/2,Ке(Б,)»-1/2, 0) <1,|h,8| < Lare( 11-58) cn m 


ML <x 
e | Here, the Mellin-Barnes transformation operator is defined as 
| 1 pe T(25,-1«&)r(-5)1 - ng" 
i М” h,B) = 2 2 t 
(4:5) {0(6)} B) = — f. rB -1) (9(5) 4. 
Further, we evaluate second interesting result in the form 
| (4.6) ред(5,2Ь)а.,-1/2-ой bu -1/2-03 Hy). (ho, hp, hay, h4) = 
| 2 (m ч (0), (ш +0), (b.), (Ha +о,), (8,0) (8,8) 
| = T (2b, -1)F (25, -1) n, (5, +1/2), (b, +1/2), nim! 
ind 


(2.— 25,,1),(1/2,1/2) 
(0,1),(1/2+m,1/2) 


А (ha) 


ME 


нэ! (1- A) (0,1),(1/2+m,1/2) 


provided that Re(b,)>-1/2,Re(b,) >-1/2,|ya <1,|h,5|< ere (f - 2) «л, and 


| ЕЕЕ 


Here, the Fox’s H-function [5] (Also see, Srivastava and Manocha [11,p.49]) is 


given by 
| (о. А,) (o. | ag Tir(5,- 25] Ir(.-; *45)7 : 
| m,n OBA EI I TD RC e oj j=l ] 
| (4.7) Н f "i 211 --i* I rep, BE) (8,-44) 


The integral in (4.7) converges absolutely if 
n p 
| m ; 0 
о=Ўв.- Ув, ATA 
(4.8) » : 2 | » о 
| and analytic in the sector 
(49) larg(z) «1/20m. ; 
$ е h-0 and 54-0 in (3.1), we get Оно os 
Srivastava and Manocha [11, p.148, eqn. (89), tor p 


= (À a (0/5, - n. 9,5.) hB, hay, h 6)Т" = 
(4.10) 50), цолугьдйггнеэ (ро, RB, Fas" 
-0 У 


valent to the result due to 
5=1] 
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M ^g tero d. aer] 
=(1 T) » п\(с,), (с), mem 


hjoT 
Е| +0,+7,0,,A+n;c, + n;h46, aT um «1. 


When, 5,—0, by direct computation of integral of (1.4) of from (2.1), we find that | 
(411) po^) (7 o hB, hyy,h,S) = 


1» 


+0,,b,; 
a^ 22 : AB [us eost tihi] 


When, 5, =0, and 5,—0, by direct computation of integral of (1.5) or from (2.1),we ; 
find that 


(4.12) p? h)enc (Аа, hb, hyy, h,8) = F, [bs *0O3,H, +0, 121753 h,y,h48] g 


When, b, =0, and 5,=0, by direct computation of integral of (1.6) or from (2.1), we 
find that 


5 г +6,,b,; ‚045 
o as (аһа) | di np n [etos : 21 
v C5; 


When, 5,—0,5,—0, and b4—0, by direct computation of integral of (1.7) or from 


(2.1), we find 


с; 


Setting 6, =0, b,=0, 55-0, b,-0, T=-S and i= =n in both sides of (3.1) and then, 
in left hand side of that on making use of (2.1), 


4 ‘ particularly, we find the generating 
relation due to Meixner [7] (Also, see Erdélyi [2, eqn. (8), р. 841) Е 


en H t0,,-n; ш T0,,-n; 
(4.15) | p np |е M һәја" 


= +з 


n=0 


(4.14) ЕЁ“ РОН, (ho. hB, hey, ho) „Е Р x Wu г 


n fa +в), (и, +о,), hh, BSS E Ш +03 t n,n—m; h BS 
n, ni(c,), (c), (1-8) es €, +n; 1-5 
E p +0, +n, n-n; h, ôS 
221 чр) 
€, +n; 1+S j 


A : n Е d 23 Ne = m = -1,Г = -S and x 25 -N in both sides of (3.6) 
and then, ш left hand side of that on makin : the 
generating relation due to Meixner [7] 8 use of (2.1), particularly, we find 


(Also, see Erdélvi 
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en И *0,,n—T|; l4 +0,,n-7; 
(4.16) Э | АВ | а "уа 


1? 23 


=(1+5)" УМ ыы [мны JE Е +0, MU 


n=o\ 70 п!(с,), (1+5)? с, +n; 1-8 


| T0, л, п-т; 
E 4:594 N; 7,6 | 


C, t n; 1-5 


Setting b, =0, b,—0, and b,=0 in both sides of (3.7), we find a result substantially 
t equivalent to Meixner’s result [7, eqn. (19b), p. 345] 


= (A), [de +0,7-; +0,,A +n; 
(4.17) Ураны! B 7212 ү "күн 


с; 23 
ES T 
(1-T) AES t 05,13 1551 
5. Conclusions and Discussions. From the results obtained in the 


previous sections, we conclude that if the sequence of quadruple hypergeometric 
function K,,(.) due to Exton [3,4] is taken in the form of the sequence 


(5.1) (К, (а,а,аза,ь (п) Б, (8) (пзе (н), (в), (8) (n) Put) 
п = 0,1,2,... 

where, n is independent of x,y,z,t,/,,h.h3, and h,,and b; (0) = 5, Vi 2 1,2,3,4, and 

c; (0) 2 c, Vj =1,2. 

H(D) is the compact set of quadruple hypergeometric functions 

(5.2) К„(а,а,а,а,Ь, (п),Ь, (л),Ь, (л), (п)зе (п), (72) 60 (п),с, (п) ња oo Put) » 

when n = 0,1,2,... 


D is the region given by 


; t 
<2<ү,0<1<5 and 085424 54 - 1. 
(53) 0<х<0,0<7у<[,0<2<ү,057 = ам y © 


Then, there exists a linear degeneration map 

PH CY (ha, haB, hy, h8): H (D) ә Rf or max {|0158} < 1,max (Бү 8|] <1, 
defined by 

(54) P^* ЇЕ, (a,a,a,a,b, (n),b, (n),5, (n)b, (n);¢ (п),с, (п), (n), (1) us hay 12223) 
(ha, hoB, hay, ha 5) = р" b m) (m). (п) lahele ( ha, hf, hy, һб) = 


n); C; t); 131, 9 : 
2 Ен, ехо Б (n), (se (n) Pa [Hs + Sorte Sata (0/44 halik : | 
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for а n=0,1,2,.. 
When, n=0, (5. 4) becomes equivalent to (2.1) such that 


(5.5) p^ [К (ossa, bn Pye Ca Њо ЊУ Ља) (ноа) = (E 
poo „babs bi) (ее) (ha, h. B, hy, h 49) = | 
| 
F. З +0, H + 65,0, ,023¢,;2,0, h,p]F. 3 + Оз, + Oo, ,5,,b,5c,; h зҮ, б]. 
Again, then for any complex апа [r| «1, (5.4) generates various generating | | 
relations and transformation formulae by the formula | 
(А) Т" s | 
(5.6) TUNE Р" [a,a,a,a,b, (п),Ь, (п),Ь, (п)Ь, (п); с, (п), 
К, (с, (п), с, (п), с, (n); а, уһ) (hia, h B, hY, hò) Ч 
< A) T” а, л),б5(1),03(л njc An, nJ)Ht,o [ 
- en шй) ее" Ср о р В, Љу, 7,5) | 0 
[ 
Example 1. When b, (n) = 5,5, (n) = b, = n,b, (n) = b,,b, (n) = b, – тС) (п) = апд 
с, (п) = с,, then making an appeal to (3.1) and (5.6), we find the generating relation [i 
SA) T” pue 
ч) » n! P Es (a,a,a,a,b,,b, — 1, b,,b, —n,0,6, 165505 x; 22223) | 
[i 
(hyo hytsh,B) = т) $ Atto), (He +04), | nut бат 3 
n=0 т\(с,) (с), Gan); [ 
F. [ 
— ОМ c +n; hon 7) | [ 
[ 
Е at +0; +}, + | 
9 зт б, t л, +0, nb, b, e тс, +п,с, +n, €, tn; 2211 | 
provided that |hya|<r,,|h,8|< ъ|<ъ- EAE «(вт 3 3 
ТҮ ЙЛ? such that r +r, = 1/2 | 
and |Р] <т,.]һ,8|<т, = Sla] |: T 
ЗГ Т 1 < T, such that T; + T, = гі. | | 


Ехашрїе 2. W = 
hen b,(n)=1+n,6, (n)=6,,b,(n)=2+n,b, (n) -bsa (n) =e» ап 


c, (n)=c,, then making an a 
formula рреа to (3.11) and (5.6), we find the transformation 
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SOT" petk | ) 

(5.8) 2, ai 2,0,0,0, +п,Ь,,). nbi ee, cs lx, hy iz, h T) 

(ho hP, Љу) = (1- 7) ^ X б), о), (14-64), [hp | pUsssnasauemas 
CIEGO oa] 


[Kul ombre та tme tme tme enl һу J ht | 


22 
-T UNT 
hB h,d ` 
| (h , 3 

(ro TEA eal 


} <1,|T|<1,max [Ау] «1. 
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ABSTRACT 
In this paper, we introduce product of generalized hypergeometric functions 
like I-function, Mittag-Leffler function and Bessel function by using Marichev- 
Saigo-Maeda fractional integral operator and we find out some image formulae by 
employing generalized hypergeometric function and generalized fractional integral 
operator due to Marichev-Saigo-Maeda involving Appell’s function РЁз(.). 
| 2010 Mathematics Subject Classification: 33D60 
| Key words and Phrases: Marichev-Saigo-Maeda fractional integral operator, 
Fox’s H-function, I-function, Mittag-Liffler function, Generalized hypergeometric 
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1. Introduction. Let iio тє С( (п) > 0) and х»0, then the 


generalized fractional integral operator is defined as : 


-н Ї 2 2 
ma d er аа Еа: 
(1.1) (ge ye KES 1, С zu 30 : 


aH EIE (oat Е 
(1.2) 12:55:22) 2-2 Ї(с-3) Е [ноти 2s Чи ) ? 


where Р.) is Appell’s function (also see Erdélyi [2]) defined by 
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„ (0) (w) (0) (v), x" y" 
(1.3) F, (ju, o, 05052, y) E. Wh RE : m!n! (max (hs, 3 1) 5 | 
The generalized fractional integral operator defined through equations (1.1) ang 
(1.2) has been introduced by Marichev [8] and later extended and studied by Saigo 
and Maeda [11]. These operators are known as the Marichev-Saigo-Maeda operators 
and properties of these operators were studied by Saigo and Maeda [11]. These 
fractional integral operators have many interesting applications in various sub. 
fields of mathematical analysis. Kim et.al. [6] have demonstrated applications 
related to certain class of complex analytic functions. For instance ([7],[8],[9] and 
[14]) can also be referred for some basic results on fractional calculus. 

In this paper, following results due to Saigo and Maeda [8] will be frequently 
used: 


S T(p)r(p*u-nu-u-v)F(p*v'-m) penny | 
[rs veni 1 x)= x! n-u-j'-1 
aa б Yo) Г(р+о)г(р+т=- = ш) Г(р+п- -v) , € 


where R(n) > 0/9 (p) » max (0,9 (a. +0 –), 51067-07) and 


(5) (Inest) - CP- YEU p nir w)r (L7 pev 7m) „з 
à T(1-p)F(1-p*u*p'*v'-n)F(1-p*u-v) 


where R(n) »0;0 « 9t(p) < max [9t (-v),9t(1 +p’ =n), R(u +v- n)} 


2. I-Function. The J-function is the generalization of H-function of Mathai- 
Завсар (7) and which was introduced Бу Saxena [12], while solving a dual integral , 
equation involving sum of H-functions as kernel. This J-function is defined as: | 


аА) (aA; 
(22) IG) - 155, [x] - 177. Sentus, -.[x(s)a"ds 
(8B) (5.85)... 2mi L 


Ji 
where 


(6 Bis] с-а — A. 


2 П Г(1-5,-8,8| TI 2273 


isl 2 
J=n+1 


with all other conditions already detailed in [12] 


The series representation of I-function is as follows : 


1 
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22522297 -SDO 
(fF a Ea) an Ез VIE, K 


(9.3) I(x)9I75,|2]- I5 5| 2 


ind 
igo where &-(e, -1—-v)/E,, h=1,2,...,m and 
ors 
p " ИШ? «F&JL]r (1-ё,-8д) 
ub. Oa o — 
x Ўта 8) е, 
ind 1-1 | Ј=т+1 Ј=п+1 
2. Marichev-Saigo-Maeda Operatoional Results Involving the 
Шу Product of two I-Functions. 
Theorem 1. Let x>0,%(p)20, and the parameters и,и,о,У,л,р e C,b,c e R such 
that 9t(B) > 0,9 (8) > 0,9t (n) > 0,9t (p) > max (0,9 (i *v-n),9t (i —v)], 
9(1+1/(1-о))>0 and i=1,2,...r. Then 
| Ё (eE, ie 12 Xo. тл f (854, ), n (аА, js E 
ах Тал BY (ъв 
: ( pF; „м? бий на ( J? Dm й” A I mela, 
аг (и), (v) (v) M,N 5 0:07:02:277 
al voie $ CO GL QD CP ма | (fF 
js 0420 Ril! (22), „(2 n 


I" л bx" 


Pi dii 


(1-р+1+ * 85 ),(аА,), „(ел Ал), | 
(;,B,), , (bi Ba) (<р ч +056) : 


? jie jt 


(2.1) and (2.3), we get 


(a; ? A, Ihe (а, 2 А; ВЭ | 


Proof. Making an appeal to (1.1), (1.3), 


pon en OS jee bt? 
a Рея) ла 


: : ^J (о), 097, (1-t/x) d 
тв - [e (к peu : - 
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у J DEN) LOCUM 22 
гаар ШЕ, ctn 2ni 7 
Now changing the order of integration and summation, under the given conditions, 
we derive 
еЕ,) v (e Es), 


пол | 0-1 PMN 5 
Ty: t In Qn c Bi dir 


лам m 
b 


pF; "OUS M+1,Q; 


(2), 0), (0), (9), SI C1) X (9) OBL го(8)ь- 
h=1 v=0 2ni +! 


Ep (т), AU! v! E,c* 


(imme s - n" ats) (1-а11) at Jas 


E ye 2). (09, (0) (9, r(n 1) Bg C1) XE) (4) 


kl=0 2 (1 ) mani hal 0-0 0!Е, 
©(s)(bx ) = Г(р-1-17-88-05) 
mt (5)( тя y -6E- Ba) 


Hence we get the required result (8.1). 
Special Cases. 


(i) Putting l=0,k= 0,n=0,n'=0,v= 0,0 =0,n = 1+n/(1- a) and reple 


xbyx/(a(1-a)) in the equation (3. D, we get the following result for ]-fanction: 


(e, E; ys; (s; E; a Р, 


1 8 (6 ша 


Пе bt! (2,А;),, 10:22:22 
22809 67 › 5, Bj Jean 


Pdr 


(8:2) Bi) лун СЯ 


j ls (eB, яр 


(f, 3 0) 5 112277 


“гар ee) eT 
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(1-р +55,8)(а,,4,), 5 485545) 


p 
Т" bx n*l,p, 


pi+l,qi+1;r ГЭЛОН 
[a(1-a)] (5, 9), 18. D (ox D | 
В В ; Ес, 


т 


which appears to be new. 


(ii) For р; = p,q; =q,P, = P,Q, =Q (i=1,2,3,...) and r=1 the result (3.2) reduces 
to well known result due to Srivastava, Lin and Wang [13, p.163, eqn (12)] for H- 


function: 
pho PN c (e, .E,) Нт^ bt? (a,,A,) 
) (6,В,) 


(8.3) 


1222 


bx? (1—р + б®,р).(а,,А,),....(а„, А„) 
HON Бэ 75-11) 
| «(1-0 (b, B. (ho Bo) (b 6 e) E 


4. Mittag-Leffler Function. According to Erdélyi et al. [3] (also see 
Prabhaker [10]) the Swedish Mathematician Gosta Mittag-Leffler introduced the 


function E,(z) defined by 


2 z" 


ёл) BeAT » 


1-0 


where z is a complex variable, Г(.) is a gamma function and B20. The Mittag- 


Leffler function is the direct generalization of the exponential function to which it 


B<1, itis interpolated between the pure exponential and 


reduces for B=1. For 0< 
Leffler function naturally occurs as 


a hypergeometric function 1/(1-z). Mittag- 
integral equations. 
Wiman [14] studied the generalization of E, (г), that is given by 


;p €C,R(p)>0,R(B)>9 


: k 
241732 
(4.2) Ep9(2)= 2, Г(ВЁ+р) 
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which is known as Wiman’s function. 


Prabhakar [10] investigated the function Ej, (z) as 


2 y k à 
(4.3) E;,, (2) = iyu" € C,R(p) 2 0, R(B) > 0, R(y) > 0 E 


The Mittag-Leffler type functions belong to family of /-functions, since they can 
be represented in terms of H -function: 


en ГИ 
11 шим : 


(44)  Ej(z)-H i 


(45) Е (2)= вм) 08: Ли єС,Е(р)»0, 


(С) =) | = E » 
«6 5.0) 2 ЧОЙ" 


1 ji ГѕГ(у-ѕ) 


E mily 27 Г(8-ав) (2) ";B,p,y€C,R(B) > 0 C 


T i i ight? : 
he relation connecting the Wright's function pW, (2) ((15],[16]) and the H-function 


was given for the first time in the monograph of Mathai and Saxena [7, p.11, eqn. j 


1.7.81] as: 
E - Le - 


5. Marichev-Saigo- 
Product of Lf ES Saigo-Maeda Operational Results Involving the 
‘function and Mittag-Leffler Function 


Th 
eorem 2. Let х» 0,9t(p) 20, and the parameters u 


(4.7) “| 
| ( 


b, B,)...(b,,B,) 


(1-2,,4,),...(1-a,, A,) 
(0,1),(1-5,,B,),...(1-55,Bs) | 


be W’, V, VN, p, y EC and b,c e? 
such th 
at Pr)» 0912) > 0,9%(1+n/(1-a)) > max 0, R(p)], 9t(p) > тах{® 


9 Г x ГД , 
(n+p + п), (ш -о)},8>0,8>0 and і=1,9,.. ғ. Th 
,2,...г. Then 


(5.1) Иса о) uie Es), 
CC-0. Gufukul Хө Ур Лий ра. п Extr), tive 


hasa 
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(e; E; Js 1027 El 


l ym u-p'+p- (хе ) e ш), (о), (0), Jut cu 


7a PITT P.Q;R 
an 
| (0-1-17-66,8),(1,1) p 
MU iia e | 


Proof. Making an appeal to (1.1), (1.8), (2.8) and (4.6), we get 
(6) к 102577 
(P) 20:12 


129Р:02077 


pu Um РТМ 


Ет, (82) 


5 
НО ct 


L! 


Suc =a 2-629) ҮЛ (зза (1-x/t) 


kI West 
1-y,1 
1212) fe 


(3838 (2) Р “ээн 


Now changing the order of integration and summation, under the given conditions, 


| we derive 
ion 
pE) „ев Ёк) 

qn. Гул со ct? (e; Ji ( а Ej, (bt*) 

{ (f, Е) (б) А 

(9,0, (9, eee XE) |_1_ рта) 
-A T 23 01ЕЁ,с" 547 pim r(B-os) 

һе te (x-4) (1-4/ху (L-x/t)' dias, 
Г БОРОО 55 УСЫ 
| TUN A (nH MC) 
0, 


s T(p-1-p-85-B5) ш, 


1 [MUT 


2nily © r(B-os) rT(p*n*k-E- 58-08) 
ЭЛ see Eius 
Dose (HO), OO) мд ЧЭ дЕ) ал 
"T È k'l! ENa (12:25:53 


r(Y) 


3-0 
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(1-р+1+ u 85 B), (1-7,1) 
(0,1),(1-0,B) (1-9 ^n - & +p’ + 865p) | 
Hence we get the required result (5.1). 

Special Cases. 


ни E 


(i) Putting 1=0,k=0,n=0,p' =0,v=0,v' 20,1214 n(1-a) and replacing x 


byx/(a(1-«)) in relation (5.1), we get the following result for I function: 


e;,E,) (ej, E; р 
CYPR ("Ign „|“ | ү Ej, (60?) ES — (т) 
Pog US us Y[a(1-o)] l-a 


bx! 


M, (e E;), (Es), (p-88,B),(y 4) ү 
B 
-8ь8 (40-01 


; йл 2V2 
30-31 =], Hm (В) M ү (р ‚В), Шил 


(ii) Now for p, = p,q; 4,8 = P,Q, =Q (i=1,2,3,...) and r=1 the result (5.2) reduces 
to well known result due to Srivastava, Lin and Wang [13, p.165, eqn 21] for H- 


function: 
(e, E,) 
(222) Tale) ў 


(5.3) ro any E 


(0-88 8),(y) E 
БӨ өнг -8а 8-0) 


АР 
6. B i 
essel Function. The Besse] function of the first kind J (x) is defined 
for complex xe C,x «0 and veC,R( ) 


о)»-1 by 
(64) 4,(5)- E 3l КЕЗЕР 
аа (22) 
and its H-function : 
bee Cone 51635814 Miwar. An eGangotri Initiative 


ed 
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{ 
Í 


ea 9-5) r oncon eroa | 


For the asymptotic behavior and other properties in detail see for instance ([1],[5] 
and [6]). 

7. Relations of the Product of I-Function and Bessel Functions 
using Marichev-Saigo-Maeda Operational Results Involving the Product 
of I-Function and Bessel Functions. 


Theorem 3. Let х» 0 and the parameters и,и,0,0,1,ү,0є C and ceR be such 


that R(n) > 0,9 (у+о) > 0,9 (12-m/(1- a)) » 0,8» 0,8» 0 and i=1,2....,r. Then 


PE pim t 40297:02277 
(np o Gee | pee 45) J, (t) 
: 2) 18485818 2 (f, Ej), Ел), о 


саг 
1V2 


(ej E;), м 22277 
(78) м 17227 


Quo |" 210 E Ч Ril! 2 
x 
(10,1), (y 80:51 *k-y -852) 4 
Proof. Making an appeal to (1.1), (1.3), (2.3) and (6.2), we derive 


: 1 (en E; 
os B YN. 45 а! (e; E; dene (е ёдэ p (t) 
E 2 


uo Ёс (7955) ад 
-p Fi 7, Ў о), 1-4/xy (1-х/ ) 
DAT „узе но la 2 


kl=0 "n k+1 


2 


pus (s Jo] 


(nepos 


(п), „810 =: Ес 


оор 


“гу 
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2ni +T(1+a-s 


grat To ^ 1: 7 ! E 55 


ни 


Thus we establish the required result (7.1). 
Special Cases. 


(e; Ej), y se E; E 


(f, ‚Е, Nee 175297 


(l-y-a+1+p' + 86,2) | 
4 |(0,1),(-0,1),(1-y-a-n-k +p’ + 96,2) 


(i) Putting /=0;k=0,n=0,p'=0,v=0,v'=0,n=1+n/(1-a) and replacing x by 


x/(a(1-a)) in the equation (7.1), we get the following result for I-function: 


| (eE;), (e E;),. 


2 (ү? ЫГ 
(1.2) Tek ) В Devote 48 


Ja (t) 
(755), 122277 
ИШ 0) a cs? (6), y 102297 


—————— 


2. [а (1 5s a) | Ң.9;8 [2а (1 -о)] 


122 1229 M 


du M 
1T2 

(1+0,1),(1+y+0+n/(1-a)-3¢,2) 
which appears to the new 

(ii) For Di = р,д; =4q,P. = P,Q; =Q @= 


1,2,3,...) and r=1 the result (7.2) reduces 10 
well known r 5 : = 
esult due to Srivastava, Lin and Wang [13 p.166 eqn (24)] for H-function: 


(e, E») 
[к Б, ie) 
r(1 + п/(1-о))хт=п 


== М 5 
2" а(1-а)| "54 Пээ 
СС-0. Gurukul Kangri Gollection 1460), 


4а? (1 = a) > 


(7.3) Pla) G НАА | c 


(е,, E, | 
(for Bed otri Initiative 


[1 


[1 


[15] 


116) 
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ү (y+a-86,2) =x" 
192 (1+0,1),(1+y+a+n/(1-a)—8€,2) 4a’ (1-a) | 
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